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ABSTRACT 


The  Tropical  Cyclone  Vortex  Tracking  Program  is  used  to  identify  vortices  in  the 
western  North  Pacific  from  the  Navy  Operational  Global  Atmospheric  Prediction  System 
(NOGAPS)  analyses  and  forecasts  during  May  -  October  2002  and  2003.  Based  on  the 
NOGAPS  analyses,  several  parameters  are  different  between  the  23  vortices  that 
developed  into  storms  during  2002  according  to  the  Joint  Typhoon  Warning  Center 
(JTWC)  and  the  231  vortices  that  did  not  develop.  After  eliminating  127  vortices  that  did 
not  persist  at  least  24  h,  this  left  104  non- developing  cases.  For  the  developing 
circulations,  the  average  850- mb  relative  vorticity  value  at  the  first  JTWC- warning  time 
was  5.0  X  10'^  s"\  with  an  easterly  deep  layer  wind  shear  of  -1.8  m  s The  average  850- 
mb  relative  vorticity  maximum  for  the  non- developing  cases  was  3.3  x  10'^  s'\  with  a 
westerly  vertical  shear  of  4.1  m  s The  NOGAPS  model  tends  to  over- forecast  relative 
vorticity  prior  to  formation  time  for  both  developers  and  non- developers.  Especially  for 
the  72-h  and  96-h  forecasts,  the  over- forecasting  tendency  leads  to  non- developing 
vortices  meeting  the  threshold  vorticity  value  of  the  developing  vortices.  The  tendency 
for  NOGAPS  to  forecast  the  non- developing  deep  layer  wind  shear  to  become 
increasingly  easterly  with  time  is  considered  to  be  a  major  factor  in  these  over- forecasts 
of  formation.  Some  adjustments  in  the  cumulus  parameterization  heating  and  moistening 
plus  convective  momentum  transport  may  improve  these  forecasts  of  tropical  cyclone 
formation. 
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represents  analyzed  average  850-mb  relative  vorticity  for  non- developing 
vortices  (at  -tOO),  and  red  line  represents  average  forecast  850-mb  relative 
vorticity  for  the  indicated  forecast  time.  Average  analyzed  850-mb 
relative  vorticity  for  this  subset  of  non- developing  vortices  at  bb  (light 
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(at  -tOO),  and  red  line  represents  average  forecast  deep  wind  shear  for  the 
indicated  forecast  time.  Average  analyzed  deep  wind  shear  at  N3  (light 

dashed  black  line)  was  4.1  m  s'^  Vo  (light  solid  black  line)  is  -1.8  m  s'^ . 61 

Figure  3.24.  Average  analyzed  and  forecast  sea  level  pressure  (SEP)  (mb)  for  non¬ 
developing  vortices  relative  to  No,  as  in  Eigure  3.18.  Heavy  black  line 
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Eigure  3.25.  Average  analyzed  and  forecast  925-mb  wind  speed  (m  s'^)  for  non¬ 
developing  vortices  relative  to  No,  as  in  Eigure  3.18.  Heavy  black  line 
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Eigure  3.26.  Average  analyzed  and  forecast  vapor  pressure  (500-700  mb  average.  Pa)  for 
non- developing  vortices  relative  to  N3,  as  in  Eigure  3.18.  Heavy  black 
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for  the  indicated  forecast  time.  Average  analyzed  vapor  pressure  at  No 
(light  dashed  black  line)  was  3.4  Pa.  E)  vapor  pressure  threshold  (light 
solid  black  line)  is  3.6  Pa. 
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Figure  3.27.  Example  of  forecast  tracks  for  an  over- forecast  non- developing  vortex. 

Circles  indicate  successive  12  h  forecasts  of  the  circulation  center,  with 
the  size  of  the  circle  indicating  the  intensity.  Black  circles  indicate 

analyzed  values,  and  colored  circles  indicate  separate  forecast  times . 69 

Figure  3.28.  Average  analyzed  and  forecast  850- mb  relative  vorticity  (x  10'^  s'^)  for 
over- forecast  non-developing  vortices  as  in  Figure  3.18.  Heavy  black  line 
represents  analyzed  average  850- mb  relative  vorticity  for  over- forecast 
non- developing  vortices  (at  -tOO),  and  red  line  represents  average  forecast 
850-mb  relative  vorticity  for  the  indicated  forecast  time.  Average 
analyzed  850-mb  relative  vorticity  for  this  subset  of  non- developing 
vortices  at  No  (light  dashed  black  line)  was  2.4  x  10'^  s ?o  (light  solid 

black  line)  is  5.0  x  10'^  s  ' . 70 

Figure  3.29.  Average  analyzed  and  forecast  deep  wind  shear  (200-850  mb)  for  non¬ 
developing  over-forecast  vortices  relative  to  Nq.  Heavy  black  line 
represents  analyzed  average  deep  wind  shear  for  non- developing  over¬ 
forecast  vortices  (at  -tOO),  and  red  line  represents  average  forecast  deep 
wind  shear  for  the  indicated  forecast  time.  Average  analyzed  deep  wind 
shear  at  No  (light  dashed  black  line)  was  3.2  m  s  '.  Vo  (light  solid  black 
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Figure  3.30.  Example  of  forecast  tracks  for  a  Group  2  non- developing  vortex.  Circles 
indicate  successive  12  h  forecasts  of  the  circulation  center,  with  the  size  of 
the  circle  indicating  the  intensity.  Black  circles  indicate  analyzed  values, 
and  colored  circles  indicate  separate  forecast  times . 73 


Figure  3.31.  Average  analyzed  and  forecast  850-mb  relative  vorticity  as  in  Figure  3.18, 
except  for  the  Group  2  subset  of  NDl  vortices.  Heavy  black  line 
represents  average  analyzed  850-mb  relative  vorticity  for  the  Group  2 
subset  of  vortices,  and  heavy  red  line  represents  average  forecast  850-mb 
relative  vorticity  for  the  time  indicated.  The  light  solid  line  represents 
average  analyzed  850-mb  relative  vorticity  at  FO  for  developing  vortices 
(5.0  X  10-5  s-1),  and  the  light  dashed  line  represents  average  analyzed 

850-mb  relative  vorticity  at  NO  (3.7  x  10-5  s-1) . 74 

Figure  3.32.  Average  analyzed  and  forecast  deep  layer  wind  shear  (200-850  mb)  (m  s  ') 
for  non- developing  vortices  that  were  forecast  to  exceed  F),  as  in  Figure 
3.18.  Heavy  black  line  represents  analyzed  average  deep  layer  wind  shear 
for  this  subset  of  non- developing  vortices  (at  -tOO),  and  red  line  represents 
average  forecast  deep  layer  wind  shear  for  the  indicated  forecast  time. 
Average  analyzed  deep  layer  wind  shear  850  mb  relative  vorticity  for  this 
subset  of  non- developing  vortices  at  Nq  (light  dashed  black  line)  was  2.8 

m  s  '.  Vo  (light  solid  black  line)  is  -1.8  m  s  ' . 75 

Figure  3.33.  Focations  of  NOGAPS  model  false  alarms  in  the  western  North  Pacific 
during  1  May  -  31  October  2002.  Circles  indicate  successive  12-h 
analyses  of  each  identified  model  false  alarm.  The  size  of  each  circle  is 
related  to  the  analyzed  intensity  of  the  circulation  with  a  larger  circle 
indicating  a  more  intense  circulation. 
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Locations  of  JTWC  TCFA  false  alarms  in  the  western  North  Pacific  during 
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1.  INTRODUCTION 


Recent  improvements  in  operational  global  models  have  lead  to  improvements  in 
tropical  cyclone  track  forecasting.  Accurate  72-h  track  forecasts  are  now  common,  and 
120-h  forecasts  have  been  issued  beginning  in  2003  because  of  their  utility  in  military 
and  civilian  operations.  The  requirement  to  accurately  forecast  to  120  h  makes  it 
necessary  that  the  models  be  able  to  forecast  the  entire  lifecycle  of  a  tropical  storm  within 
a  single  numerical  model  run.  These  global  models  must  be  able  to  forecast  tropical 
cyclone  formation,  because  a  cyclone  may  form  and  intensify  to  a  damaging  storm  within 
the  120-h  forecast  period.  Accurate  tropical  cyclone  formation  forecasts  wjuld  greatly 
enhance  the  ability  of  the  U.S.  Armed  Forces  to  evaluate,  plan,  and  execute  any  actions 
necessary  to  minimize  damage  or  loss  of  assets  or  personnel. 

An  assessment  of  the  current  capabilities  of  operational  global  (dynamical) 
models  is  thus  necessary  to  determine  the  potential  for  an  accurate  tropical  cyclone 
formation  forecast  product.  Such  a  product  would  provide  a  useful  tool  for  determining 
which  tropical  cyclone  precursors  have  the  most  predictive  value  in  discriminating 
between  developing  and  non- developing  circulations.  The  first  step  towards  developing  a 
tropical  cyclone  formation  product  is  to  analyze  the  atmospheric  variables  at  the  time  of 
formation. 

A.  FORMATION  DEFINITION 

The  lack  of  a  single  unambiguous  definition  of  tropical  cyclone  formation 
provides  a  significant  challenge  to  both  operational  forecasters  and  researchers. 
Researchers  may  cite  formation  as  occurring  when  an  organized  rotary  circulation 
develops  a  warm  core  thermal  structure  and  tangential  winds  decreasing  with  height  (i.e., 
Enagonio  and  Montgomery  2001),  or  after  the  development  of  “a  mesoscale,  warm  core 
vortex  capable  of  self- amplification  through  air-sea  interaction”  (Davis  and  Bosart  2003). 
Both  of  these  definitions  lack  the  necessary  objectivity  required  for  use  in  operations,  and 
require  significant  subjective  interpretation  by  the  forecaster  or  analyst. 

Operationally,  the  Joint  Typhoon  Warning  Center  (JTWC)  designated  a 
circulation  as  a  tropical  cyclone  when  a  specific  tangential  wind  speed  threshold  is 
reached  ^  25  knots  (kt)  in  the  western  North  Pacific).  Elsberry  (2003)  adds  more 
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stringent  requirements,  requiring  “a  non- frontal,  cyclonic  circulation  in  the  tropics  that  is 
closed  (ground  relative  westerly  on  the  equatorward  side)  with  maximum  10- minute 
averaged  surface  sustained  winds  of  at  least  25  kt  (12  m  s"^)  that  is  accompanied  by  deep 
(throughout  most  of  the  troposphere)  convection  and  a  radius  of  maximum  winds  such 
that  the  Rossby  number  is  at  least  one.”  The  formation  definition  used  here  follows 
Elsberry  (2003),  with  the  additional  requirements  of  a  closed  outer  contour  of  850- mb 
relative  vorticity,  and  a  minimum  duration  or  persistence  in  the  model  analysis  of  at  least 
24  hours.  Due  to  the  ambiguity  of  determining  the  actual  formation  time  for  each 
developing  cyclone,  the  specific  formation  time  for  a  numbered  circulation  is  assumed  to 
be  the  time  of  the  first  warning  issued  by  JTWC. 

B.  TROPICAL  CYCLONE  FORMATION  MECHANISMS 

Tropical  cyclone  formation  has  long  been  associated  with  a  set  of  necessary  large- 

scale  environmental  conditions  (e.g..  Gray  1968).  These  conditions  can  be  summed  up  as 
“high  sea- surface  temperature  with  a  relatively  deep  ocean  mixed  layer,  large  values  of 
absolute  vorticity  in  the  lower  troposphere,  weak  vertical  wind  shear  over  a  pre-existing 
disturbance,  and  mean  upward  motion”  (Harr  et  al.  1996).  In  the  presence  of  such 
environmental  conditions,  a  disturbance  often  experiences  tropical  cyclone  formation. 
Scientists  understand  that  maintenance  and  enhancement  of  the  necessary  conditions  in 
the  large-scale  tropical  atmosphere  is  necessary,  yet  the  actual  mechanisms  that  cause  or 
lead  to  the  development  of  a  warm  core  remain  uncertain. 

Davis  and  Bosart  (2001)  cite  the  importance  of  a  pre-existing  large-scale 
disturbance  capable  of  organizing  convection  in  the  tropical  cyclone  formation  process. 
This  large-scale  disturbance  can  be  a  synoptic  baroclinic  zone,  the  monsoon  trough,  or  a 
monsoon  depression.  By  concentrating  potential  vorticity,  this  large-scale  disturbance 
allows  the  growth  of  the  incipient  tropical  seed  disturbance. 

The  process  by  which  the  seed  disturbance  becomes  a  coherent  mesoscale 
circulation  may  not  be  a  single  mechanism,  but  a  series  of  processes  or  events  (Davis  and 
Bosart  2001)  that  concentrate  mesoscale  vorticity.  Possible  processes  capable  of  this 
undertaking  that  have  been  explored  are  the  vortex  interaction  theories  of  Simpson  et  al. 
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(1997)  and  Harr  et  al.  (1996),  and  formation  from  a  wave  in  the  tropical  easterlies 
(Molinari  et  al.  2000). 

The  two  prevailing  formal  theories  of  tropical  cyclogenesis  are  conditional 
instability  of  the  second  kind  (CISK)  and  wind- induced  surface  heat  exchange  (WISHE) 
(Molinari  et  al.  2000).  Both  of  these  theories  can  explain  processes  during  the  mature 
stage  of  a  tropical  cyclone,  but  do  not  sufficiently  explain  the  interaction  between  the 
background  synoptic  environment  and  a  mesoscale  seedling  disturbance. 

Hypotheses  for  development  of  a  tropical  circulation  generally  include  some  sort 
of  interaction  between  mesoscale  and  synoptic -scale  systems.  While  synoptic-  and 
planetary- scale  features  can  be  diagnosed  using  relatively  coarse  satellite  imagery,  and 
global  model  analyzes,  the  inability  to  verify  these  hypotheses  provides  a  significant 
challenge  to  scientists  due  to  the  lack  of  mesoscale  observations.  Improved  spatial  and 
temporal  resolution  of  satellite  imagery  over  the  tropical  oceans  adds  value  to  the  forecast 
process.  However,  the  paucity  of  regular  in  situ  observations  continues  to  hinder 
verification  of  a  single  plausible  mechanism  for  tropical  cyclone  formation.  Given  the 
challenge  of  actually  verifying  formation  mechanism  hypotheses,  past  research  has  often 
focused  on  tropical  cyclone  formation  after  a  symmetric  low-  level  circulation  has  already 
been  established  (Davis  and  Bosart  2001). 

The  ability  of  global  models  to  accurately  represent  the  thermodynamic 
conditions  necessary  for  tropical  cyclone  formation,  i.e.,  high  sea- surface  temperature 
and  conditional  instability  with  a  moist  mid-troposphere,  s  simplified  by  the  lack  of 
substantial  daily  variation  of  these  variables.  The  global  numerical  models,  even  with 
relatively  coarse  spatial  resolution,  are  capable  of  accurately  resolving  these  slowly 
changing  features.  By  contrast,  the  dynamic  variables  associated  with  tropical  cyclone 
formation,  i.e.,  low-level  cyclonic  vorticity  and  upper- level  anticyclonic  vorticity  plus 
minimum  vertical  wind  shear,  vary  on  much  shorter  temporal  and  spatial  scales  and  are 
thus  more  challenging  to  accurately  resolve  in  the  global  numerical  models. 

The  current  predictability  limit  for  mid- latitude  synoptic -scale  circulations  is 
about  five  days.  In  the  more  slowly  varying  tropics,  predictability  of  thermodynamic 
features  may  exceed  the  five-day  mid- latitude  threshold.  However,  given  that  a  tropical 
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disturbance  can  develop,  intensify,  and  dissipate  in  less  than  five  days,  the  ability  of  the 
global  model  to  accurately  forecast  this  process  is  questioned. 

C.  TROPICAL  CYCLONE  FORMATION  IN  THE  WESTERN  NORTH 

PACIFIC 

Adding  another  significant  challenge  to  both  tropical  cyclone  forecasters  and 
modelers  in  the  western  North  Pacific  is  the  large  number  of  potential  seed  disturbances 
that  form  year-round  in  the  region.  Prior  to  extensive  analysis,  over  300  potential  seed 
disturbances  were  identified  in  the  Navy  Operational  Global  Atmospheric  Prediction 
System  (NOGAPS)  model  in  the  western  North  Pacific  during  2002.  Only  31  of  these 
potential  seedlings  developing  into  tropical  circulations  that  were  later  warned  on  by 
JTWC. 

In  the  western  North  Pacific,  monsoon  conditions  dominate  the  favored  tropical 
cyclone  formation  regions.  These  include  formation  within  and  along  the  periphery  of  the 
monsoon  trough,  and  from  monsoon  depressions  and  gyres.  Non- monsoonal  formations 
do  occur  north  of  the  monsoon  trough,  but  with  far  less  frequency  than  monsoon- type 
formations. 

I.  Monsoon  Trough 

A  significant  number  of  tropical  cyclone  formations  occur  in  the  western  North 
Pacific  Ocean  basin  within  the  monsoon  trough  (Elsberry  2003)  during  the  Northern 
Hemisphere  summer.  The  monsoon  trough  is  manifest  as  a  cyclonic  shear  zone  between 
the  equatorial  westerlies  and  trade-wind  easterlies  (Harr  et  al.  1996),  and  is  characterized 
by  high  sea- surface  temperatures,  a  conditionally  unstable  atmosphere  with  a  moist  mid¬ 
troposphere,  and  a  region  of  weak  vertical  wind  shear  between  strong  westerly  shear  to 
the  north  and  easterly  shear  to  the  south.  Within  the  favorable  background  cyclonic 
vorticity  of  the  synoptic- scale  monsoon  trough,  cloud  clusters  easily  develop.  These 
cloud  clusters  are  a  known  precursor  to  tropical  cyclone  formation  (Simpson  et  al.  1997), 
and  may  potentially  trigger  tropical  cyclone  formation  on  the  mesoscale. 

Low-level  convergence  in  the  monsoon  trough  produces  clusters  of  deep 
convection.  Latent  heat  release  within  the  deep  convection  helps  to  enhance  the  strength 
of  these  cloud  clusters,  and  may  allow  for  the  formation  of  embedded  mesoscale 


4 


convective  systems  (MCSs).  As  outlined  by  Elsberry  (2003),  mid- level  mesoscale 
cyclonic  vortices  (MCVs)  tend  to  spin  up  in  the  stratiform  rain  region  behind  the  deep 
convection  of  a  long-lived  MCS  as  a  result  of  vertical  stretching  of  the  atmospheric 
column  between  ascent  in  the  upper-  level  cloud  and  subsidence  below  the  cloud  resulting 
from  evaporation  of  the  stratiform  precipitation.  The  maximum  amplitude  of  the  MCV 
occurs  in  the  middle  troposphere  between  the  700  and  300  millibar  (mb)  levels  where 
vertical  stretching  is  maximized  between  lo\w  level  subsidence,  divergence  aloft,  and  the 
resulting  horizontal  convergence.  The  MCV  is  warm  core  above  the  level  of  maximum 
intensity,  and  cold  core  below  it.  Key  requirements  for  the  downward  translation  of  the 
mid- level  warm  core  vortex  include  a  localized  potential  vorticity  maximum  and  strong 
convergence  in  the  low  levels  below  the  vortex,  and  an  associated  atmospheric  warm 
anomaly.  The  actual  mechanism  for  translating  the  warm  core  vortex  from  the  mid¬ 
troposphere  down  to  the  surface  remains  uncertain  (Davis  and  Bosart  2001),  but  once  the 
warm  core  vortex  is  translated  down  to  the  surface,  it  can  tap  into  the  latent  heat  energy 
of  the  ocean  (Nielsen- Gammon  1996),  and  CISK  and  WISHE  may  then  processes  begin 
to  intensify  the  developing  circulation.  Harr  et  al.  (1996)  and  Simpson  et  al.  (1997) 
propose  that  MCS  interactions  within  the  monsoon  trough  can  lead  to  merger  of  MCVs, 
downward  extension  of  the  merged  circulation,  and  eventual  development  of  the  eye  and 
inner  rain  band  from  these  MCVs. 

2.  Eastern  Monsoon  Trough 

The  eastern  end  of  the  monsoon  trough  is  also  highly  favored  for  tropical  cyclone 
formation  (Elsberry  2003).  In  this  region,  the  confluence  of  tradewind  easterlies  and 
equatorial  westerlies  produces  favorable  low-level  convergence  and  cyclonic  vorticity, 
and  deep  convection  is  favored.  Ritchie  (1995)  and  Holland  (1995)  associate  the 
preferential  formation  in  this  region  as  a  result  of  a  wave -energy  accumulation  from 
westward- propagating  waves  in  the  easterly  flow.  The  actual  presence  of  low-level 
easterly  waves  in  the  western  North  Pacific  remains  uncertain  though  the  regression 
analysis  and  band-pass  filtering  studies  of  Sobel  and  Bretherton  (1999)  indicate  the 
presence  of  such  waves  in  the  easterly  flow.  With  no  other  obvious  synoptic- scale 
features  present  at  the  time  of  tropical  cyclone  formation,  mesoscale  features  likely  play  a 
substantial,  though  as  yet  unobserved,  role  in  formation. 
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D.  PLAN  FOR  THESIS 

The  2002  NOGAPS  model  analysis  and  forecast  fields  will  be  analyzed  for 
developing  and  non- developing  circulations  to  determine  threshold  values  for  specific 
forecast  variables  that  distinguish  the  developing  from  non- developing  vortices.  These 
thresholds  will  then  be  applied  to  define  false  alarms  during  2002.  The  false  alarm  is 
defined  when  the  NOGAPS  model  forecast  development  to  surpass  the  threshold  for 
circulations  that  were  not  eventually  numbered  or  warned  on  by  JTWC.  The  threshold 
values  will  also  be  applied  to  storms  that  developed  in  2003,  and  the  NOGAPS  forecasts 
of  non- developing  circulations  during  the  same  year. 

The  NOGAPS  model  data  and  the  analysis  procedures  will  be  described  in 
Chapter  II.  Developing  and  non- developing  circulations  identified  in  Chapter  II  will  be 
presented  in  Chapter  III.  Developing  circulations,  identified  as  those  numbered  by  the 
Joint  Typhoon  Warning  Center  (JTWC)  from  1  May  -  31  October  2002,  are  examined  as 
a  whole,  and  also  broken  into  several  subsets  to  highlight  potential  formation  indicators 
and  proxies  for  formation  among  several  selected  atmospheric  variables.  These  subsets 
include  formation  subregion  and  formation  time  relative  to  a  formation  time  window 
determined  for  each  storm  by  the  JTWC.  Non- developing  vortices  are  examined  in  a 
similar  fashion  to  the  developing  circulations.  Additionally,  several  subsets  of  the  non¬ 
developing  circulations  dataset  are  examined:  model  false  alarms  and  model  over¬ 
forecasts.  Developing  and  non- developing  circulations  from  2003  are  also  presented,  as 
well  as  two  case  studies  from  2002. 
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n.  METHODOLOGY 


A  large  number  of  analyzed  and  foreeast  fields  must  be  examined  to  assess  the 
potential  for  tropical  cyclone  formation.  The  technique  summarized  here  uses  the 
Tropical  Cyclone  Vortex  Tracking  Program  (TCVTP)  developed  by  Professor  Patrick 
Harr.  This  program,  which  is  summarized  in  Figure  2.1,  searches  model- analyzed  and 
forecast  vorticity  fields  to  define  circulation  centers  that  may  or  may  not  develop  into 
tropical  cyclones.  It  provides  an  objective  method  to  detect  vorticity  centers  at  the  850- 
mb  level  and  matches  forecast  circulation  centers  with  analyzed  circulation  centers  at  the 
verification  time.  For  each  circulation,  a  series  of  formation-related  environmental 
parameters  relative  to  each  circulation  center  is  extracted  throughout  its  life  cycle. 

A.  MODEL  DATA 

The  model  fields  used  in  the  analysis  (Table  2.1)  are  Navy  Operational  Global 
Atmospheric  Prediction  System  (NOGAPS)  version  4.0  (Hogan  and  McClune  2002) 
analyses  and  forecasts  at  12,  24,  36,  48,  60,  72,  84,  96,  and  120  hours.  Since  the  area  of 
interest  is  the  western  North  Pacific  Ocean,  the  spatial  domain  covers  105°E-  180°E,  and 
from  the  Equator  to  30°N.  Model  resolution  is  one  degree  latitude  and  longitude  and  the 
time  resolution  is  12  hours  (00  UTC  and  12  UTC)  for  the  analysis  fields  and  six  hours  for 
the  forecast  fields.  The  period  of  study  is  from  1  May  -  31  October  2002. 


Table  2.1.  NOGAPS  fields  used  in  TCVTP  analysis. 

850  mb  relative  vorticity  (10'^  s‘^) 

Shallow  layer  vertical  wind  shear  (500  -  850  mb)  (m  s'^) 
Deep  layer  vertical  wind  shear  (200  -  850  mb)  (m  s'^) 
Geopotential  height  thickness  (200  -  1000  mb)  (gpm) 
Shallow  (1000  -  500  mb)  warm  anomaly  (K) 

Surface  latent  heat  flux  (W  m"  ) 

Total  (convective  plus  grid  scale)  precipitation  (kg  m"^) 
Vertical  motion  (Pa  s  ') 

Vapor  pressure  (500  -  700  mb  average)  (Pa) 

Sea- level  pressure  (SEP)  (mb) 

925-mb  wind  speed  (m  s"^) 

700-mb  wind  speed  (m  s'^) 

500-mb  wind  speed  (m  s"^) 
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B.  ANALYSIS  PROCEDURES 

Figure  2.1  is  a  summary  of  the  steps  in  the  TCVTP  algorithm.  The  TCVTP  total 
system  is  designed  to  identify,  track,  catalog,  and  summarize  tropical  circulations  within 
numerical  model  analyzed  and  forecast  fields.  Figure  2.1a  is  a  summary  of  the  databases 
and  catalogs  generated  by  the  TCVTP  process,  and  Figure  2.1b  is  an  extension  of  the 
TCVTP  algorithm  that  extracts,  dentifies  and  catalogs  tropical  circulations  within  the 
model  analysis  and  forecast  fields.  Individual  steps  within  the  TCVTP  algorithm  are 
indicated  by  the  number  on  the  left  side  of  Figure  2.1b. 


Figure  2.1a.  Summary  of  steps  used  in  the  TCVTP  algorithm. 

1.  Analyzed  Circulation  Identification 

As  part  of  the  first  step  in  Figure  2.1b,  850- mb  relative  vorticity  from  each 
NOGAPS  initial  analysis  is  computed  to  identify  each  circulation.  All  relative  vorticity 
maxima  greater  than  1.5  x  10'^  s'^  are  identified  as  trackable  circulations.  In  previous 
research  by  Dorics  (2002)  for  the  North  Atlantic,  the  minimum  vorticity  threshold  was  set 
to  1.0  X  10'^  s The  threshold  of  1.5  x  10'^  s'^  is  used  for  western  North  Pacific  vortices 
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Figure  2.1b.  Expansbn  of  steps  contained  in  TCVTP  block  of  Figure  2.1a. 


due  to  the  more  cyclonic  monsoon  trough  environment.  The  850- mb  relative  vorticity 
field  for  the  current  model  analysis  is  examined  first.  An  ellipse  is  fit  to  each  trackable 
circulation  with  an  outer  closed  vorticity  contour  of  at  least  1.5  x  10'^  s‘*  (Figure  2.2). 
The  center  of  the  ellipse  is  defined  at  the  position  of  the  relative  vorticity  maximum.  The 
ellipse- fitting  routine  is  based  on  a  bivariate  normal  probability  distribution  and  spans  the 


9 


Figure  2.2.  Example  of  ellipses  (red  lines)  fit  to  an  analyzed  vorticity  field  (x  10'^  s‘^; 
solid  lines)  by  the  TCVTP.  The  large  dots  indicate  the  TCVTP- identified  circulation 
centers.  Shading  defines  cloud -top  temperatures  derived  from  infrared  satellite  imagery. 


0.95  probability  levels  of  the  distributions.  Thus,  the  size  of  the  circulation  is 
characterized  by  the  area  enclosed  by  the  ellipse.  Each  circulation  is  given  a  unique 
designator  related  to  the  time  and  location  and  is  then  saved  into  the  Analyzed  Database 
(right  box  in  Figure  2.1a).  This  vorticity  analysis  process  is  repeated  for  each  model 
forecast  time. 

2.  Identifying  Tracks  of  Analyzed  Circulations 

To  generate  a  vortex  track,  circulations  identified  in  the  current  model  analysis 
using  the  process  outlined  in  the  Step  1  are  compared  to  circulations  identified  in  the 
previous  12- hour  analysis.  This  comparison  forms  the  basis  of  step  2  of  Figure  2.1b,  and 
determines  if  a  recently  identified  circulation  can  be  matched  to  a  previously  identified 
circulation  stored  in  the  Analyzed  Database  of  the  TCVTP  (right  box  in  Figure  2.1a). 
This  directory  contains  all  previously  analyzed  circulations  that  are  currently  active.  The 
distance  and  direction  of  the  new  circulation  relative  to  each  previous  circulation  are  used 
to  match  circulations  in  the  current  analysis  with  previously  analyzed  circulations. 

The  distance  and  direction  criteria  vary  based  on  the  translation  speed  of  the 
analyzed  circulation  (Table  2.2).  If  the  prior  translation  speed  is  small  (less  than  5  kt). 
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the  allowable  direction  orientation  to  the  new  circulation  is  relaxed  to  a  circle  to  allow  for 


a  stalled  situation.  When  a  circulation  from  the  current  analysis  is  matched  to  a  pre¬ 
existing  circulation,  it  becomes  the  next  point  of  the  circulation  track.  If  the  circulation 
cannot  be  matched  to  a  pre-existing  circulation,  it  is  stored  as  a  new  vortex  to  be 
potentially  matched  in  the  next  analysis. 


Table  2.2.  Summary  of  speed  and  distance  criteria  used  by  the  TCVTP  algorithm  to 
potentially  match  forecast  circulations  to  analyzed  circulations. 


Speed 

Distance  change  in  6  h 

Direction  change 

Slow  (<  5  kt) 

<  50  km 

360° 

Medium  (6-15  kt) 

51-150  km 

+!-  90° 

Fast  (>  15  kt) 

>  150  km 

+!-  45° 

Information  used  to  identify  and  characterize  each  circulation  via  the  ellipse 
parameters  is  given  in  Table  2.3.  Analyzed  circulations  are  assigned  a  unique  identifier 
(wpcyyyymmddhh_ll_nnn),  where  wpc  designates  a  western  North  Pacific  circulation, 
yyyy  is  the  year,  mm  the  month,  dd  the  day,  hh  the  time  in  UTC,  11  the  initial  latitude,  and 
nnn  the  initial  longitude  corresponding  to  the  time  and  location  of  the  first  appearance  of 
the  circulation  in  the  NOGAPS  analysis  (-1-00).  Tracks  are  thus  identified  by  the 
designation  given  to  the  first  analyzed  circulation  position  in  the  series.  Each  circulation 
is  also  described  by  the  circulation  size  (number  of  grid  points  within  the  ellipse  fit  to  the 
outer  closed  vorticity  contour),  shape  (lengths  of  major  and  minor  axes,  angle  of  major 
axis  with  respect  to  north),  and  orientation  of  the  ellipse  with  respect  to  the  bivariate 
distribution  (Table  2.3).  These  data  are  used  to  regenerate  the  ellipse  during  post¬ 
analysis.  The  data  in  Table  2.3  are  then  used  to  generate  a  history  file  unique  to  each 
analyzed  circulation. 


Table  2.3.  Information  used  to  identify  each  circulation  and  characterize  the  ellipse  fitted 
to  each  analyzed  circulation  (after  Dorics  2002  Table  2.2). 


Vortex  identification  &  model  prognosis 

Vortex  (ellipse)  specification 

Name 

wpcyyyymmddhh_ll_nnn 

Size  (number  of  enclosed  grid  points) 

Model  DTG 

yyyymmddhh 

Shape  (ratio  of  major/minor  axes) 

Forecast  time  (tau) 

ttt 

Ellipse  major  axis  (km) 

Fatitude  (°N) 

11 

Ellipse  minor  axis  (km) 

Fongitude  (°E) 

nnn 

Ellipse  angle  (relative  to  North) 

Bivariate  distribution  correlation 
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In  addition  to  the  identification  and  ellipse  characteristics,  the  NOGAPS  fields 
listed  in  Table  2.1  are  used  to  calculate  average  values  over  the  ellipse  of  other 
atmospheric  variables  that  characterize  each  circulation.  As  listed  in  Table  2.4, 
maximum  or  minimum  and  the  average  value  for  each  variable  are  calculated  for  the 
entire  ellipse,  or  for  each  quadrant.  These  values  then  characterize  the  vortex 
environment  for  each  model  analysis  and  forecast  time  in  terms  of  variables  that  are 
commonly  associated  with  tropical  cyclone  formation.  As  a  circulation  is  tracked  in  the 
analyses  or  forecasts,  a  history  file  of  the  variables  in  Table  2.4  is  created,  with  one  line 
per  analysis  or  forecast  time.  Comparisons  of  the  analyzed  and  forecast  history  files  will 
form  the  basis  of  the  analysis  of  the  model  forecast  accuracy. 

3.  Forecast  Circulation  Identification 

The  same  ellipse- fitting  process  described  in  section  II. A.  1  is  applied  to  all  850- 
mb  relative  vorticity  forecast  fields  from  the  current  model  integration  (step  3  of  Figure 
2.1b).  An  ellipse  is  fit  to  each  forecast  relative  vorticity  maximum  that  meets  the 
threshold  aiterion  to  define  circulations  from  each  forecast  time  in  the  current  model 
integration.  All  variables  listed  in  Tables  2.2  and  2.4  are  assigned  to  forecast  circulations 
as  well. 

After  the  circulations  in  the  6- hour  forecast  fields  are  identified,  they  are  matched, 
if  possible,  to  circulations  that  were  identified  in  the  current  (-tOO)  analysis.  Similarly, 
circulations  tracked  in  the  12- hour  forecast  field  are  matched,  if  possible,  to  circulations 
in  the  6-hour  forecast  field.  This  process  continues  to  the  120-hour  forecast  and  defines 
the  track  of  each  forecast  vortex.  The  track  is  continued  until  a  vortex  can  no  longer  be 
identified  in  the  forecast  vorticity  fields.  In  the  vortex  history  file,  a  line  is  added  for 
each  forecast  field  (i.e.,  -i-06,  -i-12,  ... -1-120)  to  summarize  the  location,  ellipse  parameters, 
and  model  forecast  environmental  parameters. 

If  a  vortex  is  identified  in  a  forecast  field,  but  was  not  matched  with  a  previously 
analyzed  vortex,  it  is  classified  as  a  “forecast  vortex.”  Forecast  vortices  are  named  based 
on  the  DTG  of  the  model  integration,  location,  and  forecast  time  in  which  they  first 
appear.  Although  the  history  files  for  these  cases  are  written  as  defined  with  analyzed 
vortices,  the  files  are  stored  as  “unclaimed  forecasts”  in  the  Forecast  Database  in  Figure 
2.1a. 
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Table  2.4.  Values  tracked  by  the  TCVTP  for  each  atmospheric  variable. 


Quadrant 

Ellipse 

Variable 

Average 

Average 

Maximum 

Minimum 

Relative  vorticity  (10'^  s'^) 

X 

Shallow  layer  vertical  wind  shear 
(500  -  850  mb)  (m  s'^) 

X 

Deep  layer  vertical  wind  shear  (200 
-  850  mb)  (m  s'^) 

X 

Geopotential  height  thickness  (200 
-  850  mb)  (gpm) 

X 

X 

X 

Warm  anomaly  (1000  -  500  mb) 
(K) 

X 

X 

X 

Surface  latent  heat  flux  (W  m'"^) 

X 

X 

X 

Total  precipitation  (kg  m‘^) 

X 

X 

X 

Vertical  motion  (Pa  s'^) 

X 

X 

X 

Vapor  pressure  (500  -  700  mb 
average)  (Pa) 

X 

X 

X 

Sea- level  pressure  (SLP)  (mb) 

X 

X 

X 

925  mb  wind  speed  (m  s'^) 

X 

X 

700  mb  wind  speed  (m  s'^) 

X 

X 

500  mb  wind  speed  (m  s'^) 

X 

X 

4.  Linking  Forecast  Circulations  from  Previous  Model  Integrations  with 
Analyzed  Circulations  from  the  Current  Model  Integration 

The  unclaimed  forecasts  catalog  in  the  Forecast  Database  (left  box  of  Figure  2.1a) 
contains  all  forecasts  not  yet  matched  to  an  analyzed  circulation.  In  step  4  of  Figure  2.1b, 
all  unmatched  forecast  circulations  from  previous  model  integrations  are  next  compared 
to  the  analyzed  circulations  from  the  current  model  integration.  When  the  matching 
criteria  are  first  met,  that  unclaimed  forecast  circulation  is  attached  to  the  analyzed 
circulation  and  becomes  the  formation  forecast  for  the  analyzed  circulation. 

5.  Finalization  of  Tracks  and  Forecasts 

At  the  completion  of  the  TCVTP  process  shown  in  Figure  2.1b,  only  two 
outcomes  are  possible:  finalization  of  a  tracked  circulation,  or  failure  to  successfully 
match  a  forecast  circulation  with  an  analyzed  circulation.  A  tracked  circulation  is 
finalized  when  no  subsequent  model  analyses  contain  a  circulation  that  can  be  matched 
with  an  existing  track.  When  this  occurs,  the  track  is  finalized  and  stored  in  the  Final 
Circulation  Catalog  (bottom  right  box  of  Figure  2.1a).  The  circulations  stored  in  the  final 
circulation  directory  represent  the  data  available  for  further  analysis.  When  a  forecast 
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circulation  cannot  be  successfially  matched  to  an  analyzed  eireulation,  the  foreeast 
eireulation  is  stored  as  a  potential  false  alarm  (bottom  left  box  of  Figure  2.1a). 

C.  QUALITY  CONTROL  AND  POST-PROCESSING 

As  with  any  automated  process,  quality  eontrol  measures  are  neeessary  when 

analyzing  model  fields  with  the  TCVTP.  The  program  identifies  only  eireulations  that 
meet  the  speeified  threshold  eriteria  (e.g.,  outer  elosed  relative  vortieity  eontour). 
Additionally,  individual  eireulations  are  matehed  to  form  traeks  only  when  the  translation 
speed  and  traek  orientation  threshold  eriteria  are  met.  Data  gaps  due  to  non-availability 
of  NOGAPS  data  fields  also  may  result  in  traeking  errors.  Sinee  the  TCVTP  program  is 
presently  not  eoded  to  deteet  sueh  data  gaps,  a  review  of  the  objeetively  analyzed  traeks 
is  required  to  eorreet  for  sueh  errors.  Additionally,  human  interaetion  is  required  to 
assoeiate  TCVTP  vortiees  with  speeifie  named  storms  if  eomparisons  of  model  foreeasts 
and  observed  storm  traeks  are  required. 
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m.  ANALYSIS  AND  RESULTS 


A.  TRACKED  VORTICITY  CIRCULATIONS 

The  TCVTP  was  used  to  track  all  circulations  meeting  the  threshold  8 50- mb 

relative  vorticity  criterion  in  the  tropical  western  North  Pacific  from  1  May  -  3 1  October 
2002.  The  resulting  254  tracked  circulations  were  further  sampled  to  remove  all 
circulations  that  did  not  survive  at  least  one  diurnal  cycle,  which  requires  existence  of  the 
circulation  center  in  three  consecutive  model  analyses  separated  by  12  h.  Hennon  and 
Hobgood  (2003)  used  a  similar  duration  criterion,  which  eliminates  short-term 
convectively- driven  circulations  and  one-time  spurious  circulations  spun  up  by 
anomalous  wind  observations.  By  using  this  minimum  duration  criterion,  the  number  of 
potential  vortices  was  decreased  to  127  (columns  one  and  two  in  Table  3.1).  The  average 
duration  of  these  127  vortices  was  4.2  days.  Vortices  that  did  not  meet  the  minimum 
duration  criterion,  but  did  satisfy  the  over- forecast  criteria  (ND2)  (discussed  later  in  this 
section),  are  included  in  column  3  of  Table  3.1.  Additionally,  each  circulation  identified 
by  the  TCVTP  was  assigned  to  one  of  three  formation  regions  based  on  the  initial  latitude 
less  than  30°N  and  longitude  (see  subregion  definitions  in  Table  3.1). 

The  vortices  that  met  the  minimum  duration  criterion  were  further  separated  into 
developing,  non- developing,  and  false  alarm  categories.  Storms  were  considered 
“developers”  if  at  least  Tropical  Depression  (TD)  strength  (Table  3.2)  was  reached,  and 


Table  3.1.  Circulations  tracked  in  NOGAPS  analyses  from  1  May  -  31  Oct  2002 
numbered  by  the  JTWC  (column  1),  non- developing  circulations  meeting  the  minimum 
duration  criterion  of  at  least  24  h  (column  2),  those  non- developing  circulations  not 
meeting  the  minimum  duration  criterion  but  meeting  the  over- forecast  criteria  (column 
3),  and  those  that  did  not  meet  the  minimum  duration  criterion,  but  had  850- mb  relative 
vorticity  forecasts  exceeding  5.0  x  10'^  s'^  (column  4),  or  judged  to  be  model  false  alarms 
(column  5)  in  each  of  the  tropical  subregions. 


Formation  Subregion 

Numbered 

TCs 

ND  >  24  h 
(NDl) 

ND  <  24  h 
(NDOE) 

ND  >  24  h 
(NDG2) 

Model 

Ealse 

Alarms 

South  China  Sea  (105E-124E) 

5 

34 

16 

19 

4 

Philippine  Sea  (125E-159E) 

10 

49 

28 

21 

10 

East  Monsoon  Trough 
(160E-180E) 

8 

21 

16 

10 

1 

Total 

23 

104 

60 

50 

15 

15 


the  circulation  was  thus  given  a  number  by  the  Joint  Typhoon  Warning  Center  (JTWC). 
During  the  period  of  study,  only  23  circulations  reached  at  least  TD  strength,  with  six 
each  reaching  STY  and  TY,  eight  reaching  TS  status,  and  three  reaching  only  TD  (Tables 
3.3  -  3.5).  Twenty-two  of  these  developing  cyclones  were  present  in  the  NOGAPS 
analyses  prior  to  the  first  JTWC  warning  time  (intensity  >  25  kt).  One  “surprise”  cyclone 
was  not  present  in  the  NOGAPS  analysis  prior  to  the  first  warning  time.  That  is,  the 
NOGAPS  analysis  did  not  include  the  storm  (TD  27)  until  a  synthetic  vortex  was  inserted 
upon  receipt  of  the  JTWC  warning  message.  Vortices  were  considered  “non- developers” 
if  they  did  not  reach  at  least  TD  strength  as  defined  by  JTWC. 

All  circulations  were  sorted  into  two  categories:  circulations  that  met  the 
minimum  duration  criterion  (NDl)  and  circulations  that  did  not  (ND2).  Circulations  that 
did  not  meet  the  minimum  duration  criterion  are  the  subset  of  NOGAPS-forecast  vortices 
that  could  be  matched  with  a  vortex  in  the  NOGAPS  analysis  for  less  than  the  required  24 
h,  but  were  forecast  to  live  much  longer.  If  the  NOGAPS  model  forecast  the  vortex  to 
exist  longer  than  24  hours,  but  the  analyzed  circulation  did  not  persist  for  24  hours,  then 
the  forecast  vortex  is  categorized  as  an  over- forecast  (NDOF). 

Model  false  alarms  (column  5  of  Table  3.1)  were  identified  as  forecast  vortices 
that  could  not  be  matched  to  a  circulation  in  the  verifying  analysis.  In  the  TCVTP 
system,  these  are  tracked  forecast  circulations  in  the  forecast  database  (Figure  2.1a)  that 
never  get  moved  to  the  analyzed  circulation  database.  The  number  of  false  alarms  in 
each  subregion  is  listed  in  column  five  of  Table  3.1.  Cheung  and  Elsberry  (2002)  also 
found  that  NOGAPS  was  more  skilllul  (fewer  false  alarms)  in  the  South  China  Sea  (SCS) 
and  Eastern  Monsoon  Trough  (EMT)  than  in  the  Philippine  Sea  (PS). 

Table  3.2.  Tropical  cyclone  designations  used  by  the  JTWC,  where  10-minute  average 
wind  speeds  are  used. _ 


Description 

Maximum  Sustained  Wind  Speed 

Tropical  Depression  (TD) 

12- 17  m  s' ^  (25-34  kt) 

Tropical  Storm  (TS) 

18-32  m  s' ^  (35-63  kt) 

Typhoon  (TY) 

33-66  m  s' ^  (64-129kt) 

Super  Typhoon  (STY) 

>  67  m  s"^  (130  kt) 
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1.  South  China  Sea  Formation 

The  presence  of  the  monsoon  trough  in  this  subregion  (McBride  1981)  makes  it  a 
favored  region  for  potential  tropical  cyclone  seedling  development  (Figure  3.1),  with 
30.7%  (39  of  127)  of  the  analyzed  vortices  during  the  study  period  present  in  the  SCS 
subregion.  The  frequent  occurrence  of  a  broad  area  of  background  cyclonic  vorticity 
between  cross- equatorial  westerlies  and  easterlies  to  the  north  contributed  to  the 
relatively  large  number  of  non- developing  circulations  in  the  NOGAPS  analyses. 
However,  the  number  of  vortices  that  developed  to  TD  strength  or  greater  was  only  five 
of  23  (21.7%)  (Table  3.3).  Consequently,  only  five  of  the  39  (12.8%)  potential  tropical 
cyclone  seedlings  in  the  NOGAPS  analyses  developed  to  warning  status.  Conversely, 
many  more  (87.2%)  of  the  potential  seedlings  in  the  SCS  in  the  NOGAPS  analyses  never 
achieve  TD  status.  These  may  be  weak  monsoon  depressions  (perhaps  with  considerable 
precipitation  and  winds  in  rain  bands  at  large  radii),  or  may  be  spurious  vortices  that 
persisted  in  the  analyses  for  24  h.  The  average  duration  of  both  developing  and  non¬ 
developing  vortices  in  this  basin  was  3.7  days,  which  is  less  than  in  the  other  basins  due 
to  the  small  size  and  limited  fetch  of  this  subregion. 


Figure  3.1.  Tracks  of  39  vortices  in  the  NOGAPS  analyses  that  formed  h  the  SCS 
(105°E-  125°E)  subregion  from  1  May  -  31  October  2002. 


17 


Table  3.3.  Maximum  intensity  (see  definitions  in  Table  3.2)  during  the  life  cycle  and 


formation  month  of  numbered  tropical  cyclones  in  the  SCS  subregion. 


May 

Jun 

Jul 

Aug 

Sep 

Oct 

Total 

TD 

1 

0 

0 

0 

0 

0 

1 

TS 

0 

0 

1 

0 

2 

0 

3 

TY 

0 

1 

0 

0 

0 

0 

1 

STY 

0 

0 

0 

0 

0 

0 

0 

Total 

1 

1 

1 

0 

2 

0 

5 

2.  Philippine  Sea  Formation 

The  Philippine  Sea  (PS)  subregion,  which  is  also  characterized  by  a  broad  area  of 
cyclonic  relative  vorticity  (McBride  1981)  associated  with  the  monsoon  trough,  was  the 
most  active  of  the  three  subregions  studied,  with  59  of  127  (46.5%)  circulations  (Figure 
3.2).  During  the  study  period,  10  of  the  23  vortices  (43.5%)  that  developed  to  TD 
strength  or  greater  formed  in  the  PS  region  (Table  3.4).  However,  these  10  named 
vortices  represent  only  a  16.9%  formation  rate  in  terms  of  the  circulations  in  the 
NOGAPS  analyses.  As  in  the  SCS,  many  of  the  non- developing  vortices  in  the  NOGAPS 
analyses  are  probably  monsoon  depressions  that  never  achieve  TD  status  with  central 
convection  and  high  winds.  The  average  duration  of  vortices  in  this  region  was  4.1  days. 


Figure  3.2.  Tracks  of  59  vortices  in  the  NOGAPS  analyses  that  formed  in  the  PS  (125°E 
-  160°E)  subregion  from  1  May  -  31  October  2002. 
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Table  3.4.  Maximum  intensity  by  formation  month  as  in  Table  3.3,  except  for  numbered 


tropical  cyclones  in  the  PS  subregion. 


May 

Jun 

Jul 

Aug 

Sep 

Oct 

Total 

TD 

0 

0 

0 

0 

0 

0 

0 

TS 

0 

0 

1 

2 

0 

1 

4 

TY 

0 

1 

1 

1 

0 

1 

4 

STY 

1 

1 

0 

0 

0 

0 

2 

Total 

1 

2 

2 

3 

0 

2 

10 

3.  East  Monsoon  Trough  Formation 

The  east  monsoon  trough  (EMT)  subregion,  which  extends  from  160°E  to  180°E, 
had  29  of  the  127  (22.8%)  circulations  tracked  by  the  TCVTP  in  the  NOGAPS  analyses 
during  the  study  period  (Eigure  3.3).  During  the  study  period,  only  eight  of  the  23 
(34.8%)  vortices  that  developed  to  TD  strength  or  greater  formed  in  the  EMT  subregion 
(Table  3.5).  These  eight  JTWC-numbered  circulations  represent  27.6%  of  the  29  tracked 
vortices  ^erifying  in  the  NOGAPS  analyses.  The  average  duration  of  vortices  within 
general  area  of  the  east  end  of  the  monsoon  trough  was  5.0  days,  which  is  the  longest  of 
the  three  subregions  studied.  Although  this  subregion  is  also  relatively  small,  vortex 
duration  is  greater  than  in  the  other  two  regions  due  to  the  large  distance  from  land. 
Perhaps  for  the 


Eigure  3.3.  Tracks  of  29  vortices  in  the  NOGAPS  analyses  that  formed  in  the  EMT 
(160°E-  180°E)  subregion  from  1  May  -  31  October  2002. 


19 


same  reason,  the  cyclones  that  developed  within  the  EMT  subregion  also  achieved,  on 
average,  a  greater  maximum  intensity  than  cyclones  that  developed  in  the  other  two 
regions.  One  factor  favoring  formation  and  intensification  within  this  region  is  the 
extensive  area  of  cyclonic  vorticity  associated  with  the  oceanic  monsoon  trough  when 
tropical  cyclones  form  at  the  eastern  end. 


Table  3.5.  Maximum  intensity  and  formation  month  as  in  Table  3.3,  except  for  numbered 
tropical  cyclones  in  the  EMT  subregion. 


May 

Jun 

Jul 

Aug 

Sep 

Oct 

Total 

TD 

0 

0 

1 

0 

0 

1 

2 

TS 

0 

0 

1 

0 

0 

0 

I 

TY 

0 

0 

0 

1 

0 

0 

I 

STY 

0 

0 

2 

1 

1 

0 

4 

Total 

0 

0 

4 

2 

I 

I 

8 

B.  TROPICAL  CYCLONE  FORMATION  ALERT  MESSAGES 

In  this  section,  an  overall  summary  of  tropical  cyclone  formation  is  given  with 

respect  to  official  forecasts  and  tropical  cyclone  formation  alert  (TCEA)  messages  issued 
by  the  JTWC.  This  summary  is  provided  as  a  baseline  by  which  the  model  forecast 
performance  of  tropical  cyclone  formation,  which  is  discussed  in  the  next  section,  may  be 
compared. 

The  JTWC  issues  a  TCEA  based  on  wind  and  pressure  analyses  at  the 
surface/gradient  level,  500  mb,  and  200  mb,  and  a  subjective  Dvorak  analysis  of  satellite 
imagery  (see  Appendix  A).  Such  formation  alerts  are  valid  for  24  hours,  and  they  may  be 
either  extended  or  cancelled  at  the  end  of  the  valid  period  depending  on  the  current  status 
of  the  circulation  center.  If  the  circulation  develops  to  tropical  depression  strength  (Table 
3.2),  the  first  warning  issued  for  the  system  replaces  the  TCEA.  All  TCEA  messages 
were  gathered  from  the  Automated  Tropical  Cyclone  Eorecast  (ATCE)  database 
maintained  at  the  Naval  Research  Eaboratory  -  Monterey. 

The  TCEA  applies  within  a  formation  box  specified  as  an  area  on  either  side  of  a 
line  defined  by  two  grid  points.  Alternately,  a  formation  circle  may  be  specified  by  a 
single  grid  point  and  radial  distance.  The  size  of  the  formation  box  or  circle  is 
subjectively  determined  by  the  forecaster’s  confidence  in  the  objective  aids.  The 
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temporal  window  during  which  tropical  cyclone  formation  is  expected  to  occur  is 
provided  in  each  TCFA.  Additionally,  an  estimate  of  the  wind  speed,  translation  speed, 
and  heading  of  the  developing  circulation  is  provided,  along  with  the  location  of  the 
center  of  the  circulation  determined  from  satellite  imagery.  Latitude  and  longitude 
positions  contain  a  “checksum”  at  the  end  of  each  string.  The  checksum  is  calculated  by 
summing  the  numbers  in  the  latitude  or  longitude  value,  and  then  the  last  digit  of  the 
value  is  listed  at  the  end  of  the  string,  e.g.,  latitude  147.7E9  sums  to  19,  and  9  is  encoded 
as  the  checksum.  A  sample  of  the  text  portion  of  the  TCFA  (for  STY05  -  Hagibis) 
follows: 

FORMATION  OF  A  SIGNIFICANT  TROPICAL  CYCLONE  IS 
POSSIBLE  WITHIN  160  NM  EITHER  SIDE  OE  A  LINE  EROM 
3.5N7  147. 7E9  TO  9 . 4N3  145. 6E6  WITHIN  THE  NEXT  06  TO 
24  HOURS.  AVAILABLE  DATA  DOES  NOT  JUSTIEY  THE 
ISSUANCE  OE  NUMBERED  TROPICAL  CYCLONE  WARNINGS  AT  THIS 
TIME.  WINDS  IN  THE  AREA  ARE  ESTIMATED  TO  BE  15  TO  22 
KNOTS.  METSAT  IMAGERY  AT  I3I430Z2  INDICATES  THAT  A 
CIRCULATION  CENTER  IS  LOCATED  NEAR  4 . 7NI  148. 9E2.  THE 
SYSTEM  IS  MOVING  NORTHWESTWARD  AT  10  KNOTS. 

During  1  May  -  31  Oct  2002,  the  JTWC  issued  32  TCFAs  (excluding  the  TCFAs 
for  TD17  and  TD27).  For  six  storms  (STY05,  TY07,  TY09,  TD15,  TS24,  and  TY26), 
two  TCFAs  were  issued  with  nominal  changes  between  the  first  and  the  second.  In 
general,  changes  were  made  to  extend  the  formation  forecast  window  if  the  formation  had 
not  occurred  within  the  specified  time,  but  was  still  expected  to  occur.  The  formation 
box  or  circle  was  also  slightly  modified  to  account  for  motion  of  the  developing 
circulation.  Six  of  these  32  TCFAs  were  false  alarms  (two  in  August,  and  four  during 
one  week  September).  For  these  six  TCFAs,  the  circulation  center  dissipated  or 
weakened  significantly,  and  therefore  a  warning  was  never  issued.  In  two  cases  (STYIO 
and  TSI6),  no  prior  TCFA  had  been  issued  before  the  first  warnings  were  issued. 

1.  Lead  Time 

Fead  time  (?  t)  was  defined  as  the  time  between  the  issuance  of  a  TCFA  and  the 
first  warning  time.  The  average  ?t  for  successful  TCFAs  was  12.1  hours.  Maximum  and 
minimum  ?t  values  ranged  from  zero  hours  for  STY05  to  42  hours  for  TY09  and  TD15 
(Figure  3.4).  TCFAs  were  reissued  for  all  TCFAs  that  did  not  verify  within  24  h  to 
extend  the  formation  spatial  and  temporal  windows  (R.  Feejoice,  JTWC,  personal 
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Figure  3.4.  Lead  time  (?t)  of  the  TCFAs  for  developing  storms.  Bars  indieate  the 
number  of  hours  between  the  TFCA  issue  time  and  the  first  JTWC  warning  time.  The  ?t 
for  the  first  TCFA  for  STY05  and  STY25  and  seeond  TCFA  for  TY  07  was  less  than  1 
hour.  TCFA  messages  were  not  issued  for  STYIO  and  TS16.  Repeated  numbers  indieate 
multiple  TCFAs  issued  for  a  single  storm. 

communication).  No  ?t  was  calculated  for  STYIO  or  TS16  since  no  TCFA  was  issued 
prior  to  the  first  warning  for  either  the  storm. 

2.  Formation  Window 

A  formation  window,  which  indicates  when  a  circulation  is  mqrected  to  reach 
tropical  depression  strength,  is  specified  in  each  TCFA.  This  window  is  generally  6-24 
hours  but  can  range  from  8-24,  or  12-24  hours.  The  onset  of  the  formation  window  is 
more  a  function  of  the  individual  forecaster  than  of  model  certainty  (S.  Vilpors,  JTWC, 
personal  communication).  Verification  of  the  formation  window  (Table  3.6),  which  is 
scored  as  a  success  if  the  issuance  of  the  first  JTWC  warning  occurred  within  that 
window,  which  only  occurred  in  11  (41%)  of  the  27  available  TCFAs.  Of  the  16  (59%) 
TCFAs  that  failed  to  result  in  a  formation  within  the  specified  time  window,  11  (61%) 
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first  warnings  were  issued  prior  to  the  onset  of  the  formation  window  (i.e.,  formed  too 
early),  and  five  (19%)  first  warnings  occurred  more  than  24  hours  after  the  TCFA  was 
issued  (i.e.,  formed  too  late).  All  five  of  these  late  forecasts  were  reissued  to  extend  the 
formation  window.  Of  the  11  TCFAs  that  verified  prior  to  the  onset  of  the  formation 
window,  nine  (82%)  of  the  first  warnings  were  less  than  4  hours  from  the  time  of  TCFA 
issue. 


Table  3.6.  Summary  of  formation  times  (first  JTWC-issued  warning)  relative  to  TCFA- 
specified  formation  window. 


Formation  Time 

Number  of  Occurrences 

Outside  of  specified  window 

16 

Early 

11 

Late 

5 

During  specified  window 

8 

Total 

24 

C.  VERIFICATION  OF  FORECAST  VARIABLES 

Thirteen  variables  listed  in  Table  2.1  were  recorded  based  on  the  vortex 

representation  defined  by  the  TCVTP.  Of  those  13  variables,  five  were  selected  as 
having  some  measure  of  utility  for  forecasting  the  formation  of  the  developing  vortices. 
Those  six  variables  were:  850- mb  relative  vorticity  (VOR),  200-850  mb  or  deep  wind 
shear  (DSH),  sea- level  pressure  (SLP),  925-mb  wind  speed  (925),  and  500  -  700  mb 
average  vapor  pressure  (VPR).  Whereas  the  database  contained  three  wind  speed 
variables  (at  925,  700,  and  500  mb)  that  were  potentially  useful  in  forecasting  tropical 
cyclone  formation,  the  700-  and  500-mb  wind  speed  variables  did  not  show  a  distinct 
signal  for  identification  of  developing  vortices,  and  so  were  not  added  to  the  list  of 
selected  variables.  Several  other  variables  were  not  selected  since  they  represent  derived 
quantities  that  are  based  on  model  formulations  rather  than  analyzed  variables.  These 
model- diagnostic  variables  did  not  appear  to  contain  as  much  information  regarding  the 
formation  condition  as  those  variables  based  on  analyzed  quantities.  Other  measured 
variables,  i.e.,  200-850  mb  thickness,  and  500-850  mb  wind  shear,  showed  little  or  no 
signal  for  distinguishing  developing  from  non- developing  vortices,  and  so  were  not 
included  in  the  list  of  selected  variables. 
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The  five  selected  variables  were  then  analyzed  to  determine  their  utility  at 
forecasting  the  potential  development  or  non- development  of  the  circulation  dataset.  The 
impact  of  the  selected  variables  on  the  developing  vortices,  i.e.,  those  numbered  by 
JTWC  during  1  May  -  31  Oct  2002,  were  analyzed  first.  Each  selected  variable  was 
examined  at  12-h  intervals  extending  from  120  h  prior  to  and  120  h  after  the  first  JTWC 
warning  and  first  Best- Track  times.  For  each  time  step,  the  range,  mean,  and  standard 
deviation  of  each  variable  were  calculated.  Threshold  values  for  each  variable  at  the  first 
JTWC  warning  time  and  first  Best- Track  time  were  determined,  and  will  be  used  as  a 
benchmark  for  the  comparison  of  non- developing  vortices. 

These  statistics  will  be  displayed  based  on  the  NOGAPS  forecast  verification  time 
relative  to  the  first  JTWC  warning  (e.g..  Figure  3.5)  or  the  first  Best- Track  time  (e.g.. 
Figure  3.6).  Since  these  are  the  first  of  many  such  displays,  some  extended  comments  on 
Figure  3.5  will  be  made  to  introduce  the  figures.  A  representative  sampling  of  forecast 
times  (i.e.,  24  h,  48  h,  72  h,  and  96  h)  are  presented  in  the  four  panels.  For  each  forecast 
interval,  the  average  analyzed  (at  -i-OO)  value  is  plotted  against  the  average  forecast  value 
for  the  indicated  forecast  verification  time,  with  the  range,  mean,  and  standard  deviation 
at  each  time  interval  relative  to  the  first  JTWC  warning  time  in  Figure  3.5.  The  average 
forecast  value  (red  line)  is  only  displayed  if  five  or  more  cases  were  available  at  each 
time  interval.  This  minimum  sample  size  criterion  will  be  relaxed  for  a  handful  of 
variables  in  which  fewer  than  five  cases  were  available  for  analysis.  Statistical 
significance  was  calculated  for  the  differences  in  the  means  of  the  forecast  and  analyzed 
variable  using  a  two-tailed  t-test,  and  significant  differences  are  indicated  in  the 
following  figures  by  a  dark- shaded  box  plot. 

Consider  first  the  24-h  forecast  verifications  in  the  upper  left  panel  of  Figure  3.5. 
A  difference  between  the  mean  forecast  (heavy  red  line)  and  analyzed  (heavy  black  line) 
850-mb  vorticity  of  the  circulations  is  first  noted  at  36  h  prior  to  the  first  JTWC  warning 
time.  These  mean  values  then  are  the  averages  of  all  24-h  forecasts  verifying  at  36  h 
prior  to  the  first  warning  time,  which  means  that  these  24-h  forecasts  were  initiated  at  60 
h  prior  to  the  first  warning  time.  Fooking  at  the  sample  size  at  -36  h,  only  18  such  24-h 
forecasts  are  available  from  the  23  developing  storms,  which  means  that  in  five  cases, 
even  a  24-h  forecast  initiated  60  h  prior  to  the  first  warning  time  did  not  exist  for 
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Figure  3.5.  Average  analyzed  and  foreeast  850- mb  relative  vortieity  for  developing 
vortiees  relative  to  the  first  JTWC  warning  time  (Fq).  Heavy  black  line  represents 
analyzed  average  850- mb  relative  vortieity  (at  +00),  and  red  line  represents  average 
forecast  850- mb  relative  vortieity  for  the  indicated  forecast  time.  Average  850- mb 
relative  vortieity  at  Fo  (light  solid  black  line)  was  5.0  x  10'^  s'\  which  is  then  a  reference 
for  the  vortieity  magnitudes  before  (to  the  left)  and  after  (to  the  right)  the  first  warning 
time.  The  sample  sizes  available  for  comparison  of  the  analyzed  and  forecast  values  at 
each  forecast  verification  time  are  listed  near  the  top  of  each  panel  for  24- h  (upper- left), 
48-h  (upper- right),  72-h  (lower- left),  and  96-h  (lower- right)  forecasts  of  these  developing 
storms. 


verification  36  h  prior  to  that  first  warning  time.  Even  larger  fall-offs  in  sample  size  are 
noted  for  the  24- h  forecasts  that  are  available  for  comparison  with  the  analyzed  values  at 
the  earlier  forecast  verification  times.  For  the  sample  of  24-h  forecast  in  the  upper-left 
panel,  only  four  cases  were  available  for  verification  120  h  prior  to  the  first  warning  time, 
because  these  forecasts  would  have  been  initiated  144  h  prior  to  the  first  warning  time 
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when  either  that  pre-tropical  cyclone  circulation  could  not  be  tracked  in  the  forecast  or  in 
the  analysis  so  early  in  the  life  cycle. 

Although  the  24- h  forecast  vorticity  value  is  larger  than  the  analyzed  vorticity 
value  at  36  h  prior  to  the  first  warning  time,  the  overlapping  box  plots  at  this  time  suggest 
that  the  difference  is  not  significant.  Notice  also  that  the  range  of  forecast  values  is 
relatively  large  compared  to  the  difference  between  the  means  of  the  24- h  forecast 
vortices  and  the  corresponding  analyzed  vortices  at  36  h  prior  to  the  first  warning  time. 
An  outlier  forecast  vorticity  of  about  8  x  10'^  s'^  is  indicated  by  the  open  circle  at  the  -36 
h  verification  time  in  the  upper- left  panel  of  Figure  3.5,  which  means  that  one  of  the  24-h 
forecasts  initiated  60  h  prior  to  the  first  warning  time  grossly  over-forecast  the  850-mb 
vorticity  compared  to  the  average  threshold  value  of  5.0  x  10"^  s'^  for  all  23  developing 
vortices  at  the  first  warning  time. 

Consider  next  the  24-h  forecast  verifications  at  00  h  in  the  upper  left  panel  of 
Figure  3.5.  Here,  the  24-h  analyzed  and  forecast  vorticity  values  are  almost  identical, 
which  means  that  the  NOGAPS  forecasts  initiated  only  24  h  prior  to  the  first  warning 
time  are  quite  accurate.  A  more  accurate  24-h  forecast  verifying  at  the  first  warning  time 
compared  to  verifying  36  h  prior  to  the  first  warning  time  (as  discussed  above)  might  be 
expected  if  the  initial  conditions  at  24  h  prior  to  warning  time  are  more  accurate  than  the 
initial  conditions  at  60  h  prior  to  the  first  warning  time.  However,  the  size  of  the  box 
plots  and  range  of  24-h  forecast  vortices  verifying  at  00  h  are  relatively  large,  which 
means  considerable  variability  exists  among  the  22  cases  at  this  verification  time. 

The  only  statistically  significant  difference  between  the  mean  24-h  analyzed  and 
forecast  vorticity  in  the  upper- left  panel  in  Figure  3.5  occurs  for  forecasts  verifying  24  h 
after  the  first  warning  time.  At  this  time,  the  mean  24-h  forecast  vorticity  is  much 
smaller  than  the  verifying  vorticity,  and  for  the  21  forecasts  verifying  at  this  time  the 
difference  in  the  means  is  statistically  significant  according  to  the  two-tailed  t  test. 

The  above  discussion  indicates  the  amount  of  information  that  can  be  extracted 
from  the  24-h  forecasts  in  the  upper- left  panel  of  Figure  3.5.  Similar  interpretations  can 
be  made  for  the  48- h  (upper- right  panel),  72-h  (lower- left  panel),  or  96- h  (lower- right 
panel)  forecast  verifying  at  the  various  times  relative  to  the  first  JTWC  warning  time. 
These  other  panels  illustrate  the  capability  of  the  NOGAPS  model  to  forecast  these 
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developing  vortices  at  longer  and  longer  forecast  intervals,  and  thus  how  useful  the 
NOGAPS  model  will  be  in  forecasting  the  formation  defined  as  occurring  at  the  first 
JTWC  warning  time.  Similar  interpretations  can  be  made  in  Figure  3.6  for  an  earlier 
definition  of  forecast  time  as  being  the  first  Best- Track  time.  In  these  two  examples,  the 
interpretations  are  with  respect  to  the  850- mb  relative  vorticity  for  just  the  developing 
vortices,  but  other  variables  that  might  characterize  the  differences  between  developing 
and  non- developing  vortices  will  also  be  presented  in  similar  plots  in  the  following 
sections. 

Appendix  B  contains  a  more  detailed  explanation  of  the  data  contained  in  each  of 
the  different  types  of  plots  presented  in  the  following  sections.  Readers  who  are  not 
familiar  with  such  forecast  verification  procedures  are  advised  to  read  Appendix  B  before 
continuing. 

1.  Developing  Vortices 

The  time  of  the  first  warning  issued  for  each  numbered  storm  (designated  Fo)  is 
one  possible  definition  of  formation  time  and  was  determined  from  the  2002  Annual 
Tropical  Cyclone  Report  (ATCR).  An  alternate  formation  time  is  the  time  of  the  initial 
Best- Track  position  for  each  numbered  storm  (designated  F*o)  and  was  also  determined 
from  the  2002  ATCR.  The  initial  Best-Track  position,  intensity,  and  time  are  determined 
during  post-storm  analysis  by  the  JTWC,  and  correspond  to  the  first  model  analysis, 
subjective  analysis  or  satellite  image  in  which  the  developing  circulation  is  apparent.  It 
should  be  noted  that  both  Fo  and  F*o  occur  prior  to  the  first  insertion  of  any  synthetic 
tropical  cyclone  observations  into  the  NOGAPS  model,  which  occurs  when  the  maximum 
analyzed  wind  speed  first  reaches  25  kt.  Threshold  values  for  these  two  potential 
definitions  of  formation  time  were  determined  to  be  the  mean  values  of  each  variable  in 
the  NOGAPS  analyses  at  the  times  of  Fb  and  F*o-  For  a  developing  storm,  F*o  must 
either  precede  or  equal  Fq. 

a.  850-mb  Relative  Vorticity 

(1)  Entire  Basin  Assessment.  Low-level  (850-mb)  relative  vorticity 
(?)  provided  the  strongest  signal  for  identification  of  developing  versus  non- developing 
vortices.  A  positive  value  indicates  a  low-level  cyclonic  circulation,  and  as  expected, 
higher  850-mb  relative  vorticity  is  analyzed  at  Fq  and  F*o  for  those  vortices  that 
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developed  into  numbered  circulations  than  for  those  that  did  not  develop,  which  will  be 
shown  in  a  later  section. 

Prior  to  Fq,  the  NOGAPS  forecasts  of  850- mb  relative  vorticity  at  the  various 
forecast  intervals  shown  by  the  panels  in  Figure  3.5  were  consistently  larger  than 
analyzed,  whereas  after  Fo  the  vorticity  forecasts  were  less  than  analyzed.  This  transition 
between  over-  and  under- forecasts,  which  occurred  within  the  first  12  hours  after  F),  is 
likely  due  to  the  insertion  of  the  synthetic  tropical  cyclone  observations  into  the 
NOGAPS  analysis.  When  the  first  warning  for  a  developing  storm  is  issued  by  the 
JTWC,  an  automated  message  is  received  at  the  Fleet  Numerical  Meteorology  and 
Oceanography  Center  (FNMOC).  The  Navy  data  assimilation  system  (NAVDAS) 
searches  hourly  for  these  automated  messages.  If  one  is  received,  and  it  meets  certain 
quality  control  criteria,  then  a  pre-set  synthetic  tropical  cyclone  wind  and  sea- level 
pressure  observations  at  13  points  centered  on  the  reported  storm  position  is  inserted  into 
the  file  for  use  in  the  next  NOGAPS  update  cycle.  The  location,  size,  and  magnitude  of 
the  synthetic  wind  observations  are  based  on  the  data  contained  in  the  original  warning 
from  the  JTWC.  Although  the  NOGAPS  update  cycle  is  run  hourly,  and  changes  are 
frequently  made  to  the  model  analysis,  the  first  time  the  analysis  will  include  the  effect  of 
the  synthetic  observations  is  in  the  next  analysis  after  the  first  warning  is  issued.  For 
example,  if  JTWC  issues  the  first  warning  on  a  developing  storm  at  0200  UTC,  the 
synthetic  observation  is  ingested  into  the  0300  UTC  update  cycle,  and  will  first  be 
apparent  in  the  0600  UTC  model  analysis  (B.  Strahl,  FNMOC,  personal  communication). 

The  mean  analyzed  850- mb  relative  vorticity  at  Fq  (?o)  was  5.0  x  10"^  s'\  which  is 
indicated  by  the  thin,  solid  horizontal  line  in  Figure  3.5.  Although  the  slope  of  the 
analyzed  ?o  curve  changed  after  Fq  due  to  the  insertion  of  the  synthetic  tropical  cyclone 
observations,  ihe  slope  of  the  forecast  ?  (?i2,  ?24,-.-)  curves  generally  do  not  have  a 
marked  change  in  slope  relative  to  Fq.  However,  a  tendency  is  noted  for  the  slope  of  the 
curves  to  flatten  for  forecasts  initiated  prior  to  the  insertion  of  the  synthetic  tropical 
cyclone  observations. 

The  impact  of  the  synthetic  tropical  cyclone  observations  starts  from  the  first 
analysis  after  Fq  when  JTWC  has  numbered  the  storm  and  issued  its  first  warning.  The 
insertion  of  the  synthetic  tropical  cyclone  observations  is  obvious  in  the  forecast  relative 
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vortkity  curves  in  Figure  3.5.  The  slope  of  the  forecast  ?  curve  in  the  upper- left  panel  of 
Figure  3.5,  which  represents  the  average  forecast  vorticity  value  of  +24  h  forecasts,  is 
almost  flat  from  F)  to  Fo-i-24  h,  and  then  begins  to  increase  after  24  h.  The  dramatk 
increase  in  slope  about  12  h  after  the  specified  forecast  time  interval  represents  the 
response  in  the  NOGAPS  forecasts  to  the  synthetic  tropical  cyclone  observations. 
Forecasts  made  after  the  synthetic  tropical  cyclone  observations  were  inserted  are  more 
accurate  when  compared  to  the  analyzed  relative  vorticity  values.  Notice  that  the  amount 
by  which  the  forecast  ?  curves  fall  below  the  analyzed  curve  after  Fq  generally  increases 
with  larger  forecast  range.  Statistically  significant  differences  between  the  analyzed 
means  and  forecast  means  are  indicated  in  Figure  3.5  by  darkened  box  plots. 

This  model  response  hypothesis  is  apparent  for  shorter-range  forecasts,  and  is 
easily  seen  in  the  upper-  and  lower- left  panels  of  Figure  3.5.  At  longer  time  ranges,  the 
model  response  time  is  delayed  by  about  12  hours.  The  slight  decrease  in  the  -i-48-h 
forecast  ?  in  the  upper-right  panel  of  Figure  3.5  prior  to  the  insertion  of  the  synthetic 
tropical  cyclone  observations  is  attributed  to  the  NOGAPS  model  already  beginning  to 
intensify  the  developing  tropical  circulation  prior  to  the  insertion  of  the  synthetic  tropical 
cyclone  observations.  The  relative  failure  of  the  response  hypothesis  at  extended  forecast 
times  could  be  due  to  decreased  dynamic  predictability  of  the  tropical  atmosphere  at 
longer  forecast  ranges,  or  the  sample  size  is  too  small. 

Because  the  NOGAPS  vorticity  forecasts  about  12  hours  after  Fq  are  consistently 
less  than  the  average  analyzed  vorticity,  this  indicates  a  tendency  for  the  NOGAPS  model 
to  under- forecast  developing  vortices  after  synthetic  tropical  cyclone  observations  are 
included  in  the  analysis.  This  decrease  in  the  magnitude  of  vorticity  may  be  due  to  the 
structure  of  the  synthetic  tropical  cyclone  observations.  Prior  to  insertion  of  the  synthetic 
observations,  the  vorticity  ellipse  may  be  too  large,  and  in  addition  may  be  too  strong. 
The  synthetic  tropical  cyclone  observations  force  the  circulation  in  the  model  to  decrease 
in  size  and  magnitude  by  concentrating  the  relative  vorticity  in  the  NOGAPS  forecast. 

The  average  time  difference  between  Fq  and  the  first  Best- Track  position  F*o  was 
28.4  hours,  and  ranged  from  6  hours  to  84  hours,  with  a  median  value  of  24  hours.  This 
time  difference  between  when  the  analyzed  curve  (heavy  black  solid)  (Figure  3.6)  crosses 
the  F*o  (light  dashed)  vorticity  threshold  and  when  it  crosses  the  Fq  (light  solid)  vorticity 
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Figure  3.6.  Average  analyzed  and  forecast  850- mb  relative  vorticity  for  developing 
vortices  as  in  Figure  3.5,  except  relative  to  the  first  Best- Track  time  (F*o)-  Heavy  black 
line  represents  analyzed  average  850-mb  relative  vorticity  (at  -i-OO),  and  red  line 
represents  average  forecast  850-mb  relative  vorticity  for  the  indicated  forecast  time. 
Average  850-mb  relative  vorticity  at  Fo  (light  solid  black  line)  was  5.0  x  10"^  s'*,  and  4.3 
X  10'^  s'*  (hght  dashed  black  line)  at  F*o- 


threshold  appears  to  be  about  18  hours  for  the  short-range  forecasts.  The  10  h  difference 
between  the  graphical  determination  and  calculation  of  timing  difference  is  attributed  to 
the  varying  F*o  times  for  each  developing  circulation.  Some  circulations  were  evident  in 
the  NOGAPS  model  analyses  for  42  h  or  more  before  the  JTWC  issued  the  first  warning. 
Removing  the  two  greatest  F*o  times  (54  and  84  h)  reduces  the  mean  to  24.6  h  and  the 
median  to  18  h,  which  is  roughly  what  was  estimated  from  Figure  3.6.  The  point  at 
which  the  short-range  forecast  ?  (heavy  red  solid)  curves  transition  from  over-  to  under¬ 
forecasts  is  approximately  12  to  24  hours  after  F*o-  Given  that  in  Figure  3.5  forecasts  of 
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relative  vorticity  transitioned  from  over-  to  under- forecast  at  Fq,  this  interval  is  close  to 
the  expected  value  given  the  roughly  24  hours  separating  Fq  and  F*o- 

As  in  Figure  3.5,  the  850- mb  relative  vorticity  (?*)  forecasts  in  Figure  3.6  made 
prior  to  insertion  of  the  synthetic  tropical  cyclone  observations  (at  Fq)  are  consistently  too 
large,  whereas  after  F)  the  forecast  vorticities  are  less  than  the  analyzed  vorticity.  The 
mean  analyzed  850- mb  relative  vorticity  atF*o  (?*o)  was  4.3  x  10'^  s‘\  which  is  indicated 
by  the  thin  dashed  horizontal  line  in  Figure  3.6.  As  expected,  ?*o  is  less  than  ?o,  which 
reflects  the  increase  in  vorticity  with  time  due  to  the  development  of  the  circulation  that 
occurs  during  the  roughly  24  hours  separating  Fq  and  F*o.  As  in  Figure  3.5,  the  insertion 
of  the  synthetic  tropical  cyclone  observations  is  evidenced  by  the  flattening  and  then 
sudden  increase  in  the  slope  of  the  forecast  ?*  curves  (Figure  3.6).  The  time  separation 
between  Fq  and  F*o  accounts  for  the  temporal  lag  in  the  increase  of  the  forecast  ?*  curve. 
Since  the  forecast  ?*  slope  is  not  as  dramatic  at  longer  forecast  ranges,  it  is  possible  that 
the  NOGAPS  model  is  less  adept  at  correctly  developing  tropical  circulations  96  h  prior 
to  the  eventual  first  Best- Track  time. 

To  examine  whether  the  model  forecast  trends  in  vorticity  might  be  due  to  a 
systematic  error  in  the  forecast  vortex  size  due  to  the  relatively  coarse  resolution  in 
NOGAPS,  forecast  size,  which  is  measured  by  the  number  of  grid  cells  within  the 
TCVTP  ellipse,  was  examined  (Figure  3.7).  Based  on  the  analyzed  (at  -i-OO)  values  of 
vortex  size  (heavy  black  line  in  Figure  3.7)  while  the  vortex  is  intensifying  prior  to  Fq,  it 
is  also  growing  in  size.  In  the  6  hours  prior  to  Fq,  the  size  decreases  before  beginning  to 
increase  again  about  36  hours  after  Fq.  This  size  decrease  prior  to  Fo  may  be  due  to  a 
concentration  of  vorticity  as  the  circulation  becomes  more  organized,  or  may  be  due  to  a 
refocusing  of  the  TCVTP  ellipse  on  the  most  intense  portion  within  a  broader  area  of 
relative  vorticity.  In  the  12  hours  immediately  following  F  (Figure  3.7),  the  slight 
decrease  in  analyzed  vortex  size  is  attributed  to  the  insertion  of  the  synthetic  vortex. 
Thirty- six  hours  after  Fq,  the  vortex  size  increases  as  the  vortex  intensifies  as  a  numbered 
system. 

Since  some  of  the  NOGAPS  forecasts  of  vortex  size  are  larger  and  some  smaller 
than  analyzed  (Figure  3.7),  no  significant  trends  toward  over-  or  under- forecast  of  size 
relative  to  Fo  are  evident,  although  the  effect  of  the  synthetic  tropical  cyclone 
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Figure  3.7.  Average  vortex  size  for  all  developing  vortices  relative  to  Fq.  Average  size  at 
Fo  (light  solid  black  line)  was  60.  Heavy  black  line  represents  average  analyzed  size  (at 
+00);  heavy  red  line  represents  average  size  at  forecast  time  indicated.  The  display  of 
boxes,  ranges  and  outlier  values  are  similar  to  Figure  3.5,  and  the  same  forecast  intervals 
are  shown  in  the  various  panels. 


observations  is  again  apparent.  However,  the  short-range  (less  than  60  h)  NOGAPS 
forecasts  of  vortex  size  tend  to  be  larger  than  analyzed  until  about  36  hours  after  Fq, which 
may  be  attributed  to  the  NOGAPS  model  tendency  to  inaccurately  represent  the  size  of 
the  circulation  because  of  the  coarse  model  resolution  (Goerss  and  Jeffries  1994).  The 
short-term  forecasts  of  vortex  size  dip  to  the  analyzed  vortex  size  threshold  (60  grid  cells) 
between  24  and  36  hours  after  Fo,  which  may  indicate  a  delay  in  the  model  response  to 
the  inserted  synthetic  tropical  cyclone  observations. 


32 


Timing  error  (?  ti)  was  calculated  as  the  amount  of  time  that  lapsed  between  Fb 
and  when  the  forecast  850-mb  vorticity  curve  crossed  the  Fo  vorticity  threshold  in  Figure 
3.5  (Table  3.7).  Since  average  values  from  Figure  3.5  are  used  to  generate  the  timing 
error  in  the  second  column  of  Table  3.7,  ?|  is  an  average  accurate  to  within  +  6  hours. 
As  discussed  above,  the  forecast  850-mb  relative  vorticity  exceeded  the  Fo  threshold  prior 
to  Fo,  which  indicates  that  the  model  over- forecast  the  relative  vorticity  by  forecasting  the 
circulation  to  reach  the  relative  vorticity  threshold  determined  for  developing  vortices 
(5.0  X  10'^  s'^)  prior  to  when  it  actually  reached  that  value.  Although  the  timing  error  for 
all  forecast  times  was  -16.7  hours,  the  average  timing  error  for  the  short-term  forecasts  (< 
60  hours)  was  just  over  -6  hours,  which  is  approximately  the  resolution  of  the  data.  As 
expected,  timing  error  increases  for  increasing  forecast  range,  although  due  to  the  small 
sample  size  the  decrease  at  -1-120  h  may  not  be  representative. 


Table  3.7.  Timing  error  (?|)  between  the  time  when  a  forecast  850-mb  relative  vorticity 
line  crosses  the  Ft)  vorticity  threshold  from  Figure  3.5  and  It).  If  the  forecast  curve 
crossed  the  Fo  threshold  between  x-axis  values,  the  mean  time  between  those  two  x-axis 
values  was  assigned.  Average  timing  error  is  -16.7  hours.  A  negative  number  indicates 
that  the  forecast  vorticity  curve  crossed  the  Fq  vorticity  threshold  prior  to  Fq  (-i-OO  on  x- 
axis  of  Figure  3.5).  No  value  is  calculated  for  the  -i-OO  forecast  time  because,  by 
definition,  the  threshold  value  is  determined  by  the  value  of  the  analyzed  curve  (-i-OO)  at 
Fo. 


Forecast  Time 

?ti 

-too 

-- 

-1-12 

-6 

-1-24 

-6 

-t36 

-6 

-1-48 

-12 

-t60 

-6 

-1-72 

-6 

-1-84 

-18 

-t96 

-54 

-1-120 

-36 

A  second  timing  error  (?fe)  was  calculated  as  the  difference  between  when  the 
forecast  vorticity  curve  crossed  the  Fb  vorticity  threshold  and  when  it  crossed  the  F*o 
vorticity  threshold  in  Figure  3.6  (Table  3.8).  The  average  ?b  was  29.4  hours,  which 
indicates  that  the  F*o  threshold  was  crossed  prior  to  the  F)  threshold  (Figure  3.6),  as 
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expected  based  on  the  definitions  of  Fq  and  F*o.  The  average  time  difference  between  the 
forecast  curves  crossing  the  Fb  threshold  and  F*o  was  12  hours,  and  the  average  time 
difference  between  the  forecast  vorticity  curves  crossing  the  F*o  vorticity  threshold  and 
F*o  was  -16.2  hours.  These  values  indicate  that  while  the  forecast  vorticity  curves 
exceeded  the  F*o  threshold  roughly  16  hours  prior  to  F*o,  they  did  not  exceed  the  F) 
threshold  until  12  hours  after  F*o.  One  possible  explanation  for  the  delay  between 
analyzed  and  forecast  curves  crossing  the  vorticity  thresholds  is  the  incorrect 
representation  of  the  circulation  in  the  NOGAPS  model,  with  the  correction  occurring 
after  the  synthetic  tropical  cyclone  observations  are  inserted  into  the  model  analysis. 


Table  3.8.  Timing  error  (?  t2)  between  the  time  when  a  forecast  vorticity  line  crosses  the 
Fo  vorticity  threshold  and  when  it  crosses  the  F*o  threshold  in  Figure  3.6.  Average 
timing  error  is  29.4  hours.  A  negative  number  indicates  that  the  forecast  time  curve 
crossed  the  Fq  or  F*o  threshold  prior  to  F*o-  A  positive  number  indicates  that  the  forecast 
time  curve  exceeded  the  threshold  value  after  F*o.  No  F*o  data  value  is  calculated  for  the 
+00  forecast  time  because,  by  definition,  the  threshold  value  is  determined  by  the  value 
of  the  analyzed  curve  (+00)  at  F*o- 


Forecast  Time 

Fo  cross  -  F*o 

F*o  cross  -  F*o 

?t2 

+00 

18 

-- 

18 

+12 

24 

-6 

30 

+24 

18 

-12 

30 

+36 

18 

-6 

24 

+48 

12 

-6 

18 

+60 

12 

6 

6 

+72 

18 

-42 

60 

+84 

30 

-42 

72 

+96 

-12 

-42 

30 

+120 

-18 

-12 

6 

(2)  Subregion  Assessments.  To  assess  whether  the  NOGAPS  model 
produced  more  accurate  forecasts  for  specific  regions  within  the  western  North  Pacific, 
850-mb  relative  vorticity  was  examined  by  subregion.  The  general  trend  of  transitioning 
from  over-  to  under-forecasts  within  6  hours  of  It)  was  not  readily  apparent  in  each 
subregion.  This  may  partially  be  due  to  the  small  sample  size  (a  total  23  vortices  in  the 
three  subregions). 
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The  trend  for  the  NOGAPS  model  transition  from  over-  to  under- forecasts  of  850- 
mb  relative  vorticity  apparent  in  Figure  3.5  is  not  evident  in  Figure  3.8,  which  represents 
only  those  vortices  that  developed  within  the  South  China  Sea.  This  is  likely  due  to  the 
small  sample  size  (five  circulations)  and  limited  number  available  for  analysis.  To  plot 
the  forecast  curves  in  Figure  3.8,  the  minimum  sample  size  was  decreased  from  five  cases 
per  time  interval  to  one  case.  A  NOGAPS  analysis  of  850- mb  relative  vorticity  was  not 
available  at  -120  h,  and  model  relative  vorticity  forecasts  were  not  available  prior  to  -60 
h. 

From  60  h  prior  to  Fb  to  48  h  after  Fd,  the  analyzed  relative  vorticity  curve  in 
Figure  3.8  has  a  gradual  increase  as  expected  for  slowly  developing  tropical  cyclones. 
The  notable  decrease  in  vorticity  after  -i-60  h  reflects  the  relatively  short  duration  of  SCS 
vortices  compared  to  other  subregions.  The  average  vorticity  value  at  Fq  for  developing 
vortices  in  the  SCS  subregion  is  4.8  x  10“^  s‘\  which  is  similar  to  the  value  (5.0  x  10'^  s'^) 
for  all  developing  vortices. 

In  general,  850- mb  relative  vorticity  was  under- forecast  by  NOGAPS  at  all  times 
for  vortices  that  developed  within  the  SCS.  The  only  forecast  curves  that  exceeded  the 
analyzed  vorticity  curve  were  -i-72  h  and  -i-96  h,  and  for  forecasts  made  prior  to  It). 
Shorter- range  forecasts  (<  60  h),  while  under- forecast  at  each  temporal  increment,  were 
more  accurate  than  the  longer-range  forecasts.  The  shorter-range  forecasts  more 
accurately  depicted  the  slower  increase  and  decrease  in  relative  vorticity  that  is  in  the 
analyzed  vorticity  curve.  No  statistically  significant  errors  (no  dark  shaded  boxes  as  in 
Figure  3.5)  occurred  until  48  hours  after  Fq,  and  the  errors  that  did  occur  were  for  relative 
vorticity  less  than  analyzed. 

Average  analyzed  and  forecast  850- mb  relative  vorticity  values  for  vortices  that 
developed  within  the  Philippine  Sea  subregion  are  shown  in  Figure  3.9.  As  in  Figure  3.5, 
the  general  trend  of  the  NOGAPS  model  forecasts  is  for  increasing  relative  vorticity 
values  with  time,  which  indicates  that  development  occurred  relatively  steadily,  and  the 
model  forecasts  of  vorticity  did  not  deviate  substantially  from  the  average  analyzed 
values.  While  the  forecast  ?  curves  follow  the  same  trend  of  increasing  with  time,  there 
is  no  distinctive  time  at  which  forecast  curves  transition  from  over-  to  under- forecast  as 
in  Figure  3.5.  The  -i-24  h  NOGAPS  relative  vorticity  forecast  curve  (upper- left  panel  of 
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Figure  3.8.  Average  analyzed  and  forecast  850- mb  relative  vorticity  as  in  Figure  3.5, 
except  for  vortices  that  developed  within  the  SCS  subregion.  Heavy  black  line  represents 
analyzed  average  850- mb  relative  vorticity  (at  -i-OO),  and  red  line  represents  average 
forecast  850-mb  relative  vorticity  for  the  indicated  forecast  time.  Average  analyzed 
vorticity  relative  to  Fq  for  developing  vortices  in  the  SCS  (light  dotted  black  line)  was  4.8 
X  10'^  s'*.  The  light  solid  black  line  represents  the  ?o  threshold  (5.0  x  10'^  s'*)  determined 
in  Figure  3.5. 


Figure  3.9)  begins  to  increase  72  h  prior  to  Fo  and  exceeds  the  ?o  threshold  determined  for 
developing  vortices  in  all  subregions  (5.0  x  10'^  s'*)  about  24  h  prior  to  F).  The  -i-24  h 
forecast  curve  reaches  a  local  maximum  value  at  F)  and  then  decreases  in  the  12  h 
immediately  following  F)  before  beginning  to  steadily  increase  again  after  +12  h.  The 
increase  prior  to  Fo  is  associated  with  NOGAPS  developing  the  vortex  prior  to  the 
insertion  of  the  synthetic  tropical  cyclone  observation.  The  decrease  after  Fq  and  the 
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Figure  3.9.  Average  analyzed  and  forecast  850- mb  relative  vorticity  as  in  Figure  3.5, 
except  for  vortices  that  developed  within  the  PS  subregion.  Heavy  black  line  represents 
analyzed  average  850- mb  relative  vorticity  (at  -tOO),  and  red  line  represents  average 
forecast  850-mb  relative  vorticity  for  the  indicated  forecast  time.  Average  analyzed 
vorticity  relative  to  Fq  for  developing  vortices  in  the  PS  (light  dotted  black  line)  was  5.5  x 
10'^  s Light  solid  black  line  represents  the  ?o  threshold  (5.0  x  10'^  s'^)  determined  in 
Figure  3.5. 


subsequent  increase  in  vorticity  following  Fo  is  attributed  to  the  more  accurate 
representation  of  the  relative  vorticity  after  the  synthetic  vortex  was  inserted. 

The  average  vorticity  for  vortices  developing  within  the  PS  region  was  5.5  x  10"^  s  \ 
which  is  slightly  higher  than  the  average  for  all  developing  vortices.  This  may  be  due  to 
the  higher  background  cyclonic  vorticity  associated  with  the  monsoon  trough. 

The  analyzed  and  forecast  850-mb  relative  vorticity  curves  for  the  EMT  region 
(Figure  3.10)  most  closely  resemble  the  relative  vorticity  curves  for  all  vortices  in  Figure 
3.5.  That  is,  these  storms  have  a  general  trend  towards  increasing  relative  vorticity 
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values  with  increasing  time.  As  in  Figure  3.5,  the  NOGAPS  model  tends  to  over- forecast 
relative  vorticity  in  this  subregion  prior  to  Fq,  and  under- forecast  it  after  Fq.  The  average 
analyzed  850-mb  relative  vorticity  for  vortices  developing  in  the  EMT  subregion  at  Fq  is 
4.7  X  10'^  s'\  which  is  not  significantly  different  from  ?o  based  on  all  storms.  Almost  all 
vorticity  forecasts  prior  to  Fq  are  greater  than  the  analyzed  vorticity  at  each  forecast  time 
interval.  The  model  transition  to  under- forecasts  of  vorticity  occurs  within  12  hours  of 
Fo,  as  expected  with  the  insertion  of  the  synthetic  tropical  cyclone  observations. 
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Figure  3.10.  Average  analyzed  and  forecast  850-mb  relative  vorticity  as  in  Figure  3.5, 
except  for  vortices  that  developed  within  the  EMT  subregion.  Heavy  black  line 
represents  analyzed  average  850-mb  relative  vorticity  (at  -rOO),  and  red  line  represents 
average  forecast  850-mb  relative  vorticity  for  the  indicated  forecast  time.  Average 
analyzed  vorticity  relative  to  Fb  for  developing  vortices  in  the  EMT  (light  dotted  black 
line)  was  4.7  x  10'^  s Eight  solid  black  line  represents  the  ?o  threshold  (5.0  x  10'^  s'^) 
determined  in  Eigure  3.5. 
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Interestingly,  the  NOGAPS  model  predicts  maximum  intensity  too  early  in  this 
basin  for  short-range  forecasts.  Maximum  relative  ■\orticity  is  forecast  earlier  in  the  +\2 
h  forecasts  than  in  the  -i-36  h  forecasts  (not  shown).  The  -i-48  h  forecasts  and  +12  h 
forecasts  (Figure  3.10)  were  more  accurate  on  the  timing  and  magnitude  of  maximum 
relative  vorticity.  This  deficiency  in  the  shorter-range  forecasts  may  indicate  a  tendency 
of  the  NOGAPS  model  to  intensify  developing  vortices  in  the  EMT  subregion  too 
quickly. 

(3)  Assessment  Relative  to  the  TCFAs.  To  examine  whether  model 
guidance  accuracy  may  be  related  to  the  accuracy  of  the  placement  of  the  TCFA  forecast 
window  relative  to  the  actual  formation  time,  relative  vorticity  was  examined  based  on 
formation  time  relative  to  the  TCFA-defined  formation  window  for  each  JTWC- 
numbered  circulation.  By  identifying  formation  time  as  the  time  of  he  first  JTWC 
warning,  the  TCFA  formation  windows  for  developing  circulations  were  verified.  A 
positive  verification  occurred  if  the  first  JTWC  warning  was  issued  during  the  formation 
window  defined  in  the  TCFA  for  that  storm. 

Eleven  of  27  (40.7%)  TCFAs  verified  prior  to  the  onset  of  the  formation  window. 
The  average  analyzed  and  forecast  850- mb  vorticity  curves  for  those  11  circulations  that 
verified  prior  to  the  TCFA  formation  window  are  plotted  in  Figure  3.11.  Prior  to  Fo,  the 
analyzed  and  forecast  curves  for  this  subset  of  vortices  are  very  similar  to  Figure  3.5, 
which  contains  curves  for  all  developing  vortices.  However,  the  amount  of  under¬ 
forecasting  following  Fb  is  much  larger  for  this  subset  of  vortices  than  in  Figure  3.5  as 
shown  by  the  amount  of  deviation  between  the  analyzed  and  forecast  ?  curves.  The 
change  in  slope  of  the  analyzed  (-i-OO)  vorticity  curve  in  Figure  3.11  is  more  dramatic  for 
this  subset  of  developers,  especially  in  the  shorter-range  forecasts.  The  combination  of 
an  early  formation  time  relative  to  the  TCFA  formation  window,  the  dramatic  increase  in 
slope  and  severe  under- forecasting  errors  after  It)  suggest  that  these  vortices  may  have 
undergone  a  period  of  rapid  intensification. 

The  first  warning  was  issued  during  the  TCFA-defined  formation  window  for 
34.8%  (8)  of  the  23  developing  vortices.  Analyzed  and  forecast  850- mb  relative  vorticity 
for  these  circulations  are  shown  in  Figure  3.12.  Given  this  relatively  small  sample,  the 
minimum  number  of  cases  has  to  be  relaxed  in  this  comparison  (see  samples  sizes  near 
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Figure  3.11.  Average  analyzed  and  forecast  850- mb  relative  vorticity  displays  as  in 
Figure  3.5,  except  for  11  developing  vortices  that  were  warned  on  PRIOR  to  the  TCFA- 
defined  formation  window.  Heavy  black  line  represents  analyzed  average  850-mb 
relative  vorticity  (at  +00),  and  red  line  represents  average  forecast  850-mb  relative 
vorticity  for  the  indicated  forecast  time.  Average  analyzed  vorticity  relative  to  F)  for 
vortices  verifying  prior  to  the  TCFA  formation  window  (light  dotted  black  line)  was  4.9  x 
10'^  s'\  and  5.0  x  10'^  s‘*  for  all  vortices  (light  solid  black  line). 


the  tops  of  the  panels),  and  thus  the  differences  should  be  regarded  as  tentative.  For  this 
subset  of  developers,  the  forecast  vorticity  characteristics  prior  to  Fq  are  similar  to  the 
total  set  of  all  developing  vortices  in  Figure  3.5.  However,  the  transition  from  over-  to 
under- forecasts  begins  as  early  as  12  hours  prior  to  Fq,  and  continues  until  36  hours  after 
Fo  (not  shown).  The  absence  of  a  drastic  slope  change  in  the  24  hours  surrounding  F) 
indicates  that  rapid  intensification  did  not  occur  for  this  subset,  and  suggests  that  the 
NOGAPS  model  may  perform  better  during  the  formation  and  intensification  process  for 
this  subset. 
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Figure  3.12.  Average  analyzed  and  forecast  850- mb  relative  vorticity  displays  as  in 
Figure  3.5,  except  for  developing  vortices  that  were  warned  on  DURING  the  TCFA- 
defined  formation  window.  Heavy  black  line  represents  analyzed  average  850-mb 
relative  vorticity  (at  -I-OO),  and  red  line  represents  average  forecast  850-mb  relative 
vorticity  for  the  indicated  forecast  time.  Average  analyzed  vorticity  relative  to  Fb  for 
vortices  verifying  during  the  TCFA  formation  window  (light  dotted  black  line)  was  4.9  x 
10'^  s'^  and  5.0  x  10'^  s'^  for  all  developing  vortices  (light  solid  black  line). 


Five  of  the  24  (20.8%)  TCFAs  verified  after  the  specified  TCFA  formation 
window.  All  five  of  those  TCFAs  were  reissued  and  subsequently  four  verified  during 
the  specified  TCFA  formation  window  and  one  verified  before  the  specified  formation 
window.  These  reissued  TCFAs  were  included  in  the  subset  discussed  above,  while  the 
original  TCFAs  that  were  changed  are  included  in  the  subset  of  vortices  that  verified  after 
the  TCFA  formation  window.  The  average  analyzed  and  forecast  850-mb  relative 
vorticity  for  these  five  circulations  are  included  in  Figure  3.13.  To  show  the  forecast 
curves,  the  minimum  number  of  cases  displayed  at  each  time  interval  was  decreased  to 
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Figure  3.13.  Average  analyzed  and  forecast  850- mb  relative  vorticity  displays  as  in 
Figure  3.5,  except  for  developing  vortices  that  were  warned  on  AFTER  the  TCFA- 
defined  formation  window.  Heavy  black  line  represents  analyzed  average  850-mb 
relative  vorticity  (at  +00),  and  red  line  represents  average  forecast  850-mb  relative 
vorticity  for  the  indicated  forecast  time.  Average  analyzed  vorticity  relative  to  E)  for 
vortices  verifying  after  to  the  TCFA  formation  window  (light  dotted  black  line)  was  6.3  x 
10'^  s'^  and  5.0  x  10'^  s'^  for  all  developing  vortices  (light  solid  black  line). 


one.  The  over-  and  under-forecast  trends  in  850-mb  vorticity  relative  to  the  Fo  threshold 
evident  in  Figure  3.5  are  not  evident  in  Figure  3.13.  The  transition  between  over-  and 
under- forecasts  of  relative  vorticity  instead  occurs  for  short-range  forecasts  between  24 
and  36  h  after  Fq.  The  lack  of  a  change  in  the  slope  of  the  analysis  curve  at  Fq,  and  the 
lack  of  the  predominant  under- forecasting  that  was  evident  with  respect  to  early 
formations  (Figure  3.11),  and  to  a  lesser  degree  with  the  on-time  formations  (Figure 
3.12),  suggest  that  these  cases  developed  relatively  slowly. 
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Also  of  interest  is  the  relatively  high  average  analyzed  850- mb  relative  vorticity 
value  at  Fo  for  this  subset  of  developers.  The  average  analyzed  vorticity  for  circulations 
verifying  after  the  TCFA- specified  formation  window  was  6.3  x  10'^  s'\  which  is 
noticeably  higher  than  the  average  analyzed  vorticity  value  at  Fo  for  the  total  set  of 
developers.  This  number  is  likely  erroneously  high  due  to  the  small  sample  size  of  this 
subset.  However,  the  higher  value  could  also  indicate  that  the  circulation  is  more  intense 
at  the  time  of  the  first  warning. 

b.  Deep  Layer  Wind  Shear 

While  the  deep  layer  wind  shear  (200-850  mb)  variable  provides  insight  about  the 
relative  magnitudes  of  the  upper- level  anticyclone  and  low-level  cyclone  associated  with 
the  developing  tropical  circulation,  it  primarily  indicates  the  prevailing  direction  of  the 
environmental  flow  surrounding  the  tropical  circulation.  A  positive  shear  value  indicates 
westerly  shear,  or  an  increase  in  westerly  winds  with  height.  A  negative  shear  value 
indicates  easterly  shear,  or  in  increase  in  easterly  winds  with  height.  While  classic 
research  on  tropical  cyclone  formation  by  Gray  (1968)  indicates  that  minimum  shear  over 
the  developing  circulation  is  a  necessary  condition  for  formation,  more  recent  research 
(i.e.,  McBride  and  Zehr  1981)  points  to  the  importance  of  a  transition  zone,  with  strong 
westerly  shear  to  the  north  and  strong  easterly  shear  to  the  south  of  the  developing 
tropical  circulation.  While  the  deep  layer  wind  shear  variable,  as  measured  here,  does  not 
have  a  horizontal  gradient  of  wind  shear,  it  contains  valuable  information  about  the 
atmosphere  in  which  the  tropical  circulation  is  developing.  Kurihara  and  Tuleya  (1981) 
found  in  a  Global  Fluid  Dynamics  Laboratory  (GFDL)  study  that  easterly  shear  is  more 
favorable  than  westerly  shear  for  tropical  cyclone  formation,  which  matches  the 
westward  translation  speed  of  the  wave  with  the  deep- layer  average  easterly  winds  if 
there  is  easterly  shear  through  the  troposphere.  Consequently,  a  negative  (easterly)  deep 
shear  favors  formation. 

For  this  set  of  developing  vortices  (Figure  3.14),  the  analyzed  deep  layer  wind 
shear  transitions  from  positive  (westerly)  to  negative  (easterly)  between  48  and  36  h  prior 
to  Fo.  Interestingly,  the  shear  remains  just  slightly  negative  (easterly)  until  F),  when  it 
begins  to  steadily  decrease  (become  more  negative)  as  the  circulation  intensifies.  The 
average  analyzed  deep  shear  at  Fo  was  -1.8  m  s‘\  which  indicates  easterly  shear  over  the 
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Figure  3.14.  Analyzed  and  forecast  deep  wind  shear  (200-850  mb)  for  developing 
vortices  relative  to  Fb,  as  in  Figure  3.5.  Heavy  black  line  represents  analyzed  average 
deep  wind  shear  (at  +00),  and  red  line  represents  average  forecast  deep  wind  shear  for  the 
indicated  forecast  time.  Average  analyzed  deep  wind  shear  at  Fo  for  developing  vortices 
(light  solid  black  line)  was  -1.8  m  s 


circulation  center,  which  is  favorable  for  the  development  of  the  tropical  circulation. 
Since  the  deep  shear  variable  can  be  either  positive  or  negative  (indicating  directional 
change)  the  terms  over- and  under- forecast  must  be  more  specifically  defined.  An  under¬ 
forecast  of  deep  shear  is  then  defined  to  not  be  a  negative  value,  but  a  forecast  value  that 
is  weaker  or  closer  to  zero  than  the  analyzed  value.  Likewise,  an  over- forecast  of  deep 
shear  is  defined  here  as  a  forecast  value  that  is  stronger  or  greater  than  the  analyzed 
value.  As  shown  in  Figure  3.14,  the  forecast  deep  layer  shear  curves  transition  from  less 
than  the  analyzed  shear  to  greater  than  the  analyzed  shear  approximately  12  h  prior  to  the 
insertion  of  the  synthetic  tropical  cyclone  observations. 
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The  forecast  deep  wind  shear  prior  to  Fq  is  too  low  due  to  excessive  easterly  wind 
shear  associated  with  the  NOGAPS  model  over- forecasts  of  850- mb  relative  vorticity 
during  the  same  period.  At  each  forecast  time,  a  slight  increase  in  the  deep  shear  variable 
occurs  just  prior  to  the  insertion  of  the  synthetic  tropical  cyclone  observation.  This 
increase  in  westerly  shear  is  associated  with  the  NOGAPS  model  over- forecasting  the 
850-mb  relative  vorticity,  and  thus  strengthening  easterly  environmental  flow.  This 
coupling  of  vorticity  and  wind  shear  is  expected  in  a  hydrostatically  and  geostrophically 
balanced  model. 

Throughout  the  entire  forecast  period,  the  differences  between  the  analyzed  and 
forecast  deep  wind  shear  are  small.  This  lack  of  significant  variation  is  due  to  the  already 
weak  shear  associated  with  the  developing  tropical  circulation.  Forecasts  of  vertical 
shear  away  from  the  center  of  the  ellipse  would  likely  show  greater  deviations  from  the 
analyzed  values. 

c.  Sea-level  Pressure 

Sea- level  pressure  (SLP)  varies  only  minimally  in  the  tropics  (McBride  1981). 
However,  examining  the  SLP  variable  provides  valuable  information  about  the  intensity 
of  the  developing  circulation.  As  expected,  the  average  analyzed  and  forecast  sea-  level 
pressure  (SLP)  decreases  with  increasing  forecast  time  for  these  developing  tropical 
cyclones  (Figure  3.15).  The  mean  analyzed  SLP  at  Fq  (Pq)  was  lower  (1006.9  mb)  than 
the  1007.4  mb  at  F*o  (P*o)  as  expected  from  the  deepening  of  the  storm  between  these 
two  times.  Although  the  forecast  SLP  values  were  only  slightly  less  than  the  analyzed 
SLP  (Figure  3.15),  this  difference  is  important.  Only  small  SLP  drops  occur  during  the 
early  stages  of  tropical  cyclone  development  when  the  warm  core  is  relatively  weak.  The 
trend  of  NOGAPS  SLP  forecasts  to  be  slightly  lower  than  the  average  analyzed  value  is 
consistent  with  the  model  forecasts  of  850-mb  relative  vorticity.  This  tendency  is 
expected  as  the  NOGAPS  model  fields  should  be  in  hydrostatic  and  geostrophic  balance 
during  the  early  stages  of  tropical  cyclone  development,  and  at  these  horizontal  scales. 

One  significant  trend  in  the  NOGAPS  model  forecasts  of  SLP  is  the 
disproportionate  number  of  forecasts  that  are  too  low  when  compared  to  the  analyzed 
SLP  values.  All  of  the  extreme  values  (indicated  in  Figure  3.15  by  an  open  circle)  err  on 
the  low  side  of  both  the  analyzed  and  forecast  curves.  Additionally,  most  of  the  extreme 
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Figure  3.15.  Average  analyzed  and  forecast  sea-level  pressure  (mb)  for  all  vortices 
relative  to  Fb,  as  in  Figure  3.5.  Heavy  black  line  represents  analyzed  average  SLP  (at 
+00),  and  red  line  represents  average  forecast  SLP  for  the  indicated  forecast  time. 
Average  analyzed  SLP  at  Fo  (light  solid  black  line)  was  1006.9  mb. 


SLP  values  are  forecast  prior  to  Fo,  with  the  exception  of  the  +24  h  forecasts.  While  the 
short-term  (24  h)  average  forecast  SLP  is  consistent  with  the  analyzed  SLP,  the  NOGAPS 
model  tendency  to  err  towards  excessive  and  rapid  development  is  evident  in  Figure  3.15. 
d.  925-mb  Wind  Speed 

The  925-mb  wind  speed  is  another  indicator  of  the  intensity  of  the  developing 
tropical  circulation.  This  variable  is  closely  linked  to  the  850-mb  relative  vorticity.  As 
expected,  a  transition  from  over-  to  under- forecasts  of  925-mb  wind  speed  occurs  in  the 
first  12  hours  after  Fq  (Figure  3.16).  The  sudden  increase  in  the  slope  of  the  forecast 


46 


30 


Fq  2002  WPAC  NOGAPS  DV  al  +24 


Fq  2002  WPAC  NOGAPS  DV  at  +48 


4 111214121618202122232122221918161913 


0 - ^ - ' . . 

-96  -72  -48  -24  00  +24  +48  +72  +96 

Forecast  Verification  Time 


-96  -72  -48  -24  00  +24  +48  +72  *96 


Forecast  Verification  Time 


Forecast  Verification  Time 


-96  -72  -48  -24  00  +24  *48  *72  +96 


Forecast  Verification  Time 


Figure  3.16.  Average  analyzed  and  forecast  925-mb  wind  speed  (m  s‘^)  for  all  vortices 
relative  to  Fb,  as  in  Figure  3.5.  Heavy  black  line  represents  analyzed  average  925-mb 
wind  speed  (at  -i-OO),  and  red  line  represents  average  forecast  925-mb  wind  speed  for  the 
indicated  forecast  time.  Average  analyzed  925-mb  wind  speed  at  F)  (light  solid  black 
line)  was  10.8  m  s 


wind  speed  curve  in  each  panel  of  Figure  3.16  indicates  the  NOGAPS  model  reaction  to 
the  insertion  of  the  synthetic  tropical  cyclone  observations  into  the  model  analysis.  In 
general,  the  deviation  of  the  forecast  wind  speed  from  the  analyzed  wind  speed  increases 
with  time,  which  indicates  the  continued  intensificatbn  of  the  developing  tropical 
circulation. 


e.  Vapor  Pressure 

The  vapor  pressure  averaged  over  500-700  mb  (Pa)  provides  an  indication  of  the 
amount  of  mid- level  moisture  available  to  the  developing  tropical  cyclone.  The  slope  of 
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the  analyzed  (at  +00)  vapor  pressure  (Figure  3.17)  remains  relatively  flat  throughout  the 
development  of  the  tropical  circulation.  This  lack  of  change  in  the  slope  of  the  analyzed 
curve  makes  it  a  poor  variable  for  use  in  the  determination  of  a  formation  threshold. 

At  almost  all  forecast  times  in  Figure  3.17  the  NOGAPS  model  forecasts  of  vapor 
pressure  are  below  the  analyzed  value.  Another  interesting  factor  is  the  set  of  extreme 
forecast  values,  which  are  indicated  in  Figure  3.17  by  an  open  circle.  These  extreme 
values  are  almost  aU  substantially  below  even  the  forecast  vapor  pressure  (indicating  a 
drier  mid- troposphere  than  analyzed)  and  they  tend  to  occur  at  extended  forecast  ranges. 
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Figure  3.17.  Analyzed  and  forecast  vapor  pressure  (500-700  mb  average)  (Pa)  for  all 
vortices  relative  to  Fb,  as  in  Figure  3.5.  Heavy  black  line  represents  analyzed  average 
vapor  pressure  (at  +00),  and  red  line  represents  average  forecast  vapor  pressure  for  the 
indicated  forecast  time.  Average  analyzed  vapor  pressure  at  Fb  (light  solid  black  line) 
was  3.6  Pa. 
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Because  of  this  tendency  for  the  NOGAPS  model  to  incorrectly  forecast  this  variable, 
mid- level  vapor  pressure  is  not  likely  to  be  a  good  one  for  distinguishing  developing  from 
non- developing  vortices. 

Additionally,  the  consistent  under- forecasting,  evident  in  both  the  forecast  curves 
and  also  the  outliers  of  the  forecast  vapor  pressure  indicates  that  the  NOGAPS  model  is 
not  correctly  assessing  the  interior  vortex  processes.  This  excessive  dryness  at  the  mid¬ 
levels  is  possibly  a  result  of  the  convective  parameterization.  As  noted  in  Elsberry 
(2003),  a  dry  mid-troposphere  inhibits  formation  by  preventing  the  development  of  a 
warm  core.  However,  within  the  NOGAPS  model,  this  excessive  dryness  in  the  mid¬ 
troposphere  may  actually  limit  the  number  of  tropical  cyclone  formations,  both  accurate 
and  inaccurate,  that  the  NOGAPS  model  develops.  Formations  then  are  a  result  of  large- 
scale  processes  rather  than  the  internal  dynamics  of  the  tropical  vortex. 

/.  Other  Variables 

Seven  variables  (not  shown)  listed  in  Table  2.1  were  examined  to  determine  the 
presence  of  a  forecast  signal,  but  were  not  selected  for  further  analysis.  Of  these 
variables,  most  showed  little  to  no  slope  of  the  analyzed  variable  curve,  and  also  had 
average  forecast  curves  that  were  consistently  lower  than  the  analyzed  curve,  and  a 
distinct  trend  toward  excessively  low  forecast  values.  This  flat  slope,  which  is  consistent 
with  under-forecasting  of  developing  vortices,  and  the  lack  of  deviations  from  the 
analyzed  values  may  provide  value  information  regarding  the  atmosphere  in  which  the 
tropical  cyclone  develops,  or  information  regarding  the  developing  circulation. 
Moreover,  many  of  these  variables  were  model-derived  rather  than  directly  analyzed 
variables  and  thus  more  likely  contained  model  errors. 

Vertical  motion  was  consistently  too  large  and  positive  after  F).  This  error  is 
likely  due  to  the  NOGAPS  over- forecasts  of  precipitation.  Total  (convective  plus  grid- 
scale)  precipitation  in  the  storm- centered  ellipse  was  consistently  over- forecast  for  the 
developing  vortices  at  all  forecast  times  by  NOGAPS.  Excessive  precipitation 
forecasting  is  attributed  to  the  cumulus  parameterization  within  the  model. 

Eatent  heat  flux  from  the  ocean  was  excluded  from  the  selected  variable  list 
because  it  is  a  purely  model-derived  variable,  and  its  accuracy  is  subject  to  the  correct 
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representation  of  the  wind  speed  and  sea- surface  temperature.  It  may  also  not  be  an 
important  physical  process  at  the  early  stages  of  formation. 

Geopotential  height  thickness  (200  -  1000  mb)  had  no  strong  signal  for  the 
developing  vortices.  The  forecast  values  were  consistently  greater  than  analyzed  for 
these  developing  vortices,  which  might  be  expected  solely  due  to  the  increase  in  relative 
vorticity  and  falling  1000- mb  heights  associated  with  the  intensification  of  the  tropical 
circulation.  Likewise,  the  shallow  (1000  -  500  mb)  warm  anomaly  variable  showed  no 
distinct  signal  for  distinguishing  these  developing  vortices.  The  lack  of  a  lower 
tropospheric  signal  may  be  because  the  warm  core  first  forms  in  the  upper  troposphere  of 
developing  storms. 

The  700- mb  wind  speed  forecasts  did  not  show  a  significant  difference  from  the 
analyzed  curve  prior  to  Fq,  and  the  wind  speed  was  consistently  under- forecast  after  Fq. 
The  under- forecasting  of  this  variable  did  not  become  statistically  significant  until  after 
-1-84  h.  This  variable,  and  the  500-mb  wind  speed  variable,  were  excluded  from  analysis 
because,  in  addition  to  showing  little  prognostic  value  in  the  cases  analyzed,  they  are 
indirectly  contained  in  the  deep  layer  (200  -  850  mb)  wind  shear  variable. 

Shallow  layer  (500  -  850  mb)  wind  shear  was  also  excluded  for  lack  of  a 
prognostic  signal  for  distinguishing  developing  vortices  during  1  May  -  31  October  2002. 
Whether  this  lack  of  a  signal  was  indicative  that  low-level  shear  was  not  an  important 
physical  process  in  tropical  cyclone  formation,  or  a  reflection  that  this  variable  is  not 
well- predicted  by  the  NOGAPS  model,  can  not  be  determined  from  this  sample. 

2.  Non-  developing  Vo  rtices 

The  threshold  values  for  six  forecast  variables  determined  for  the  tropical  vortices 
that  developed  into  systems  that  were  later  numbered  and  warned  on  by  JTWC  will  serve 
as  a  reference  for  the  non- developing  vortices.  Three  sets  of  non- developing  vortices 
identified  by  the  TCVTP  algorithm  will  be  compared  to  the  thresholds  determined  for 
developing  vortices.  These  datasets  (see  definitions  in  Chapter  UFA)  were:  non¬ 
developing  vortices  that  met  the  minimum  duration  criterion  (NDl);  non- developing 
vortices  that  did  not  meet  the  minimum  duration  criterion,  but  did  meet  the  over- forecast 
criteria  (ND2);  and  both  model  and  TCFA  false  alarms. 


50 


Since  no  JTWC  warnings  exist  for  vortices  that  do  not  develop  into  tropical 
cyclones,  the  time  (No)  corresponding  to  the  maximum  analyzed  850- mb  relative 
vorticity  was  used  as  the  reference  time  for  the  non- developing  vortices.  This  time  was 
determined  during  post-storm  analysis.  The  non- developing  vortex  dataset  NDl  included 
104  vortices  that  persisted  for  at  least  24  consecutive  hours.  The  NDl  dataset  was  further 
analyzed  to  identify  model  false  alarms.  The  ND2  dataset  that  included  non- developing 
vortices  that  did  not  meet  the  minimum  duration  criterion  was  analyzed  to  determine 
which  vortices  were  model  over- forecasts.  Only  60  of  127  vortices  in  the  ND2  dataset 
met  all  four  model  over-  forecast  criteria. 

a.  Non-developing  Vortices  (NDl) 

The  next  data  set  considered  consists  of  the  non- developing  vortices  that  met  the 
minimum  duration  criteria  of  24  hours  (NDl).  Of  the  254  tracked  circulations,  104 
vortices  met  the  minimum  duration  criterion  (column  2  of  Table  3.1).  The  average 
duration  of  NDl  vortices  was  3.1  days. 

(1)  8 50- mb  relative  vorticity. 

(a)  Entire  Basin  Assessment.  In  contrast  to  the  analyzed  850- mb 
relative  vorticity  for  the  developing  vortices  (Figure  3.5),  the  analyzed  850- mb  relative 
vorticity  for  NDl  vortices  (Figure  3.18)  has  very  little  slope  change  throughout  the 
lifetime  of  the  potential  tropical  circulation.  The  average  vorticity  at  bb  represents  a 
local  maximum  in  the  analyzed  vorticity  curve.  The  average  analyzed  relative  vorticity  at 
No  (3.3  X  10'^  s'^)  is  less  than  the  average  analyzed  relative  vorticity  at  Fq  (5.0  x  10'^  s'^) 
for  developing  vortices.  The  small  standard  deviations  about  the  analyzed  means  indicate 
that  overall  these  non- developers  have  significantly  smaller  850- mb  relative  vorticity 
than  the  developing  cyclones,  which  is  favorable  for  distinguishing  the  non- developing 
vortices  if  the  forecasts  are  accurate  in  amplitude  and  trend. 

Prior  to  No  the  analyzed  vorticity  increases  only  slowly.  The  relatively  flat  slope 
of  the  analyzed  vorticity  is  attributed  to  the  near  steady  states,  or  at  most  slow 
developments,  of  these  tropical  vortices  or  monsoon  depressions.  The  slope  of  the 
analyzed  vorticity  in  Figure  3.18  decreases  in  the  36  h  following  N),  which  indicates  a 
broadening  or  weakening  of  the  circulation.  Although  the  analyzed  vorticity  value  at  No 
appears  to  be  again  reached  96  h  after  No,  the  circulation  size  is  fairly  small  by  this  time. 
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Figure  3.18.  Average  analyzed  and  forecast  850- mb  relative  vorticity  as  in  Figure  3.5, 
except  for  non- developing  vortices  relative  to  Nq.  Heavy  black  line  represents  analyzed 
average  850-mb  relative  vorticity  for  non- developing  vortices  (at  -i-OO),  and  red  line 
represents  average  forecast  850-mb  relative  vorticity  for  the  indicated  forecast  time. 
Average  analyzed  850-mb  relative  vorticity  at  No  (light  dashed  black  line)  was  3.3  x  10'^ 
s The  ?o  (light  solid  black  line)  of  5.0  x  10'^  s'*  is  for  the  developing  vortices  in  Figure 
3.5  and  provides  a  reference. 


At  forecast  intervals  of  72  h  and  less,  the  mean  forecast  850-mb  relative  vorticity 
curves  in  Figure  3.18  do  not  exceed  the  ?o  value  for  developing  storms  at  any  forecast 
verification  time.  Therefore,  in  the  mean  the  NOGAPS  model  correctly  forecasts  the 
non- development  of  these  circulations.  However,  a  few  outliers  at  -24  h  and  -i-48  h 
(relative  to  N))  do  exceed  the  5.0  x  10'^  s'*  ?o  value  and  thus  could  be  mistaken  for 
developing  storms.  In  the  48-h  forecasts  in  the  upper-right  panel  of  Figure  3.18,  the 
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average  850- mb  vorticity  forecast  values  are  significantly  larger  than  the  analyzed  values 
from  -1-12  h  to  -i-48  h  (darkened  boxes). 

Although  these  average  forecast  values  do  not  exceed  the  5.0  x  10"^  s‘^  ?o  value, 
the  upper  range  and  the  outlier  values  do  exceed  the  ‘^)  value,  which  again  indicates 
erroneous  48- h  formation  forecasts  are  contained  in  this  non- developing  sample.  This 
trend  toward  having  a  range  of  significantly  larger  forecast  vortices  than  analyzed 
vortices  continues  for  the  72-h  and  96-h  forecasts  (lower-left  and  lower-right  panels  in 
Figure  3.18,  respectively).  Although  the  mean  72-h  forecast  vortices  still  do  not  exceed 
the  5.0  X  10'^  s'^  ?o  threshold,  the  box  representing  the  upper  25%  above  the  mean  does 
exceed  this  threshold  for  the  forecasts  at  the  -i-36,  -i-48,  and  -i-60  h  verification  times.  For 
the  96-h  forecast,  even  the  mean  exceeds  the  5.0  x  10'^  s'^  ?o  threshold  at  the  -i-60  h 
through  -1-84  h  verification  times. 

Interestingly,  the  time  at  which  the  maximum  relative  vorticity  is  forecast  to  occur 
increases  with  forecast  range.  In  the  upper-left  panel  of  Figure  3.18,  the  maximum 
vorticity  is  forecast  to  occur  24  h  after  Nq.  This  coincidence  of  a  24- h  forecast  reaching 
maximum  vorticity  24  h  after  the  analyzed  circulation  reached  a  maximum  value  suggests 
the  NOGAPS  model  has  somehow  continued  the  upward  trend  in  analyzed  values  that 
occurred  over  the  previous  24  h.  In  the  48-h  forecasts  (upper-right  panel),  the  forecast 
maximum  vorticity  tends  to  occur  36  h  after  Nq.  Similarly,  the  72-h  forecasts  (lower- left 
panel)  have  maximum  vorticity  occurring  at  48  h  after  Nq.  The  trend  is  continued  for  the 
96-h  forecasts  (lower-right  panel),  where  the  maximum  vorticity  is  forecast  to  occur  72  h 
after  No.  As  noted  above,  the  average  magnitudes  of  the  maximum  vorticity  forecasts 
also  increase  dramatically,  along  with  the  statistical  significance  of  those  average  forecast 
errors. 

The  trends  in  Figure  3.18  confirm,  similarly  to  Figure  3.5,  that  in  2002  the 
NOGAPS  model  tends  to  over- forecast  relative  vorticity.  The  growth  of  errors  in  the 
longer  forecasts  suggests  a  systematic  error  that  may  be  attributed  to  an  overly  active 
convection  scheme.  If  one  assumes  these  non- developing  storms  do  not  contain  the 
mesoscale  convective  systems  (and  vortices)  that  are  assumed  to  exist  in  the  developing 
systems,  then  the  NOGAPS  parameterizations  of  convection  and  friction  may  effectively 
simulate  their  effects  in  an  excessive  manner.  An  alternate  explanation  s  that  these 
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convective  uncertainties  are  related  to  the  unpredictability  of  the  tropical  atmosphere  at 
extended  time  ranges. 

To  determine  if  inaccurate  representation  of  the  vortex  was  causing  some  of  the 
forecast  vorticity  errors,  vortex  size  was  examined  (Figure  3.19).  During  the  24  h  before 
and  after  No,  vortex  size  remains  relatively  consistent  around  the  analyzed  No  value  of 
54.6  grid  cells,  which  is  only  slightly  less  than  the  average  analyzed  size  at  Fo  (60  grid 
cells).  Otherwise,  no  distinct  trend  in  vortex  size  was  noted  in  Figure  3.19,  although  in 
the  early  stages  of  the  vortex  life  cycle  a  small  increase  in  vortex  size  may  occur  from  84 
to  48  hours  prior  to  No. 
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Figure  3.19.  Average  analyzed  and  forecast  vortex  size  (grid  cells)  for  all  non¬ 
developing  vortices  relative  to  No  as  in  Figure  3.7.  Heavy  black  line  represents  analyzed 
average  vortex  size  for  non- developing  vortices  (at  -i-OO),  and  red  line  represents  average 
forecast  vortex  size  for  the  indicated  forecast  time.  Average  analyzed  vortex  size  at  No 
(light  dashed  black  line)  was  54.6  grid  cells.  Average  size  at  Fo  was  60  grid  cells. 
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That  very  little  differenee  exists  between  the  size  of  the  developing  vortiees  and 
non- developing  vortices  at  the  first  JTWC  warning  and  maximum  vorticity  time, 
respectively,  implies  that  the  size  differences  between  the  two  circulations  are  small,  at 
least  early  in  life  cycle  of  the  circulation.  At  about  96  h  prior  to  Fq,  when  the  developing 
circulation  (Figure  3.5)  has  the  same  analyzed  relative  vorticity  as  No  (3.3  x  10'^  s‘*),  the 
corresponding  average  size  of  the  developing  circulations  (Figure  3.7)  is  about  53  grid 
cells,  which  is  comparable  to  the  size  of  the  non- developing  vortices  at  Nq.  The  primary 
difference  between  these  two  data  sets  in  the  hours  following  this  equivalence  point  is  the 
steady  increase  in  analyzed  850- mb  relative  vorticity  for  the  developing  circulations,  and 
initial  decrease  in  analyzed  vorticity  for  the  non- developing  vortices.  Since  the  sizes  are 
about  the  same,  the  spin- up  of  the  vorticity  in  the  developing  cyclones  must  be  at  interior 
points. 

Forecasts  of  vortex  size  (Figure  3.19)  were  not  preferentially  over-  or  under¬ 
forecast,  which  indicates  no  strong  trend  in  the  NOGAPS  model  forecasts.  Only  a  weak 
trend  toward  under-forecasting  the  vortex  size  was  noted  in  the  hours  following  No. 

Notice  the  trend  for  the  outliers  to  be  forecast  to  be  too  large  by  a  factor  of  2-3.  An 
undersized  vortex  is  expected  with  over-  forecast  relative  vorticity  in  a  hydrostatically  and 
geostrophically  balanced  model. 

(b)  Subregion  Assessments.  To  determine  if  the  NOGAPS 
model  preferentially  developed  vortices  in  one  of  the  three  western  North  Pacific  Ocean 
subregions,  the  NDl  vortices  were  separated  based  on  formation  location.  There  were 
more  non- developing  vortices  in  the  PS  subregion  (44%)  than  in  the  SCS  (28%)  or  EMT 
(28%)  subregions  (column  two  of  Table  3.1).  Non- developing  vortices  in  the  SCS 
subregion  persisted  slightly  longer  (3.6  days)  than  vortices  in  either  the  EMT  (2.9  days) 
or  PS  (2.7  days)  subregions.  The  monsoon  depressions  in  the  SCS  may  drift  slower  due 
to  opposing  flows  on  either  side  of  the  monsoon  trough  leading  to  a  small  steering  effect. 
In  the  PS,  the  westward  motion  may  be  larger  and  the  circulations  either  develop  or  move 
out  of  the  region.  The  longer  duration  of  non- developing  vortices  in  the  SCS  subregion  is 
possibly  associated  with  an  earlier  identification  of  the  vortex  in  the  NOGAPS  model 
analyses,  due  to  better  observations  or  a  more  organized  circulation  in  the  NOGAPS 
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model  analysis.  It  is  possible  that  these  vortices  took  longer  to  dissipate  once  the 
development  ceased  than  in  other  subregions.  Another  possible  explanation  is  that  the 
non- developing  vortices  in  the  SCS  experienced  several  growth  and  development  pulses 
prior  to  dissipation,  as  might  be  inferred  from  the  repeated  increases  and  decreases  in  the 
slope  of  the  analyzed  vorticity  in  Figure  3.20. 

The  forecasts  of  relative  vorticity  the  SCS  subregion  (Figure  3.20)  do  not  reflect 
any  pulses  in  intensity,  but  instead  are  consistently  greater  than  analyzed,  and  also 
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Figure  3.20.  Average  analyzed  and  forecast  850- mb  relative  vorticity  (x  10'^  s'^)  as  in 
Figure  3.18,  except  for  non- developing  vortices  within  the  SCS  subregion  relative  to  Nq. 
Heavy  black  line  represents  analyzed  average  850-mb  relative  vorticity  for  non¬ 
developing  vortices  (at  -i-OO),  and  red  line  represents  average  forecast  850-mb  relative 
vorticity  for  the  indicated  forecast  time.  Average  analyzed  850-mb  relative  vorticity  for 
this  subset  of  non- developing  vortices  at  No  (light  dashed  black  line)  was  3.7  x  10"^  s ?o 
(light  solid  black  line)  is  5.0  x  10'^  s 
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steadily  increasing.  As  in  the  overall  sample  (Figure  3.18),  the  range  of  forecast  values 
increases  with  increasing  forecast  range,  as  does  the  time  of  maximum  forecast  vorticity. 
However,  the  sample  sizes  quickly  become  small  at  extended  forecast  intervals. 

The  duration  of  the  vortices  in  the  PS  (2.7  days)  was  slightly  less  than  the  average 
duration  for  all  non- developing  vortices  (3.1  days).  Less  fluctuation  about  the  analyzed 
vorticity  in  the  early  stages  of  development  is  observed  in  the  PS  (Figure  3.21)  than  in  the 
SCS  (Figure  3.20),  with  the  profile  of  analyzed  relative  vorticity  for  the  PS  region  subset 
more  closely  resembling  the  analyzed  vorticity  profile  for  all  non- developers  (Figure 
3.18).  The  average  analyzed  vorticity  for  this  subset  of  non- developers  was  3.2  x  10“^ 
s‘\  which  is  just  slightly  less  than  the  average  analyzed  vorticity  for  all  non- developers 
(3.3  X  10“^  s'^). 

As  in  the  SCS  subregion,  little  deviation  between  the  analyzed  and  forecast  values 
was  noted  in  the  short-range  (<  48  h)  forecasts  for  the  PS  subregion  in  the  hours  leading 
up  to  No.  This  lack  of  significant  error  implies  that  the  NOGAPS  model  accurately 
forecasts  the  early  stages  of  development  for  these  vortices.  In  the  hours  following  No 
the  analyzed  vorticity  decreases  while  the  forecast  vorticity  increases.  This  deviation  is 
associated  with  the  tendency  of  the  NOGAPS  model  to  continue  de’selopment  of  a 
circulation  that  is  actually  beginning  to  dissipate.  The  amount  of  deviation  between  the 
analyzed  and  forecast  relative  vorticity  increases  with  larger  forecast  range,  as  does  the 
forecast  time  of  maximum  analyzed  relative  vorticity.  However,  the  difference  between 
forecast  and  analyzed  vorticity  after  No  is  not  as  large  over  the  PS  as  it  is  over  the  SCS. 

Interestingly,  the  analyzed  850-mb  relative  vorticity  increases  again  48  h  after  No 
in  Figure  3.21.  At  this  secondary  maximum,  which  exists  only  within  the  PS  subregion, 
the  850-mb  relative  vorticity  is  roughly  4.0  x  10“^  s'\  which  is  greater  than  the  average 
analyzed  vorticity  at  No  for  both  the  subregion  (3.2  x  10'^  s'^)  and  all  the  non- developing 
vortices  (3.3  x  10“^  s'^).  This  secondary  maximum  could  be  associated  with  the 
reintensification  of  the  non- developing  vortex.  However,  since  analyzed  vorticity  again 
decreases  after  this  time,  the  vortex  does  not  undergo  further  development,  which  implies 
this  reintensification  may  be  associated  with  a  short-term  increase  in  convection. 
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Figure  3.21.  Average  analyzed  and  forecast  850- mb  relative  vorticity  (x  10'^  s  ')  for  non¬ 
developing  vortices  relative  to  ISf)  as  in  Figure  3.18,  except  within  the  PS  subregion. 
Heavy  black  line  represents  analyzed  average  850-mb  relative  vorticity  for  non¬ 
developing  vortices  (at  -i-OO),  and  red  line  represents  average  forecast  850-mb  relative 
vorticity  for  the  indicated  forecast  time.  Average  analyzed  850-mb  relative  vorticity  for 
this  subset  of  non- developing  vortices  at  No  (light  dashed  black  line)  was  3.2  x  10'^  s  '. 
?0  (light  solid  black  line)  is  5.0  x  10'^  s  '. 


The  average  duration  of  non- developing  vortices  in  the  EMT  subregion  was  2.7 
days,  which  is  only  slightly  less  that  the  duration  for  vortices  developing  in  the  PS 
subregion.  The  shorter  duration  of  these  non- developing  EMT  vortices,  in  direct  contrast 
to  the  duration  of  the  EMT  developing  vortices  (5.0  days),  may  be  associated  with  the 
shorter  life  cycle  of  vortices  not  associated  with  the  main  body  of  the  monsoon  trough. 

The  analyzed  850-mb  relative  vorticity  for  the  non- developing  vortices  in  the 
EMT  subregion  (Eigure  3.22)  closely  resembles  the  analyzed  vorticity  for  both  the  PS 
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Figure  3.22.  Average  analyzed  and  forecast  850- mb  relative  vorticity  (x  10'^  s'^)  for  non¬ 
developing  vortices  relative  to  No  as  in  Figure  3.18,  except  for  non- developing  vortices 
within  the  EMT  subregion.  Heavy  black  line  represents  analyzed  average  850- mb 
relative  vorticity  for  non- developing  vortices  (at  -tOO),  and  red  line  represents  average 
forecast  850- mb  relative  vorticity  for  the  indicated  forecast  time.  Average  analyzed  850- 
mb  relative  vorticity  for  this  subset  of  non- developing  vortices  at  No  (light  dashed  black 
line)  was  3.0  x  10'^  s ?o  (light  solid  black  line)  is  5.0  x  10'^  s 


subregion  (Figure  3.21)  and  the  entire  basin  (Figure  3.18).  However,  somewhat  less 
drastic  deviations  from  the  average  analyzed  vorticity  at  N)  are  noted.  As  in  the  SCS 
subregion  (Figure  3.20),  the  analyzed  vorticity  for  the  EMT  non- developing  vortices 
appears  to  fluctuate  through  several  minor  cycles  of  intensification  and  weakening,  which 
may  partially  be  a  result  of  the  sample  size. 

Although  the  vorticity  forecasts  consistently  exceed  the  analyzed  vorticity  values, 
the  24-h  forecasts  contain  the  general  increase  and  subsequent  decrease  in  vorticity 
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relative  to  No,  and  have  minimal  deviations  from  the  analyzed  mean  values.  The 
fluctuations  in  forecast  relative  vorticity  that  are  missing  in  the  48-,  72-  and  96-h 
forecasts  are  likely  associated  with  small  sample  size.  Although  the  previous  trend  is  for 
increasing  forecast  error  with  increased  forecast  range,  the  magnitudes  are  not  as 
excessive  in  the  EMT  as  in  the  overall  sample  (Figure  3.18). 

Also  of  interest  is  the  decrease  in  the  number  of  extreme  forecasts  (open  circles) 
in  the  EMT  subregion,  and  the  smaller  standard  deviation  of  the  analyzed  values  (light 
green  bars  at  each  time  interval  in  Figure  3.22.  In  addition  to  being  associated  with  the 
smaller  sample  size,  this  is  also  due  to  a  smaller  spread  in  forecast  vorticity  values. 
While  a  smaller  forecast  spread  does  not  directly  indicate  greater  forecast  accuracy,  it 
does  increase  the  probability  of  generating  a  more  accurate  forecast.  This  potentially 
increased  analysis  and  forecast  accuracy  in  the  EMT  may  occur  because  the  EMT 
vortices  can  be  more  easily  distinguished  from  the  background  environment  than  those 
that  occur  in  conjunction  with  the  monsoon  trough  in  the  PS  and  SCS  subregions. 

(2)  Deep  Eayer  Wind  Shear.  As  for  the  developing  vortices,  the 
terminology  for  over-  and  under- forecast  must  be  outlined  since  the  deep  layer  wind 
shear  variable  can  be  either  positive  or  negative  (indicating  direction).  Negative  shear 
indicates  increasing  easterlies  with  height,  and  positive  shear  indicates  the  opposite,  with 
increasing  westerlies  with  height.  The  term  over- forecast  is  applied  to  forecasts  of 
decreasing  westerly  wind  shear  if  they  become  less  positive,  or  to  decreasing  easterly 
shear  if  the  forecasts  become  less  negative. 

While  the  slope  of  the  analyzed  shear  with  time  for  the  developing  vortices  in 
Figure  3.14  decreases  steadily  throughout  the  life  cycle  of  the  tropical  circulation,  this 
trend  is  not  apparent  for  the  non- developing  vortices  that  met  the  minimum  duration 
criterion  (Figure  3.23).  The  analyzed  deep  wind  shear  value  at  No  for  the  NDl  data  set 
was  4.1  m  s'\  compared  to  the  -1.8  m  s'^  analyzed  at  Fq  for  developing  vortices.  Thus, 
the  developing  storms  experienced  easterly  shear  that  is  considered  to  be  favorable  for 
tropical  cyclone  formation.  However,  the  non- developing  vortices  experienced  an  almost 
equal  amount  of  westerly  shear,  which  is  considered  to  be  less  favorable.  This  easterly 
shear  remained  relatively  constant  in  time  relative  to  No  between  -60  h  and  -i-60  h  when  a 
reasonable  number  of  cases  are  available.  Thus,  this  variable  may  be  quite  useful  in 
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Figure  3.23.  Average  analyzed  and  foreeast  deep  layer  wind  shear  (200-850  mb)  for  non¬ 
developing  vortiees  relative  to  No,  as  in  Figure  3.18.  Heavy  blaek  line  represents 
analyzed  average  deep  wind  shear  for  non- developing  vortices  (at  -i-OO),  and  red  line 
represents  average  forecast  deep  wind  shear  for  the  indicated  forecast  time.  Average 
analyzed  deep  wind  shear  at  No  (light  dashed  black  line)  was  4.1  m  s'^  Vo  (light  solid 
black  line)  is  -1.8  m  s 


distinguishing  between  developing  and  non- developing  circulations  if  it  can  be  accurately 
forecast. 

Forecasts  of  deep  shear  for  the  NDl  data  set  are  almost  all  less  than  the  analyzed 
shear  values  at  each  of  the  four  forecast  intervals  in  Figure  3.23.  Forecast  deviations 
from  the  analyzed  values  increase  with  increasing  forecast  range.  Whereas  only  one  (-1-24 
h)  difference  between  mean  forecast  and  mean  analyzed  value  for  the  24-h  forecasts 
(upper- left  panel  of  Figure  3.23)  is  statistically  significant,  five,  six,  and  five  values  are 
significantly  different  for  the  48-h,  72-h,  and  96-h  forecasts,  respectively,  commensurate 
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with  the  increases  in  forecast  relative  vorticity  discussed  above.  The  average  deep  shear 
forecasts  are  also  preferentially  lower  than  the  analyzed  deep  shear,  and  do  not  transition 
from  forecasts  of  westerly  to  easterly  shear  until  approximately  24  h  after  Nq.  Forecasts 
of  easterly  shear  occur  mainly  in  the  extended  range  (>  72-h)  forecasts,  and  correspond  to 
vorticity  forecasts  that  exceed  the  ?o  threshold  (Figure  3.18).  Such  a  decreasing  forecast 
vertical  shear  toward  near- zero  values  is  considered  to  be  favorable  for  cyclonic 
development  as  the  warm  core  will  be  less  ventilated.  Thus,  this  factor  is  a  plausible 
physical  explanation  for  the  over- forecast  of  850-mb  relative  vorticity  (Figure  3.18)  for 
this  sample  of  non- developing  vortices.  What  aspect  of  the  NOGAPS  model  that  might 
lead  to  this  decreasing  westerly  shear  over  non- developers  is  not  clear.  One  hypothesis 
might  be  that  the  cumulus  momentum  transfer  is  too  larger  and  tends  to  minimize  vertical 
wind  shear. 

Comparison  of  Figures  3.5  and  3.14  for  the  developing  storms  indicates  that  the 
transition  from  westerly  to  easterly  analyzed  shear  occurs  when  analyzed  850-mb  relative 
vorticity  exceeds  about  4  x  10“^  s This  is  also  true  for  the  non- developing  vortices  in 
Figure  3.18  and  3.23.  When  the  forecast  850-mb  vorticity  in  Figure  3.18  exceeds  about  4 
X  10'^  s'^  (-1-36  h  for  the  48-h  forecasts),  the  forecast  deep  wind  shear  in  Figure  3.23 
transitions  from  westerly  to  easterly.  This  transition  from  westerly  to  the  preferential 
easterly  deep  wind  shear  does  not  occur  in  the  NDl  analyzed  curve  since  the  analyzed 
850-mb  relative  vorticity  for  this  data  set  does  not  exceed  4  x  10'^  s It  is  thus  possible 
that  this  westerly  shear  prevented  further  development  of  the  vortices  in  this  data  set,  but 
the  NOGAPS  forecast  of  a  decreasing  vertical  wind  shear  in  time  allowed  same  non¬ 
developing  cases  to  intensify  to  4  x  10'^  s 

(3)  Sea-level  Pressure.  The  analyzed  mean  sea-level  pressures  for  the 
NDl  data  set  (Figure  3.24)  showed  very  little  change  relative  to  18b  and  throughout  the 
life  cycle  of  these  non- developing  vortices.  At  all  times,  the  analyzed  SLP  for  the  NDl 
data  set  was  greater  than  the  average  analyzed  SLP  for  developing  vortices  (Figure  3.15). 
The  average  analyzed  SLP  at  No  was  1008.0  mb,  which  is  1.4  mb  higher  than  the 
analyzed  SLP  for  developing  vortices  (1006.6  mb).  This  difference,  while  small,  is 
important  in  that  it  implies  a  less  deep  trough  within  which  the  vortex  is  trying  to 
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Figure  3.24.  Average  analyzed  and  forecast  sea  level  pressure  (SLP)  (mb)  for  non¬ 
developing  vortices  relative  to  No,  as  in  Figure  3.18.  Heavy  black  line  represents 
analyzed  average  SLP  for  non- developing  vortices  (at  -i-OO),  and  red  line  represents 
average  forecast  SLP  for  the  indicated  forecast  time.  Average  analyzed  SLP  at  No  (light 
dashed  black  line)  is  1008.0  mb.  Po  (light  solid  black  line)  for  developing  circulations 
is  1006.6  mb. 


develop.  Slightly  higher  SLP  implies  a  less-developed  warm-core  aloft.  In  this  sense, 
the  SLP  may  be  more  of  a  proxy  as  to  the  degree  of  development  rather  than  a  predictor. 

The  slight  increase  in  the  slope  of  the  analyzed  SLP  line  in  the  24  h  following  No 
is  consistent  with  the  weakening  of  the  analyzed  850- mb  relative  vorticity  (Figure  3.18) 
during  the  same  time  frame.  While  the  slope  of  the  analyzed  SLP  for  developing 
circulations  (Figure  3.15)  decreased  throughout  the  120  h  prior  to  and  following  Fb,  no 
commensurate  decrease  is  found  for  the  analyzed  SLP  of  the  non- developing  vortices. 
While  part  of  the  decrease  after  lb  in  analyzed  SLP  for  the  developing  vortices  was 
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associated  with  the  hsertion  of  the  synthetic  tropical  cyclone  observations,  no  similar 
insertion  would  occur  for  the  non- developing  vortices,  which  leads  to  the  NDl  analyzed 
SLP  profile  more  closely  resembling  the  early  stages  of  the  developing  vortices  in  Figure 
3.15,  i.e.,  prior  to  intensification  of  the  developing  circulation. 

Forecasts  of  SLP  for  the  non- developing  vortices  are  typically  for  a  lower  SLP 
than  the  analyzed  SLP  (Figure  3.24).  Deviations  of  the  forecast  SLPs  from  the  analyzed 
SLPs  increase  with  increasing  forecast  range  from  24  h  through  96  h.  In  particular,  the 
forecast  SLPs  for  the  NDl  data  set  at  72  h  (lower-left  panel  in  Figure  3.24)  and  at  96  h 
(lower-right  panel)  are  consistently  at  or  below  the  SLP  threshold  (1006.6  mb) 
determined  for  the  developing  vortices.  These  decreases  in  forecast  SLP  can  be 
associated  with  the  over- forecast  850- mb  relative  vorticity  at  these  longer  forecast  ranges 
in  Figure  3.18. 

Forecast  error  also  increases  with  forecast  range  in  Figure  3.24,  which  again 
suggests  a  systematic  bias  in  the  heating  over  these  non- developing  cases.  Notice  that  the 
forecast  extrema  in  Figure  3.24  are  considerably  lower  than  the  analyzed  SLP  for  non¬ 
developing  vortices,  and  these  outliers  have  SLPs  well  below  those  of  the  mean  analyzed 
SLP  for  the  developing  cases  than  the  analyzed  SLP  for  the  non- developing  cases.  These 
low  extrema  are  expected  given  the  NOGAPS  model  tendency  to  over-  forecast  vorticity 
and  generally  to  over- intensify  tropical  circulations.  The  high  SLP  extrema  were  not 
expected.  Although  infrequent,  these  high  SLP  values  represent  NOGAPS  forecasts  that 
substantially  weaken  the  circulation,  or  perhaps  can  be  attributed  to  erroneous  tracking  of 
systems  closer  to  the  subtropical  high. 

(4)  925-mb  Wind  Speed.  The  analyzed  925- mb  wind  speed  (Figure 
3.25)  has  relatively  small  (~1-L5  m  s'^)  changes  to  No,  although  the  slope  does  increase 
steadily  from  60  h  prior  to  No,  decrease  in  the  36  h  immediately  following  No,  and  then 
increases  again  after  36  h.  The  steady  increase  in  wind  speed  prior  to  No  is  attributed  to 
the  early  developmental  stages  of  the  potential  tropical  cyclone.  While  the  value  of  the 
wind  speed  at  No  (7.0  m  s'^)  is  a  local  maximum,  it  is  not  the  maximum  analyzed  value 
for  the  entire  forecast  period.  However,  the  maximum  value  (about  8  m  s‘^)  occurs  120  h 
after  No,  but  it  is  suspect  since  it  is  based  on  fewer  than  five  cases.  The  decrease  in  wind 
speed  in  the  hours  immediately  following  No  is  consistent  with  the  decrease  in  850- mb 
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Figure  3.25.  Average  analyzed  and  forecast  925-mb  wind  speed  (m  s'^)  for  nor+ 
developing  vortices  relative  to  No,  as  in  Figure  3.18.  Heavy  black  line  represents 
analyzed  average  925-mb  wind  speed  for  nor+ developing  vortices  (at  +00),  and  red  line 
represents  average  forecast  925-mb  wind  speed  for  the  indicated  forecast  time.  Average 
analyzed  925-mb  wind  speed  at  No  (light  dashed  black  line)  is  7.0  m  s Vo  (light  solid 
black  line)  is  10.8  m  s 


relative  vorticity  (Figure  3.18)  that  is  attributed  to  a  broadening  and  weakening  of  the 
tropical  circulation.  The  increase  in  wind  speed  that  begins  36  h  after  No  follows  the 
slight  increase  in  analyzed  relative  vorticity  that  occurs  at  the  same  time. 

The  difference  of  about  4  ms'  between  these  925-mb  wind  speeds  for  non¬ 
developing  cases  and  developing  cases  potentially  offers  another  distinguishing 
characteristic  for  formation.  Since  this  wind  speed  has  been  averaged  over  the  ellipse 
centered  on  the  position,  it  is  not  surprising  that  the  development  sample  has  larger 
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values  at  later  times  than  the  non- developing  cases.  However,  the  standard  deviations  of 
the  two  samples  are  relatively  small  compared  to  the  differences  in  the  means. 

Forecasts  of  925- mb  wind  speed  for  the  NDl  data  set  are  consistently  higher  than 
the  analyzed  wind  speed  for  the  same  vortices.  However,  forecasts  of  wind  speed  do  not 
exceed  the  wind  speed  threshold  determined  for  developing  vortices  (10.8  m  s‘^).  As  was 
the  case  for  the  forecasts  of  non- developing  vortices  for  other  variables,  the  forecast 
spread  and  errors  increase  with  increasing  forecast  range,  and  both  the  magnitude  and 
occurrence  time  of  the  maximum  forecast  wind  speed  also  increase.  However,  the 
forecast  means  for  these  non- developers  do  exceed  the  mean  value  for  the  developing 
cases  (10.8  m  s'^)  even  for  the  96- h  forecasts. 

The  outlier  values  of  wind  speed  are  all  greater  than  the  analyzed  wind  speed  at 
No,  and  also  at  Fq  for  the  developing  storms,  which  is  consistent  with  the  over- forecasts 
of  850- mb  relative  vorticity.  These  extrema  indicate  a  tendency  of  the  NOGAPS  model 
to  over- forecast  925- mb  wind  speed,  as  expected  in  a  hydrostatically  and  geostrophically 
balanced  model  that  forecasts  too  large  850- mb  relative  vorticity  and  too  low  SLP. 

(5)  Vapor  Pressure.  Average  500-700  mb  vapor  pressure,  which 
indicates  the  amount  of  mid-  level  moisture  available  to  a  developing  tropical  circulation, 
has  more  variation  relative  to  the  threshold  value  for  non- developing  vortices  (Figure 
3.26)  than  for  the  developing  vortices  (Figure  3.17).  The  average  analyzed  vapor 
pressure  at  No  (3.4  Pa)  and  Fo  (3.6  Pa)  differ  by  only  0.3  Pa,  which  makes  the  vapor 
pressure  a  poor  indicator  of  potential  development. 

The  analyzed  vapor  pressure  for  non- developing  vortices  increases  steadily  from 
96  hours  prior  to  No  to  a  maximum  value  at  No-12  h,  and  then  steadily  decreases  until  36 
h  after  No.  This  trend  indicates  intensification  and  then  weakening  of  the  non- developing 
warm-core  vortex.  The  increase  in  vapor  pressure  36  h  after  No  is  attributed  to  a  weak 
and  unsuccessful  reintensification  of  the  warm  core.  That  the  analyzed  vapor  pressure 
never  exceeds  the  threshold  value  for  the  developing  circulations  indicates  that  in 
addition  to  a  drier  mid- troposphere,  the  potential  warm  core  vortex  is  somewhat  weaker. 

The  forecast  curves  for  vapor  pressure  lie  almost  entirely  below  the  analyzed 
vapor  pressure  curve,  and  for  the  most  part,  are  also  below  the  average  analyzed  vapor 
pressure  at  No.  This  small  but  consistent  under- forecasting  of  vapor  pressure  by  the 
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Figure  3.26.  Avera^  analyzed  and  forecast  vapor  pressure  (500-700  mb  average,  Pa)  for 
non- developing  vortices  relative  to  ISb,  as  in  Figure  3.18.  Heavy  black  line  represents 
analyzed  average  vapor  pressure  for  non- developing  vortices  (at  +00),  and  red  line 
represents  average  forecast  vapor  pressure  for  the  indicated  forecast  time.  Average 
analyzed  vapor  pressure  at  ISb  (light  dashed  black  line)  was  3.4  Pa.  Fb  vapor  pressure 
threshold  (light  solid  black  line)  is  3.6  Pa. 


NOGAPS  model  is  likely  a  result  of  incorrect  representation  of  the  mid- level  mesoscale 
warm  core,  and  again  shifts  formation  from  a  mesoscale  event  in  the  NOGAPS  model 
forecasts  to  a  large-scale  event. 

As  with  other  forecast  variables,  forecast  spread  and  error  increase  with 
increasing  forecast  range.  The  outlier  values,  while  not  all  high  or  low  compared  to  the 
average  analyzed  vapor  pressure  at  any  forecast  interval,  do  have  some  exceedingly  dry 
values.  Such  dry  mid -tropospheric  values  would  not  be  consistent  with  formation,  as 
noted  previously. 
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b.  Model  Over-forecasts  (NDOF) 

Model  over- forecast  vortices  were  identified  as  those  for  which  the  NOGAPS 
model  forecast  continued  development  of  a  circulation  that  did  not  meet  the  minimum 
duration  criterion  (existence  in  three  consecutive  model  analyses).  When  determining 
whether  or  not  an  analyzed  vortex  was  over- forecast,  some  simple  rules  were  applied.  If 
a  vortex  existed  only  once  in  an  analysis  panel  (-i-OO),  then  there  had  to  be  at  least  three 
corresponding  forecast  groups  for  it  to  be  classified  as  an  over- forecast.  If  the  vortex 
existed  for  two  consecutive  analyses,  then  it  had  to  have  at  least  four  consecutive  forecast 
groups  in  common  with  the  two  analyses  for  it  to  be  classified  as  an  over- forecast. 
Vortices  that  existed  in  three  or  more  analyses  were  not  included  in  this  category,  but 
were  instead  included  in  the  NDl  data  set.  Although  the  forecast  circulation  then  existed 
at  times  when  an  analogous  circulation  (of  at  least  ?  =  1.5  x  10'^  s'^)  could  not  be 
matched,  this  is  not  a  false  alarm  -  it  is  over-  forecast  in  length  of  time  and  not  necessarily 
in  magnitude  of  some  variable  such  as  SLP.  Of  the  127  vortices  rejected  for  not  meeting 
the  minimum  duration  aiterion,  60  were  classified  as  over- forecasts  (column  three  of 
Table  3.1).  Thus,  this  situation  is  a  relatively  common  occurrence  that  needs  to  be 
examined. 

An  example  of  the  analyzed  and  forecast  tracks  for  an  over-  forecast  circulation  is 
shown  in  Figure  3.27.  This  circulation  was  initially  rejected  from  the  dataset  since  it  did 
not  exist  in  a  NOGAPS  analysis  for  at  least  three  consecutive  periods.  However,  it  did 
meet  the  above  criteria  for  an  over-forecast  circulation,  since  it  existed  for  two  analyses 
and  at  least  four  consecutive  model  forecasts.  The  PS  subregion  contained  the  most  over¬ 
forecast  circulations  (28)  compared  to  the  SCS  (18)  and  EMT  (16)  subregions. 

(1)  850-mb  Relative  Vorticity.  Since  the  minimum  duration  criterion 
was  not  met  for  the  over- forecast  non- developing  vortices,  a  new  reference  time  had  to  be 
determined.  All  the  vortices  that  met  the  over- forecast  non-  developing  vortex  criteria 
existed  in  at  least  one  NOGAPS  model  analysis.  The  new  reference  time  (designated 
No*)  from  which  forecast  time  intervals  were  calculated  was  defined  to  be  the  last  time 
the  vortex  existed  in  a  NOGAPS  model  analysis,  rather  than  the  time  at  which  the 
maximum  analyzed  850-mb  vorticity  occurred.  Given  the  criteria  for  over- forecast  non¬ 
developing  vortices,  the  analyzed  850-mb  relative  vorticity  values  are  available  for  at 


68 


Figure  3.27.  Example  of  forecast  tracks  for  an  over- forecast  non- developing  vortex. 
Circles  indicate  successive  12  h  forecasts  of  the  circulation  center,  with  the  size  of  the 
circle  indicating  the  intensity.  Black  circles  indicate  analyzed  values,  and  colored  circles 
indicate  separate  forecast  times. 

most  only  two  intervals  (Figure  3.28).  On  average,  the  analyzed  vorticity  at  No*  (2.4  x 
10'^  s'^)  is  slightly  less  than  the  analyzed  vorticity  12  h  prior  to  No*,  which  implies  a 
slight  weakening  of  the  analyzed  vortex. 

The  forecasts  of  850- mb  relative  vorticity  are  the  defining  feature  of  this  subset  of 
non- developing  vortices  because  each  of  these  over- forecast  vortices  was  defined  to  exist 
in  the  forecasts  for  at  least  twice  as  long  as  it  existed  in  the  NOGAPS  analysis.  The 
forecast  vorticities  (Figure  3.28)  are  nearly  all  greater  than  the  analyzed  vorticities  with 
very  few  exceptions.  The  longer-range  forecasts  (>  72-h)  tend  to  not  only  have  longer 
durations,  but  also  to  forecast  larger  vorticity  magnitudes  than  the  shorter-range  forecasts. 
Whereas  the  analyzed  vorticity  decreases  from  No*- 12  h  to  No*,  the  vorticity  is  forecast 
to  increase  during  and  following  that  time.  This  tendency  for  the  NOGAPS  model  to 
hold  on  to  dissipating  circulations  is  consistent  with  the  general  tendency  to  over- develop 
the  tropical  circulations  in  this  sample.  Although  the  mean  forecast  vorticities  do  not 
exceed  the  mean  analyzed  vorticity  for  developing  vortices  (?o  =  5.0  x  10'^  s'^),  some  of 
the  outliers  do  cross  this  threshold  and  thus  could  be  mistaken  for  developing  vortices. 
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Figure  3.28.  Average  analyzed  and  forecast  850- mb  relative  vorticity  (x  10'^  s'^)  for 
over- forecast  non- developing  vortices  as  in  Figure  3.18.  Heavy  black  line  represents 
analyzed  average  850- mb  relative  vorticity  for  over- forecast  non- developing  vortices  (at 
-1-00),  and  red  line  represents  average  forecast  850- mb  relative  vorticity  for  the  indicated 
forecast  time.  Average  analyzed  850- mb  relative  vorticity  for  this  subset  of  non¬ 
developing  vortices  at  No  (light  dashed  black  line)  was  2.4  x  10'^  s ?o  (light  solid  black 
line)  is  5.0  X  10"^  s  '. 


(2)  Deep  Layer  Wind  Shear.  As  was  the  case  for  the  overall  non¬ 
developer  sample,  the  analyzed  deep  layer  wind  shear  is  westerly  for  these  cases.  The 
magnitude  of  the  analyzed  deep  wind  shear  for  over- forecast  non- developing  vortices 
(Figure  3.29)  does  not  change  appreciably  from  No*- 12  h  to  No*.  This  lack  of  slope 
change  is  due  to  the  short  time  elapsed  between  the  two  analysis  intervals,  which  may  be 
insufficient  to  reflect  a  change  in  wind  shear  as  a  lesponse  to  the  decrease  in  850- mb 
relative  vorticity  at  No*  in  Figure  3.28. 


70 


-96  -72  -48  -24  00  +24  +48  +72  +96 


Forecast  Verification  Time 


Nq  2002  WPAC  NOGAPS  NDOF  at  +48 


20 


10 


+r 

«  0 
-a 


-10 


1  1  1  1  1  4  7  9  19231612  9  6  0  0  0  0  0 

oO 

o® 

iiii 

CD 

O 

1 

■ 

-96  -72  -48  -24  00  +24  +48  +72  +96 

Forecast  Verification  Time 


Ng  2002  WPAC  NOGAPS  NDOF  at  +72 


1  1  1  0  0  4  3  112220281711  6  5  4  0  0  0 

20- 


-20- 


-96  -72  -48  -24  00  +24  +48  +72  +96 

Forecast  Verification  Time 


Ng  2002  WPAC  NOGAPS  NDOF  at  +96 


<0 

"O 


-96  -72  -48  -24  OQ  ■f24  +48  +72  +96 

Forecast  Verification  Time 


Figure  3.29.  Average  analyzed  and  forecast  deep  wind  shear  (200-850  mb)  for  non¬ 
developing  over-forecast  vortices  relative  to  No.  Heavy  black  line  represents  analyzed 
average  deep  wind  shear  for  non- developing  over- forecast  vortices  (at  +00),  and  red  line 
represents  average  forecast  deep  wind  shear  for  the  indicated  forecast  time.  Average 
analyzed  deep  wind  shear  at  No  (light  dashed  black  line)  was  3.2  m  s'^  Vo  (light  solid 
black  line)  is  -1.8  m  s 


The  forecast  deep  layer  wind  shear  at  No*  is  very  close  to  the  analyzed  value  (3.2 
m  s'^).  There  is  no  strong  trend  in  the  NOGAPS  model  for  the  forecasts  to  be 
substantially  over-  or  under-forecast  relative  to  N)*-  The  forecasts  generally  indicate 
positive  (westerly)  shear,  which  was  shown  above  to  be  common  among  the  non¬ 
developing  vortices.  The  trend  identified  above  for  a  transition  from  westerly  to  easterly 
shear  when  the  850- mb  relative  vorticity  is  forecast  to  be  greater  than  about  4.0  x  10'^  s  ' 
is  not  as  apparent  in  Figure  3.29  as  it  is  in  Figures  3.23.  Although  more  of  the  outlier 
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vertical  shears  are  for  large  westerly  shears,  especially  for  the  72-h  forecasts,  some 
outliers  with  large  easterly  shear  are  forecast. 
c.  False  Alarms 

False  alarms  were  separated  into  three  categories:  the  non- developing  vortices 
that  met  the  minimum  duration  criterion,  and  were  forecast  to  exceed  the  vorticity 
threshold  determined  for  the  developing  circulations;  the  NOGAPS  model  false  alarms, 
in  which  the  model  predicted  development  of  a  vortex  that  was  never  analyzed;  and 
TCFA  false  alarms,  in  which  the  JTWC  issued  a  TCFA  for  a  circulation  that,  while  it 
could  be  usually  be  matched  with  a  circulation  in  the  NOGAPS  analysis,  was  never 
warned  on 

(1)  Non- Developing  Vortices  -  Group  2  (NDG2).  Group  2  represents 
the  subset  of  non- developers  that  met  the  minimum  duration  criterion  (NDl),  but  were 
relatively  extreme  in  the  sense  that  they  were  forecast  to  exceed  the  850-mb  relative 
vorticity  value  for  developing  vortices  (5.0  x  10'^  s'^).  A  total  of  50  vortices  met  the 
criteria  for  inclusion  in  this  data  set  (column  four  of  Table  3.1).  Whereas  the  cases 
contained  in  this  data  set  are  a  subset  of  the  NDl  data  set,  the  analysis  here  is  to 
determine  if  any  single  factor  prevented  these  favorable -forecast  vortices  from 
developing  into  numbered  tropical  circulations.  An  example  of  a  set  of  Group  2  tracks  is 
shown  in  Figure  3.30.  Intensity  is  indicated  by  the  size  of  the  circle  at  each  position.  The 
Group  2  circulations  formed  preferentially  in  the  PS  subregion  (21),  although  the  number 
of  Group  2  circulations  in  the  SCS  was  only  slightly  less  (19). 

(a)  850-mb  Relative  Vorticity.  More  variability  exists  in  the 

slope  of  the  analyzed  850-mb  relative  vorticity  for  Group  2  (Figure  3.31)  than  in  the 
overall  NDl  data  set.  The  slope  of  the  line  prior  to  the  maximum  analyzed  vorticity  at  No 
more  closely  resembles  the  slope  of  the  analyzed  vorticity  for  developing  vortices  in 
Figure  3.6  prior  to  F*o.  However,  the  slope  for  the  Group  2  vortices  is  slightly  more 
shallow  than  the  analyzed  vorticity  slope  for  the  overall  NDl  data  set  (Figure  3.18).  The 
sudden  decrease  in  analyzed  vorticity  at  No  is  attributed  to  a  weakening  of  the  vortex  in 
the  model  associated  with  dissipation. 
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Figure  3.30.  Example  of  forecast  tracks  for  a  Group  2  non- developing  vortex.  Circles 
indicate  successive  12  h  forecasts  of  the  circulation  center,  with  the  size  of  the  circle 
indicating  the  intensity.  Black  circles  indicate  analyzed  values,  and  colored  circles 
indicate  separate  forecast  times. 

The  average  analyzed  vorticity  at  No  for  this  data  set  was  3.7  x  10'^  s  ',  which  is 
between  the  threshold  for  developing  vortices  (5.0  x  10'^  s  ')  and  non- developing  vortices 
(3.3  X  10'^  s  '),  but  slightly  below  the  4.0  x  10'^  s  '  threshold  identified  for  transitioning  of 
deep  wind  shear  from  westerly  to  easterly  in  the  developing  vortices. 

Forecasts  of  850- mb  relative  vorticity  for  these  Group  2  cases  (Figure  3.31)  vary  only 
slightly  from  the  analyzed  vorticity  values  in  the  48  h  leading  up  to  N),  especially  for 
shorter-range  forecasts  (24-  and  48-h).  Following  Nq,  the  NOGAPS  model  continued  to 
forecast  development  of  the  vortex  for  a  minimum  of  24  h  before  decreasing  the  forecast 
vorticity.  The  most  dramatic  deviations  between  analyzed  and  forecast  vorticities  occur 
within  36  h  of  Nq.  Although  there  are  fewer  forecast  time  intervals  in  Figure  3.30  in 
which  the  deviation  between  analyzed  and  forecast  vorticities  are  statistically  significant 
than  in  the  forecast  and  analyzed  vorticities  of  the  overall  NDl  data  set  (Figure  3.18),  the 
errors  appear  more  noticeable,  possibly  due  to  the  small  sample  size. 

(b)  Deep  Layer  Wind  Shear.  Analysis  of  the  deep  wind  shear 
curve  for  the  Group  2  subset  (Figure  3.32)  indicates  that  the  analyzed  deep  wind  shear 
never  transitioned  from  westerly  to  easterly,  as  is  expected  with  analyzed  vorticity  that 
did  not  exceed  4.0  x  10'^  s  '  (refer  to  section  C.2.a.2  of  this  chapter).  Analysis  data  were 
not  available  at  72  h  and  120  h  after  Nq.  The  average  analyzed  deep  shear  at  No  was  2.8 
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Figure  3.31.  Average  analyzed  and  forecast  850- mb  relative  vorticity  as  in  Figure  3.18, 
except  for  the  Group  2  subset  of  NDl  vortices.  Heavy  black  line  represents  average 
analyzed  850-mb  relative  vorticity  for  the  Group  2  subset  of  vortices,  and  heavy  red  line 
represents  average  forecast  850-mb  relative  vorticity  for  the  time  indicated.  The  light 
solid  line  represents  average  analyzed  850-mb  relative  vorticity  at  FO  for  developing 
vortices  (5.0  x  10-5  s-1),  and  the  light  dashed  line  represents  average  analyzed  850-mb 
relative  vorticity  at  NO  (3.7  x  10-5  s-1). 


ms  *,  which  indicates  a  westerly  wind  shear,  which  has  been  shown  to  be  unfavorable  for 
tropical  cyclone  formation.  Within  24  h  of  No,  the  analyzed  wind  shear  not  only  remains 
positive  but  also  increases  in  magnitude,  which  is  not  favorable  for  development  of  the 
potential  tropical  circulation.  The  early  decrease  in  analyzed  shear  with  time  indicates  a 
more  favorable  condition  for  the  development  and  intensification  of  the  tropical 
circulation.  However,  the  lack  of  a  pronounced  decrease  in  the  slope  of  the  analyzed 
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Figure  3.32.  Average  analyzed  and  forecast  deep  layer  wind  shear  (200-850  mb)  (m  s'^) 
for  non- developing  vortices  that  were  forecast  to  exceed  Fb,  as  in  Figure  3.18.  Heavy 
black  line  represents  analyzed  average  deep  layer  wind  shear  for  this  subset  of  non¬ 
developing  vortices  (at  +00),  and  red  line  represents  average  forecast  deep  layer  wind 
shear  for  the  indicated  forecast  time.  Average  analyzed  deep  layer  wind  shear  850  mb 
relative  vorticity  for  this  subset  of  non- developing  vortices  at  No  (light  dashed  black  line) 
was  2.8  m  s'^  Vo  (light  solid  black  line)  is  -1.8  m  s 


deep  shear  as  in  Figure  3.15  would  not  lead  one  to  expect  the  eventual  development  of  a 
JTWC-numbered  tropical  circulation. 

In  general,  the  NOGAPS  deep  shear  forecasts  (Figure  3.32)  become  more 
negative  with  increasing  time  and  forecast  range.  At  each  forecast  time,  the  NOGAPS 
forecasts  of  deep  shear  become  distinctly  negative  within  24  h  of  No,  which  indicates  a 
transition  of  the  shear  from  positive  to  negative,  and  thus  more  favorable  conditions  for 
development.  As  noted  earlier,  the  transition  from  positive  to  negative  shear  tends  to 
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occur  when  the  850- mb  relative  vorticity  is  forecast  to  exceed  4.0  x  10“^  s Since  the 
analyzed  vorticity  24  h  after  No  is  distinctly  below  this  threshold,  this  may  be  related  to 
the  fact  that  the  deep  wind  shear  remained  westerly. 

Interestingly,  the  outlier  values  in  Figure  3.32  are  almost  all  positive  while  the 
forecast  range  at  each  time  interval  is  skewed  toward  negative  shear.  The  range  of  more 
negative  forecast  values  indicates  that  the  NOGAPS  model  over-  forecast  the  shear  to  be 
smaller  than  analyzed.  However,  the  positive  deep  layer  shear  extrema  that  have  no 
corresponding  erroneously  low  vorticity  forecasts  must  occur  as  a  result  of  some  other 
model  forecast  error. 

(c)  Other  Variables.  Trends  for  the  other  selected  forecast 
variables  for  the  Group  2  subset  were  very  similar  to  those  for  the  overall  NDl  data  set. 
The  most  noticeable  difference  was  the  same  one  evident  in  the  comparison  of  850-mb 
relative  vorticity  for  Group  2  and  NDl  -  the  Group  2  analyzed  and  forecast  curves  were 
more  peaked  than  those  for  the  overall  NDl  set,  and  more  closely  resembled  the  curves 
for  developing  vortices,  although  both  the  analyzed  and  forecast  values  for  the  non¬ 
developers  were  smaller  in  magnitude  than  their  counterparts  for  developing  circulations. 

The  relatively  small  sample  size  may  have  contributed  to  the  relatively  large 
variation  evident  in  the  analyzed  curves  for  Group  2.  This  variation,  while  not  readily 
apparent  in  Figure  3.21,  was  distinctly  more  noticeable  in  the  analyses  and  forecasts  for 
some  of  the  other  variables,  especially  size,  and  vapor  pressure  (not  shown). 

(2)  Model  False  Alarms .  Model  false  alarms  were  identified  from  the 
forecast  vortices  stored  in  the  Forecast  Database  (Figure  2.1a)  that  never  corresponded  to 
a  verifying  analysis.  To  qualify  as  a  model  false  alarm,  four  criteria  had  to  be  met:  the 
vortex  had  to  exist  for  at  least  24  hours  (three  consecutive  model  runs)  in  the  forecast 
sequence;  a  potential  match  had  to  exist,  meaning  that  the  false  alarm  forecast  circulation 
had  to  be  matched  to  at  least  one  other  forecast;  the  match  had  to  be  within  a  4°  latitude 
radius  circle  of  the  vortex  it  was  being  matched  to;  and  the  match  had  to  occur  within  + 
12  hours  of  when  it  was  forecast  to  occur.  Given  these  criteria,  15  vortices  were 
identified  as  false  alarms  (column  five  of  Table  3.1).  The  850-mb  relative  vorticity  for 
the  model  false  alarms  could  not  be  compared  to  the  developing  circulation  vorticity 
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threshold  since  no  first  JTWC  warning  time  or  time  of  maximum  analyzed  vorticity 
existed. 

The  NOGAPS  predictions  had  more  false  alarm  forecasts  in  the  Philippine  Sea 
(10)  than  in  the  South  China  Sea  (4)  and  Eastern  Monsoon  Trough  (1)  combined  (Figure 
3.33).  This  regional  distribution  is  consistent  with  the  results  of  Cheung  and  Elsberry 
(2003).  That  more  false  alarms  occurred  subregions  associated  with  the  monsoon  trough 
indicates  that  the  NOGAPS  model  may  have  over- forecast  the  vorticity  associated  with 
smaller  circulations  within  the  monsoon  trough.  As  discussed  previously,  over- forecast 


Figure  3.33.  Eocations  of  NOGAPS  model  false  alarms  in  the  western  North  Pacific 
during  1  May  -  31  October  2002.  Circles  indicate  successive  12-h  analyses  of  each 
identified  model  false  alarm.  The  size  of  each  circle  is  related  to  the  analyzed  intensity  of 
the  circulation  with  a  larger  circle  indicating  a  more  intense  circulation. 

vorticity  leads  to  forecasts  of  excessively  easterly  deep  layer  wind  shear  which  is 
favorable  for  tropical  cyclone  formation.  That  more  false  alarms  do  not  occur  is  possibly 
due  to  the  excessive  drying  of  the  mid- troposphere,  as  noted  above. 

(3)  TCFA  False  Alarms.  During  the  study  period,  there  were  six 
TCFA  false  alarms.  Since  so  few  TCFA  false  alarms  existed  during  2002,  the  minimum 
sample  size  was  relaxed  for  this  dataset.  Four  of  these  six  false  alarms  could  be 
associated  with  non- developing  vortices  tracked  in  the  NOGAPS  analyses  (Figure  3.34). 
Additionally,  four  of  the  six  TCFA  false  alarms  occurred  in  the  Philippine  Sea  (PS),  and 
three  of  those  occurred  between  136  °E  and  142°E.  The  first  two  false  alarms 
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Figure  3.34.  Locations  of  JTWC  TCFA  false  alarms  in  the  western  North  Pacific  during 
1  May  -  31  October  2002.  Circles  indicate  successive  12-h  analyses  of  each  identified 
model  false  alarm. 

occurred  within  6  hours  of  one  another,  and  are  likely  for  the  same  vortex.  One  false 
alarm  in  the  PS  subregion  was  for  a  vortex  near  30°N,  and  thus  might  have  been  a 
subtropical  development.  The  other  two  TCFA  false  alarms  occurred  in  the  South  China 
Sea  within  three  days  of  each  other.  It  is  possible  that  all  six  of  these  vortices  represented 
weak  circulations  associated  with  the  monsoon  trough  or  monsoon  depressions. 

Although  never  warned  on,  the  four  of  the  vortices  corresponding  to  TCFA  false 
alarms  did  exist  in  at  least  one  analysis  panel.  Therefore,  No  could  be  determined.  The 
mean  analyzed  vorticity  (at  No)  for  this  subset  of  non- developing  vortices  (Figure  3.35) 
was  7.1  X  10'^  s'*,  which  is  about  2  x  10'^  s'*  higher  than  the  vorticity  threshold 
determined  for  the  developing  vortices  (see  Figure  3.5).  The  larger  average  analyzed 
vorticity  for  this  subset  of  non- developers  is  unexpected  for  a  series  of  vortices  assumed 
to  be  associated  with  weak,  disorganized  monsoon  depressions. 

The  slope  of  the  analyzed  vorticity  for  the  four  vortices  included  in  Figure  3.35 
dramatically  increases  from  36  h  prior  to  N)  and  then  decreases  sharply  in  the  12  h 
following  No.  The  slope  of  the  analyzed  vorticity  for  this  subset  of  non- developers  is 
steeper  than  the  slope  of  the  analyzed  vorticity  for  the  developing  vortices  in  Figure  3.5, 
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Figure  3.35.  Average  analyzed  and  forecast  850- mb  relative  vorticity  (x  10'^  s'^)  as  in 
Figure  3.5,  except  for  the  four  TCFA  false  alarm  vortices  with  corresponding  NOGAPS 
analyses.  Heavy  black  line  represents  analyzed  average  850- mb  relative  vorticity  (at 
-tOO)  for  this  subset  of  non- developers,  and  red  line  represents  average  forecast  850- mb 
relative  vorticity  for  this  subset  of  non- developers  for  the  indicated  forecast  time. 
Average  analyzed  vorticity  relative  at  the  time  of  maximum  analyzed  vorticity  No  (light 
dotted  black  line)  for  these  TCFA  false  alarm  vortices  was  7.1  x  10'^  s'\  which  is  larger 
than  the  5.0  x  10'^  s'*  for  developing  vortices  at  Fq  as  in  Figure  3.5. 


but  more  closely  resembles  the  slope  of  the  analyzed  vorticity  prior  to  Fq  for  the  subset  of 
developers  that  verified  prior  to  the  TCFA- specified  formation  window  (Figure  3.11). 
This  increase  and  sudden  decrease  is  attributed  to  temporary  spin- ups  of  vorticity  that 
may  be  associated  with  overly  active  convection  in  the  model.  However,  since  four  of 
these  alerts  corresponded  to  circulations  that  existed  in  the  NOGAPS  analysis,  which 
likely  were  also  tracked  concurrently  using  satellite  imagery,  the  NOGAPS  analysis 
probably  also  reflects  a  relatively  strong  circulation  in  connection  with  the  convection. 
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The  variability  of  forecast  range  at  each  time  interval  for  longer-range  forecasts  is 
evident  in  the  96-h  forecasts  of  this  data  set  (lower-right  panel  of  Figure  3.35).  While  the 
48-  through  72-h  forecasts  predict  vorticity  substantially  greater  than  the  analyzed  value, 
they  are  suspect  since  only  two  cases  contributed  to  the  mean  calculated  at  each  time 
interval.  However,  that  the  vorticity  forecast  at  these  three  verification  times  is  so  much 
greater  (and  a  statistically  significant  deviation  at  48  h  and  72  h)  than  the  analyzed  value 
indicates  that  at  longer  time  ranges  the  NOGAPS  model  has  again  over- forecast  the 
development,  or  more  accurately  the  apparent  dissipation  of  these  circulations. 

D.  CASE  STUDY  COMPARISON 

To  illustrate  the  effects  of  the  forecast  thresholds  determined  previously,  two  case 
studies  were  compared  to  the  average  analyzed  profiles  for  2002.  One  case  study  will 
examine  a  developing  circulation  (TS  24)  that  was  consistently  analyzed  and  forecast  to 
be  above  the  average  analyzed  values  for  all  2002  developing  vortices,  as  well  as  the  one 
model  “miss”  (TD  27)  that  occurred  in  2002. 

A  second  analysis  line  was  added  to  the  figures  for  the  two  case  studies  (Figures 
3.36  and  3.37).  This  heavy  blue  line  represents  the  average  analyzed  value  of  the 
selected  forecast  variable  for  the  case  study  circulation.  The  light  dashed  line 
corresponds  to  the  average  value  at  It)  for  the  selected  variable.  The  heavy  red  line 
represents  forecast  values  of  the  selected  forecast  variable  for  the  case  study  circulation 
for  the  time  indicated  in  the  title  of  each  forecast  panel  (24-,  48-,  72-,  and  96-h  in  the 
upper- left,  upper- right,  lower- left,  and  lower- right,  respectively).  As  in  Figure  3.5,  the 
heavy  black  line  represents  the  average  analyzed  value  (at  -i-OO)  of  the  selected  variable 
for  all  developing  circulations  in  2002,  and  the  light  green  whiskers  extending  from  this 
line  indicate  one  standard  deviation  of  variable  values  at  each  time  interval.  Box  plots 
were  not  plotted  for  the  forecast  vorticity  due  to  the  small  number  of  forecasts  available 
at  each  time  interval. 

1.  Over-forecast  Circulation  (TS  24) 

Tropical  Storm  24  during  September  2002  was  a  moderate  circulation  with 
maximum  sustained  surface  winds  of  55  kt  in  the  South  China  Sea.  The  average 
analyzed  850- mb  relative  vorticity  profile  for  TS  24  (Figure  3.36)  was  at  almost  all  times 
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greater  than  the  average  analyzed  vorticity  for  all  developing  circulations  in  2002.  The 
average  analyzed  vorticity  at  the  time  JTWC  issued  its  first  warning  (Fq)  for  TS  24  was 
8.2  X  10'^  s'\  which  is  more  than  one  standard  deviation  greater  (over  3  x  10'^  s'^  higher) 
than  the  average  analyzed  vorticity  at  F)  for  all  2002  developing  ■vortices.  Had  the 
averaged  analyzed  vorticity  value  of  5.0  x  10'^  s'*  been  used  as  an  indicator  of  formation, 
the  time  of  the  first  JTWC-issued  warning  would  have  been  about  60  h  earlier.  The  weak 
increase  in  vorticity  that  occurs  in  the  24  h  following  Fo  for  TS  24  represents  that  a  small 
amount  of  adjustment  occurred  as  a  result  of  the  synthetic  tropical  cyclone  observations 
that  were  inserted  within  12  h  of  Fq. 

The  deviation  between  the  analyzed  vorticity  for  all  developing  circulations  and 
the  analyzed  vorticity  for  TS  24  increases  with  increasing  time  after  Fo  for  TS  24.  In  the 
early  stages  of  the  development  of  TS  24,  the  analyzed  vorticity  is  very  similar  to  the 
analyzed  vorticity  for  all  developing  circulations  in  2002.  However,  this  similarity  does 
not  persist  more  than  24  h,  and  by  84  h  prior  to  Fb,  the  analyzed  vorticity  for  TS  24 
increased  almost  to  the  vorticity  threshold  for  all  developing  vortices.  That  the  analyzed 
vorticity  for  TS  24  continues  to  increase  after  this  point  partially  explains  why  the 
analyzed  vorticity  at  Fq  is  so  much  greater  than  the  average  analyzed  vorticity  at  Fq  for  all 
developing  vortices. 

The  24- h  forecasts  of  850- mb  relative  vorticity  (upper- left  panel  of  Figure  3.36) 
for  TS  24  fluctuate  between  over-  and  under- forecasts  throughout  the  120  h  before  and 
after  Fq.  The  vorticity  is  forecast  relatively  accurately  from  120  to  84  h  prior  to  Fq,  and 
again  from  36  h  to  about  12  h  prior  to  Fq.  Vorticity  is  over-forecast  (forecast  values 
exceed  analyzed  values)  between  those  two  periods  of  accurate  forecasting,  and  again 
from  roughly  12  h  prior  to  F)  until  24  h  after  F)-  The  sudden  decrease  in  forecast 
vorticity  from  12  h  after  F)  to  24  h  after  Fo  indicates  a  transition  from  forecasts  made 
prior  to  the  insertion  of  the  synthetic  tropical  cyclone  observations  to  those  made  after 
they  were  inserted.  Following  the  insertion  of  the  synthetic  tropical  cyclone 
observations,  the  forecasts  of  relative  vorticity  for  TS  24  are  consistently  less  than  the 
analyzed  vorticity.  This  tendency  of  NOGAPS  to  under- forecast  circulation  development 
after  the  synthetic  tropical  cyclone  observations  are  inserted  is  consistent  with  the  trend 
found  for  the  overall  sample  of  developing  storms. 
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Figure  3.36.  Average  analyzed  and  foreeast  850- mb  relative  vortieity  (x  10'^  s'^)  as  in 
Figure  3.5,  except  for  TS  24  (2002).  The  heavy  black  line  represents  analyzed  average 
8 50- mb  relative  vortieity  for  all  2002  developing  vortices  (at  -i-OO),  the  heavy  blue  line 
represents  average  analyzed  vortieity  (at  -i-OO)  for  TS  24,  and  the  heavy  red  line 
represents  average  forecast  850- mb  relative  vortieity  for  TS  24  for  the  indicated  forecast 
time.  Average  analyzed  8 50- mb  relative  vortieity  for  TS  24  at  Fb  (light  dashed  black 
line)  was  8.2  x  10'^  s ?o  for  developing  storms  in  2002  (light  solid  black  line)  is  5.0  x 
10-'  s-i. 


A  similar  trend  in  the  12  h  prior  to  Fq  is  apparent  in  the  48-h  forecasts  for  TS  24 
(upper- right  panel  of  Figure  3.36).  The  forecast  vortieity  at  Fo  greatly  exceeds  the 
analyzed  vortieity,  and  then  is  consistently  over-forecast  until  48  h  after  F).  Once  the 
synthetic  tropical  cyclone  observations  are  inserted,  at  the  forecast  verification  time 
corresponding  to  -i-48  h  in  the  upper- right  panel  of  Figure  3.36,  the  characteristics  of  the 
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vortex  in  the  analysis  are  corrected,  and  the  subsequent  vorticity  forecasts  are  less  than 
the  analyzed  value. 

The  reliability  of  the  forecast  vorticity  curves  in  the  extended  range  forecasts 
(bottom  panels  of  Figure  3.36)  are  questionable  given  the  very  small  number  of  forecasts 
available  at  each  of  these  forecast  times  and  time  intervals.  Data  are  missing  at  many  of 
the  time  intervals  in  the  extended  range  forecasts,  which  indicates  that  forecasts  verifying 
at  these  time  intervals  were  not  available,  or  possibly  that  the  circulation  was  not  forecast 
by  NOGAPS  at  these  times  (data  gaps). 

Statistical  significance  at  each  forecast  verification  time  could  not  be  calculated 
for  TS  24,  since  an  insufficient  number  of  cases  were  available  at  each  time  interval  and 
forecast  time.  An  ad  hoc  determination  of  statistical  significant  in  Figure  3.36  is 
determined  visually  by  noting  whether  or  not  the  forecast  value  lies  within  the  range  of 
analyzed  values  for  the  mean  curve  at  each  forecast  verification  time.  Most  of  the 
forecast  values  for  TS  24  are  greater  than  the  analysis  range,  which  suggests  statistical 
significance  using  this  rough  estimate.  However,  most  of  the  analyzed  vorticity  values 
for  TS  24  also  exceed  the  analysis  range  for  all  developing  circulations,  which  indicates 
that  TS  24  is  likely  one  of  those  circulations  that  represents  a  forecast  extreme  for  this 
year. 

2.  Under-forecast  Circulation  (TD  27) 

The  TCVTP  tracker  failed  to  detect  the  developing  circulation  that  became  TD  27 
prior  to  the  first  JTWC- warning  time.  This  error  represented  the  only  major  “miss”  of  the 
NOGAPS  analyses  during  1  May  -  31  October  2002.  TD  27  was  a  weak  tropical 
depression  with  maximum  sustained  surface  winds  of  only  30  kt  (2002  ATCR)  within  the 
Philippine  Sea. 

Analyzed  850- mb  vorticity  values  greater  than  1.5  x  10'^  s‘^  do  not  exist  for  TD 
27  until  after  the  first  warning  was  issued  by  JTWC  and  the  synthetic  tropical  cyclone 
observations  were  inserted  into  the  NOGAPS  analysis.  Consequently,  24-,  48-,  and  72-h 
forecasts  are  not  available  for  TD  27  until  after  the  synthetic  observations  were  inserted, 
and  no  120-h  forecasts  are  available  since  the  circulation  was  warned  on  for  only  2.5  days 
by  JTWC.  To  display  the  analyzed  and  forecast  relative  vorticity  (Figure  3.37)  for  TD 
27,  the  first  warning  time  had  to  be  shifted  to  12  h  after  Fo,  since  the  TCVTP  process  did 


83 


not  identify  the  circulation  in  the  NOGAPS  analysis  until  that  time.  Thus,  the  vorticity 
value  indicated  at  Fq  is  actually  the  analyzed  vorticity  value  at  Fq  +  12  h.  Similar 
displacement  of  the  forecast  curves  is  not  necessary.  The  apparent  phase  shift  in  the  start 
time  of  each  forecast  curve  in  Figure  3.37  is  due  to  the  lag  time  between  F)  and  the 
forecast  time,  i.e.,  the  24-h  vorticity  forecasts  (left  panel  of  Figure  3.37)  are  not  available 
until  24  h  after  Fo  since  no  forecasts  were  made  for  TD  27  until  Fq.  No  forecasts  were 
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Figure  3.37.  Average  analyzed  and  forecast  850- mb  relative  vorticity  (x  10'^  s'^)  as  in 
Figure  3.5,  except  for  TD  27  (2002).  The  heavy  black  line  represents  analyzed  average 
8 50- mb  relative  vorticity  for  all  2002  developing  vortices  (at  -rOO),  the  heavy  blue  line 
represents  average  analyzed  vorticity  (at  -i-OO)  for  TD  27,  and  the  heavy  red  line 
represents  average  forecast  850-mb  relative  vorticity  for  TD  27  for  the  indicated  forecast 
time.  Average  analyzed  850-mb  relative  vorticity  for  TD  27  at  Fb  (light  dashed  black 
line)  was  3.6  ?o  for  developing  storms  in  2002  (light  solid  black  line)  is  5.0  x  10'^  s 


available  at  120  h,  and  so  few  forecasts  were  available  at  72  h  that  neither  of  those 
forecast  times  are  shown  in  Figure  3.37. 

Examination  of  the  analyzed  850-mb  relative  vorticity  for  TD  27  (Figure  3.37) 
indicates  that  after  a  weak  development  phase  that  lasted  24  h,  the  circulation  began  to 
weaken.  At  all  times,  the  analyzed  vorticity  for  this  circulation  was  substantially  less 
than  the  analyzed  vorticity  for  all  developing  circulations,  and  never  reached  the  vorticity 
threshold  of  5.0  x  10'^  s'^  for  developing  circulations.  If  surpassing  this  threshold  in  the 
NOGAPS  analyses  were  the  only  criterion  for  issuing  a  first  warning,  it  would  not  be 
surprising  for  JTWC  to  have  missed  this  warning.  That  is,  evidence  was  lacking  in  the 
synoptic  analyses,  and  JTWC  had  to  rely  on  satellite  imagery  to  issue  its  first  warning. 
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After  the  brief  phase  of  increasing  vorticity,  the  analyzed  vorticity  for  TD  27  decreased 
as  the  analyzed  vorticity  for  all  developing  storms  continued  to  increase,  and  never 
underwent  a  second  pulse  in  development. 

Most  forecasts  of  relative  vorticity  that  are  available  for  TD  27  are  all  markedly  in 
excess  of  the  analyzed  vorticity  at  the  verifying  time  of  the  forecast,  especially  36  h  after 
Fq.  The  exception  to  this  is  the  24-h  forecast  (left  panel  of  Figure  3.37),  which  does  show 
a  transition  from  under-  to  over- forecast  vorticity  values  24  h  after  the  first  JTWC 
warning  time.  This  transition  is  in  direct  contrast  to  the  shift  in  forecast  values  relative  to 
the  insertion  of  the  synthetic  tropical  cyclone  observations,  but  it  does  occur  at  the  same 
time  that  any  sort  of  forecast  transition  is  expected.  Additionally,  the  time  of  the  forecast 
maximum  vorticity  increases  with  increasing  forecast  range,  and  with  a  small  decrease  in 
magnitude  with  increasing  forecast  range. 

The  apparent  phase  shift  that  occurs  between  the  analyzed  and  forecast  vorticity 
in  Figure  3.37  is  attributed  to  an  inaccurate  representation  of  the  large-scale  environment, 
and  also  that  the  vortex  grew  only  to  a  small,  weak  tropical  depression.  The  average 
analyzed  vorticity  at  the  first  JTWC-waming  time  was  3.6  x  10'^  s  ',  almost  1.5  x  10'^  s'* 
lower  than  analyzed  for  the  full  developer  data  set,  and  only  0.3  x  10'^  s'*  greater  than  the 
threshold  for  non- developing  vortices.  That  the  average  analyzed  vorticity  value  at  It) 
was  less  than  4  x  10'^  s'*  indicates,  from  previous  discussion,  that  shear  should  be 
westerly,  which  is  not  favorable  for  formation.  For  formation  to  occur,  the  shear  must  be 
less  westerly,  and  perhaps  very  close  to  zero.  It  is  possible  that  continued  westerly  shear 
existed  and  ventilated  the  developing  vortex,  and  thus  further  development  of  TD  27  was 
prevented. 

The  failure  of  NOGAPS  analysis  to  include  the  developing  TD  27  circulation 
until  it  was  artificially  inserted  into  the  analysis  suggests  that  the  formation  environment 
of  TD  27  may  have  been  substantially  disorganized,  very  weak,  or  misrepresented  in  the 
model  analysis  due  to  a  lack  of  observations.  The  NOGAPS  analysis  only  has  a  weak 
(1010  mb)  low  with  a  closed  isobar,  but  has  no  closed  850- mb  relative  vorticity  contour 
of  at  least  1.5  x  10'^  s'*.  Even  after  the  insertion  of  synthetic  observations,  the  circulation 
appears  weakly  in  the  NOGAPS  forecasts.  Perhaps  access  to  other  global  and  regional 
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numerical  models  would  allow  the  TCVTP  program  to  identify  such  weak,  disorganized 
circulations. 


E.  COMPARISON  OF  2003  VORTICES  TO  2002  VORTICES 

Since  2003  was  a  year  following  a  warm  El  Nino  (EN)  event,  tropical  cyclone 

activity  in  the  western  North  Pacific  was  forecast  to  be  below  normal,  and  above -normal 
tropical  cyclone  activity  was  forecast  in  the  South  China  Sea  (Chan  2003).  The 
anticipated  cold  phase  of  the  El  Nino-Southem  Oscillation  (ENSO),  called  a  Ea  Nina 
event,  was  much  weaker  than  expected.  Although  the  observed  flow  patterns  and  tropical 
cyclone  formations  in  the  SCS  were  consistent  with  a  Ea  Nina  year  (Chan  2003),  below- 
normal  formation  numbers  occurred  in  the  SCS  due  to  an  anomalous  anticyclonic 
circulation  over  the  SCS  from  June  -  October  2003.  This  anomalous  circulation  is  linked 
to  interdecadal  variations  in  the  strength  of  the  subtropical  high  (Chan  2003).  Although 
numbered  storms  were  fewer  in  2003,  the  total  number  of  tracked  non- developing 
vortices  in  the  NOGAPS  model  was  similar  to  2002  (see  Tables  3.1  and  3.9). 

I.  Developing  Vortices 

The  western  North  Pacific  Ocean  basin  was  less  active  during  2003  than  it  had 
been  during  2002,  with  20  circulations  numbered  by  JTWC  from  1  May  -  31  October 
2003  compared  to  24  in  2002.  As  in  2002,  the  NOGAPS  failed  to  identify  several 
circulations  prior  to  insertion  of  the  synthetic  tropical  cyclone  observations.  It  is  possible 
that  the  NOGAPS  model  was  unable  to  distinguish  these  two  weak  circulations  from  the 


Table  3.9.  Circulations  tracked  in  NOGAPS  analyses  from  1  May  -  31  Oct  2003 
numbered  by  the  JTWC  (column  1),  non- developing  circulations  meeting  the  minimum 
duration  criterion  of  at  least  24  h  (column  2),  those  non- developing  circulations  not 
meeting  the  minimum  duration  criterion  (column  3),  and  the  non- developing  vortices  that 
were  forecast  to  occur,  but  were  never  analyzed  in  NOGAPS  (column  4). 


Eormation  Subregion 

Numbered 

TCs 

ND  >  24  h 
(NDl) 

ND  <  24  h 
(NDOE) 

Model 

Ealse 

Alarms 

South  China  Sea  (105E- 124E) 

4 

35 

18 

4 

Philippine  Sea  (125E-159E) 

14 

64 

28 

18 

East  Monsoon  Trough  (160E-180E) 

2 

6 

16 

5 

Total 

20 

105 

62 

27 
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background  environment  due  to  insufficient  organization  or  relative  vorticity.  These 
circulations,  TD  22  and  TS  23,  are  thus  excluded  from  Table  3.10  and  also  from  the 
discussion  below.  Numbered  storms  again  formed  preferentially  in  the  PS  subregion  (14 
of  20),  and  percentages  of  developments  in  the  SCS  (4)  and  EMT  (2)  were  consistent 
with  2002  (Table  3.10).  It  is  worth  noting  that  formation  within  a  specific  subregion,  as 
noted  in  Table  3.10,  is  based  on  the  first  analysis  time  the  TCVTP  classified  an  ellipse  for 
each  circulation.  Given  this,  it  is  likely  that  two  of  the  developing  circulations  identified 
as  developing  in  the  PS  subregion  (near  126°E)  may  have  first  been  warned  on  in  the  SCS 
subregion. 

The  average  analyzed  850- mb  relative  vorticity  for  developing  vortices  in  2003 
(Eigure  3.38)  was  4.8  x  10'^  s'\  which  is  only  slightly  less  than  ?o  for  2002  (5.0  x  10'^  s'*). 
The  standard  deviation  at  each  time  interval  for  the  analyzed  vorticity  in  2003  was  also 
notably  larger,  especially  at  the  later  time  intervals,  than  during  2002.  This  greater 
variation  is  associated  with  greater  variability  of  the  developing  circulations,  but  may 
also  have  been  due  to  the  NOGAPS  transition  in  mid- October  2003  from  the  Multivariate 
Optimum  Interpolation  (MVOI)  data  assimilation  scheme  to  the  new  Navy  Data 
Assimilation  System  (NAVDAS). 

Goerss  et  al.  (2003)  compared  the  performance  of  simultaneous  runs  of  the 
NOGAPS  using  the  NAVDAS  scheme  and  the  MVOI  scheme  during  two  periods  in 
2002,  and  during  the  beta- test  during  the  spring  of  2003.  They  found  that  tropical 
cyclone  track  forecasts  to  72  h  initiated  from  the  two  data  assimilation  schemes  had  only 
small  differences.  While  performance  was  worse  for  96-  and  120-h  track  forecasts,  these 
errors  were  based  mainly  on  the  variability  of  the  individual  developing  circulations. 
Given  this,  it  is  unlikely  that  the  variation  between  the  2002  and  2003  vorticity  analyses 


Table  3.10.  Summary  of  tropical  cyclone  activity  in  the  western  North  Pacific  by 
subregion  for  1  May  -  31  October  2003.  Intensities  corresponding  to  each  column  are 
identified  in  Table  3.2. 


Eormation  Subregion 

TD 

TS 

TY 

STY 

Total 

South  China  Sea  (105E- 124E) 

2 

2 

0 

0 

4 

Philippine  Sea  (125E-159E) 

0 

3 

9 

2 

14 

East  Monsoon  Trough  (160E-180E) 

0 

0 

2 

0 

2 

Total 

2 

5 

11 

2 

20 
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is  due  to  the  different  data  assimilation  schemes.  It  is  also  noted  that  NAVDAS  was  only 
operationally  used  for  the  last  three  months  of  2003  and  during  only  one  month  of  this 
study. 

The  slope  of  the  analyzed  850-mb  relative  vorticity  for  developing  circulations  in 
2003  (Figure  3.38)  was  steeper  than  the  slope  of  the  analyzed  vorticity  for  the  2002 
developing  vortices  (Figure  3.5),  which  possibly  indicates  that  intensification  during 
2003  occurred  in  spurts  rather  than  gradually.  The  increase  in  the  slope  of  the  analyzed 
vorticity  in  the  first  12  h  after  Fq  in  Figure  3.38  is  substantially  steeper  than  in  2002.  This 


Fq  2003  WPAC  NOGAPS  DV  at  +24 


Fq  2003  WPAC  NOGAPS  DV  at  +48 


Forecast  Verification  Time 


Fq  2003  WPAC  NOGAPS  DV  at  +72 


Fq  2003  WPAC  NOGAPS  DV  at  +96 


Forecast  Verification  Time 


Figure  3.38.  Average  analyzed  and  forecast  850-mb  relative  vorticity  (x  10'^  s'^)  as  in 
Figure  3.5,  except  for  developing  vortices  in  2003  relative  to  F).  Heavy  black  line 
represents  analyzed  average  850-mb  relative  vorticity  for  developing  vortices  (at  -i-OO), 
and  red  line  represents  average  forecast  850-mb  relative  vorticity  for  the  indicated 
forecast  time.  Average  analyzed  850-mb  relative  vorticity  for  2003  developing  vortices 
at  Fo  (light  dashed  black  line)  was  4.8  x  10'^  s Light  solid  black  line  represents  average 
analyzed  vorticity  at  Fq  for  2002  developing  vortices  (5.0  x  10'^  s"^). 
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dramatic  slope  increase  corresponds  to  the  insertion  of  the  synthetic  tropical  cyclone 
observations,  as  discussed  previously.  That  the  analyzed  vorticity  is  so  distinctly 
different  on  either  side  of  Fb  indicates  that  developing  vortex  was  not  being  accurately 
represented  in  the  model  analyses  prior  to  the  insertion  of  the  synthetic  tropical  cyclone 
observations. 

The  forecast  curves  of  850- mb  relative  vorticity  for  developing  vortices  h  2003 
(Figure  3.38)  were  dramatically  less  than  the  analyzed  vorticity  at  almost  all  forecast 
verification  times  at  all  forecast  intervals.  Notably,  while  the  deviations  between 
NOGAPS  analyses  and  forecasts  of  relative  vorticity  at  Fq  were  not  statistically 
significant  for  the  circulations  that  developed  in  2002,  the  differences  were  substantially 
larger  for  the  24-  and  48-h  forecasts  (Figure  3.38)  during  2003.  In  the  2002  data  set,  the 
forecasts  of  vorticity  at  It)  were  consistently  at  or  greater  than  the  analyzed  vorticity, 
whereas  in  the  2003  data  set,  forecast  vorticity  was  consistently  less  than  the  analyzed 
value  at  Fq. 

The  under- forecasting  of  relative  vorticity  at  Fq  occurred  partially  due  to  the 
earlier  transition  from  over-  to  under- forecast  values  of  vorticity  in  the  2003  NOGAPS 
forecasts  than  in  2002.  This  transition  occurred  in  the  2003  sample  between  60  and  48  h 
prior  to  F)  for  the  shorter-range  forecasts,  which  is  much  earlier  than  during  2002  in 
which  the  transition  from  over-  to  under- forecasts  did  not  occur  until  after  F).  This 
difference  in  the  forecasts  is  not  attributed  to  the  data  assimilation  scheme,  or  the 
insertion  of  the  synthetic  tropical  cyclone  observations,  but  instead  to  the  inability  of  the 
NOGAPS  model  to  correctly  represent  the  early  formation  environment  of  the  developing 
vortices  during  2003. 

The  insertion  of  the  synthetic  tropical  cyclone  observations  in  the  NOGAPS 
analyses  is  less  apparent  in  the  2003  forecasts  (Figure  3.38)  than  in  the  2002  forecasts 
(Figure  3.5).  Although  the  data  assimilation  scheme  was  changed  during  October  2003, 
the  process  by  which  synthetic  tropical  cyclone  observations  were  inserted  into  the  model 
analysis  did  not  change  (B.  Strahl,  FNMOC,  personal  communication).  Given  this,  the 
expected  flattening  of  the  forecast  vorticity  slope  between  F)  and  the  forecast  time  in 
Figure  3.5  was  not  noted  in  Figure  3.38.  Instead,  the  slope  of  forecast  vorticities 
generally  increased  steadily  between  Fq  and  the  time  corresponding  to  the  forecast 
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interval  (indicated  in  the  title  of  each  panel),  and  then  increased  dramatically  after  that. 
The  forecast  vorticity  slope  between  Fb  and  the  forecast  time  appeared  to  increase  with 
increasing  forecast  range  (a  shallower  slope  appears  in  the  mean  24-  h  forecast  than  in  the 
mean  72-h  forecast).  The  apparent  growth  between  Fo  and  the  forecast  time  is  attributed 
to  continued  development  of  the  circulation  in  the  NOGAPS  forecasts. 

Additionally,  the  forecast  spread  increased  with  increasing  forecast  range,  and 
tended  more  toward  lower  values  than  analyzed,  rather  than  the  expected  over- forecasts. 
In  many  instances,  the  forecast  maximum  for  a  time  interval  identified  at  the  top  of  the 
box  plot  was  at  or  below  the  analyzed  value.  The  forecast  extrema  during  2003  were  also 
clustered  around  It),  and  included  both  over-  and  under- forecast  values.  The  under¬ 
forecast  extrema  occurred  only  in  the  shorter- range  forecasts  (<  48  h).  Combining  the 
extrema,  forecast  spread,  and  consistent  under- forecasting  evident  in  Figure  3.38,  it  is 
concluded  that  the  NOGAPS  model  did  not  accurately  characterize  either  the  developing 
tropical  circulations  or  the  background  environment. 

The  average  analyzed  deep  layer  wind  shear  for  vortices  developing  in  2003 
(Figure  3.39)  has  a  slope  similar  to  the  average  analyzed  deep  layer  wind  shear  for  the 
developing  vortices  in  2002,  although  the  average  analyzed  deep  shear  at  Fo  in  2003  was 
-1.6  m  s  ',  which  is  slightly  less  easterly  than  the  analyzed  deep  shear  Fo  in  2002  (-1.8  m 
s  ').  As  in  the  2002  data  set,  the  transition  from  westerly  to  easterly  shear  occurs  when 
the  analyzed  vorticity  exceeds  5.0  x  10'^  s  '.  In  2003,  this  transition  occurs  between 
roughly  24  and  12  h  prior  to  Fo. 

The  weaker  analyzed  wind  shear  during  2003  is  likely  a  result  of  the  \^eak  La 
Nina  conditions  that  occurred  in  the  western  North  Pacific  during  that  year  (Chan  2003). 
While  the  model  analyses  and  forecasts  reflect  this,  the  range  of  the  forecast  deep  shear  at 
each  forecast  verification  time  is  significantly  smaller  during  2003  (Figure  3.39)  than 
they  were  in  2002  (Figure  3.14).  While  impacted  by  the  smaller  sample  size  in  2003,  this 
decrease  in  forecast  range  at  each  forecast  verification  time  is  likely  a  result  of 
modifications  made  to  the  NOGAPS  Emanuel  cumulus  parameterization  scheme  during 
2003  (Hogan  et  al.  2004).  One  of  the  key  modifications  made  to  the  Emanuel  cumulus 
parameterization  scheme  was  to  increase  the  amount  of  convective  momentum  transport, 
which  then  brings  more  easterly  momentum  downward  and  is  countered  by  an  increase  in 
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Figure  3.39.  Average  analyzed  and  forecast  deep  layer  (200-850  mb)  wind  shear  (m  s'^) 
as  in  Figure  3.5,  except  for  developing  vortices  in  2003  relative  to  Fq.  Heavy  black  line 
represents  analyzed  average  deep  layer  wind  shear  for  developing  vortices  (at  +00),  and 
red  line  represents  average  forecast  deep  layer  wind  shear  for  the  indicated  forecast  time. 
Average  analyzed  deep  layer  wind  shear  for  2003  developing  vortices  at  F)  (light  solid 
black  line)  was  -1.6  m  s Light  solid  black  line  represents  average  analyzed  shear 
threshold  for  2002  developing  vortices  (-1.8  m  s'^). 


the  vertical  turbulent  mixing  by  the  parameterization  of  the  planetary  boundary  layer 
(Hogan  et  al.  2004).  The  end  result  of  this  increase  in  convective  momentum  transport  is 
a  decrease  in  the  magnitude  of  the  low-level  winds,  which  thus  decreases  the  amount  of 
easterly  deep  layer  wind  shear.  Thus,  not  only  are  deep  layer  wind  shear  forecast  ranges 
smaller  at  each  verification  time  than  in  the  2002  NOGAPS  forecasts,  but  they  are  also  on 
average  weaker  or  closer  to  zero. 


91 


NOGAPS  forecasts  of  deep  shear  were  remarkably  close  to  the  analyzed  deep 
shear  values,  especially  in  the  24-h  forecasts  prior  to  Fq  (upper- left  panel  in  Figure  3.39). 
The  24-h  forecast  closely  resembles  the  analyzed  deep  shear  for  2002  (Figure  3.14),  and 
relatively  accurately  forecasts  the  transition  from  westerly  to  easterly  deep  layer  shear  by 
12  h  prior  to  Fq.  Although  the  time  of  the  occurrence  of  the  maximum  easterly  shear  is 
inaccurately  forecast  at  24  h,  and  correctly  forecast  at  48  h,  the  forecast  maximum 
magnitude  is  within  tolerable  error  limits. 

The  range  of  forecast  values  in  Figure  3.39,  which  is  larger  farther  away  from  Fq, 
is  not  statistically  significant  in  the  2003  forecasts,  whereas  it  is  at  a  few  time  intervals 
prior  to  Fq  in  2002.  This  improvement  in  the  forecast  accuracy  of  2003  is  attributed  to  a 
smaller  range  of  analyzed  deep  layer  shear  values  in  the  developing  circulations, 
especially  within  24  h  of  Fq. 

2.  Non- developing  Vortices 

The  average  analyzed  vorticity  for  non- developing  vortices  in  2003  (Figure  3.40) 
was  similar  to  the  analyzed  vorticity  for  non- developing  vortices  in  2002  (Figure  3.18). 
The  average  analyzed  vorticity  at  No  in  2003  for  this  data  set  was  3.0  x  10"^  s'\  which  is 
only  slightly  smaller  than  the  average  analyzed  value  at  No  in  2002  (3.3  x  10'^  s'^).  This 
smaller  average  analyzed  vorticity  is  attributed  to  the  weaker  circulations  that  formed  in 
2003  when  compared  to  2002. 

In  direct  contrast  to  the  developing  vortices  for  2003,  forecasts  for  the  non¬ 
developing  vortices  are  almost  entirely  in  excess  of  the  analyzed  relative  vorticity.  While 
these  over- forecast  errors  in  2003  are  consistent  in  deviation  with  the  over- forecast  in 
2002,  the  smooth  forecast  increase  and  subsequent  decrease  relative  to  No  is  not  present 
in  2003.  Considerably  larger  fluctuations  exist  in  the  longer- range  forecasts  (bottom 
panels  of  Figure  3.40)  for  2003  than  2002,  which  may  be  partially  due  to  the  smaller 
sample  size  for  2003. 

Consistent  with  the  over- forecast  vorticity  at  almost  all  forecast  times,  the 
extrema  for  this  non- developing  data  set  all  exceed  the  analyzed  vorticity,  which 
indicates  a  tendency  for  the  NOGAPS  model  to  over- forecast  tropical  circulation 
development.  This  tendency  among  non- developing  vortices  is  in  direct  contrast  to  the 
tendency  towards  under- forecasting  for  the  developing  vortices  in  2003  (Figure  3.38).  In 
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Figure  3.40.  Average  analyzed  and  forecast  850- mb  relative  vorticity  (x  10'^  s'^)  as  in 
Figure  3.18,  except  for  non- developing  vortices  in  2003  relative  to  No.  Heavy  black  line 
represents  analyzed  average  850-mb  relative  vorticity  for  non- developing  vortices  (at 
-1-00),  and  red  line  represents  average  forecast  850-mb  relative  vorticity  for  the  indicated 
forecast  time.  Average  analyzed  850-mb  relative  vorticity  for  2003  non- developing 
vortices  at  No  (light  solid  black  line)  was  3.0  x  10'^  s Light  solid  black  line  represents 
?o  threshold  for  2002  developing  vortices  (5.0  x  10"^  s'^). 


the  first  36  h  of  the  life  cycle  of  the  developing  vortices  shown  in  Figure  3.38  (2.5  x  10'^ 
s  '),  the  NOGAPS  vorticity  forecasts  tend  to  exceed  the  analyzed  value. 

The  24- h  forecasts  of  850-mb  relative  vorticity  (upper- left  panel  of  Figure  3.40) 
most  closely  match  the  analyzed  relative  vorticity,  with  only  minor  deviations  in  the 
hours  following  Fq.  This  trend  for  accurate  forecasting  prior  to  Fq,  and  increasing 
forecast  deviation  from  the  analyzed  value  increases  with  increasing  forecast  range.  In 
the  48-,  72-,  and  96-h  forecasts,  this  deviation  becomes  statistically  significant  12  h  after 
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Fo,  as  forecasts  of  850-mb  relative  vorticity  increase  and  the  analyzed  values  steadily 
decrease.  The  2002  trend  for  increasing  forecast  time  of  the  vorticity  maximum  with 
increased  forecast  range  is  not  present  in  the  2003  forecasts,  which  could  indicate  that 
NOGAPS  better  resolved  the  non- developing  circulations  in  2003  than  it  did  in  2002. 

The  analyzed  deep  layer  wind  shear  (Figure  3.41)  for  these  non- developing  cases 
fluctuates  more  in  magnitude  in  2003  than  it  did  in  2002  (Figure  3.23),  potentially 
because  of  the  smaller  number  of  cases  in  the  2003  data  set.  Notice  that  the  vertical  shear 
is  even  stronger  westerly  during  2003  than  it  was  during  2002,  which  would  be  consistent 
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Figure  3.41.  Average  analyzed  and  forecast  deep  layer  (200-850  mb)  wind  shear  (m  s'^) 
as  in  Figure  3.18,  except  for  non- developing  vortices  in  2003  relative  to  Nq.  Heavy  black 
line  represents  analyzed  average  deep  layer  wind  shear  for  non- developing  vortices  (at 
+00),  and  red  line  represents  average  forecast  deep  layer  wind  shear  for  the  indicated 
forecast  time.  Average  analyzed  deep  layer  wind  shear  for  2003  non- developing  vortices 
at  No  (light  solid  black  line)  was  5.5  m  s  ^  Light  solid  black  line  represents  average 
analyzed  shear  threshold  for  2002  developing  vortices  (-1.8  m  s'*). 
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with  less  favorable  formation  conditions  during  2003.  This  increase  in  westerly  wind 
shear  is  likely  due  to  the  increased  westerlies  aloft  associated  with  the  weak  La  Nina 
event,  or  cold  phase  of  the  El  Nino- Southern  Oscillation,  as  noted  above  by  Chan  (2003). 

The  transition  to  easterly  analyzed  shear  at  -i-84  h  is  suspect  due  to  the  small 
number  of  cases  at  >  72  h  after  N)-  The  analyzed  wind  shear  remains  positive  and 
westerly  for  this  subset  of  non- developers,  as  in  2002,  which  is  a  primary  indicator  that 
these  vortices  will  not  likely  undergo  significant  further  development. 

Forecast  deep  layer  wind  shear  is  at  almost  all  times  weaker  (less  positive  and  thus  less 
westerly)  than  the  analyzed  shear,  which  is  consistent  with  the  NOGAPS  tendency  to 
over- forecast  tropical  circulations.  No  significant  trend  exists  in  forecast  range  compared 
to  the  forecast  time,  with  too  large  shear  errors  nearly  as  common  as  too  small  shear 
errors  at  the  longer  forecast  ranges.  However,  this  may  be  due  to  the  small  sample  size  at 
extended  forecast  intervals. 

3.  False  Alarms 

Surprisingly,  the  number  of  model  false  alarms  (27)  in  2003  increased 
significantly  from  2002  (15).  Formation  of  these  model  false  alarms  (Figure  3.42) 
occurred  preferentially  in  the  Philippine  Sea,  as  in  2002.  These  erroneous  forecasts  are 
likely  associated  with  spurious  convective  events  within  the  monsoon  trough. 


Figure  3.42.  Focations  of  model  false  alarms  in  the  western  North  Pacific  during  1  May 
-31  October  2003.  Circles  indicate  successive  12-h  analyses  of  each  identified  model 
false  alarm.  The  size  of  each  circle  is  related  to  the  analyzed  intensity  of  the  circulation 
with  a  larger  circle  indicating  a  more  intense  circulation. 
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However,  the  number  of  Group  2- type  false  alarms  decreased  to  fewer  than  five 
(from  50),  which  is  attributed  to  modifications  made  to  the  Emanuel  cumulus 
parameterization  scheme  late  in  2002  that  preferentially  decreased  the  easterly  deep  layer 
wind  shear  by  increasing  the  convective  momentum  transport.  The  2003  model  false 
alarms  were  not  similarly  decreased  by  this  change  to  the  cumulus  parameterization  since 
unlike  the  Group  2  non- developing  vortices,  they  were  never  analyzed  and  also  did  not 
have  excessive  vorticity  forecasts. 
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IV.  SUMMARY  AND  CONCLUSIONS 


The  performance  of  the  NOGAPS  model  in  predicting  the  formation  of  tropical 
circulations  in  the  western  North  Pacific  is  assessed  by  tracking  254  circulations  from  1 
May  -  31  October  2002.  The  Tropical  Cyclone  Vortex  Tracking  Program  (TCVTP) 
developed  by  Professor  Patrick  Harr  is  used  to  identify  circulations  in  the  NOGAPS  850- 
mb  relative  vorticity  analysis  and  forecast  fields  that  have  a  value  of  at  least  1.5  x  10“^  s 
An  ellipse  is  fit  to  the  outer  closed  vorticity  contour  to  define  the  size  of  the  circulation. 
Identified  circulations  are  then  matched  when  distance  and  movement  criteria  are  met  to 
define  tracks  in  the  analyses  and  forecasts. 

Any  circulations  that  formed  north  of  30°N  were  not  tracked  by  the  TCVTP. 
Additionally,  circulations  that  formed  east  of  180°E  were  also  excluded  from  the  analysis 
data  set.  The  remaining  circulations  were  categorized  according  to  whether  they  formed 
within  the  South  China  Sea  (SCS),  Philippine  Sea  (PS),  east  of  the  monsoon  trough 
(EMT).  A  minimum  duration  of  24  h  was  imposed  for  the  analyzed  vortex  in  the 
analyses.  Developing  circulations  were  analyzed  relative  to  their  first  JTWC-waming 
time,  and  relative  to  the  first  Best-Track  time.  Those  same  circulations  were  also 
analyzed  based  on  whether  they  formed  before,  during,  or  after  the  TCEA- specified 
formation  window.  Additionally,  non- developing  vortices  and  both  model  and  Tropical 
Cyclone  Eormation  Alert  (TCEA)  false  alarms  were  examined. 

Tropical  circulations  within  the  Philippine  Sea  subregion  accounted  for  46.5%  of 
the  total  tracked  vortices,  and  43.5%  of  the  JTWC- numbered  circulations  during  2002. 
Of  the  59  potential  seedlings  in  the  PS  subregion,  the  ten  vortices  that  developed  to 
warning  strength  represented  only  a  16.9%  formation  rate.  Eormations  at  the  east  end  of 
the  monsoon  trough  accounted  for  only  22.8%  of  the  total  tracked  vortices,  but  led  to 
34.8%  of  the  numbered  storms.  The  eight  numbered  storms  that  formed  in  the  EMT 
subregion  represented  27.6%  of  the  29  circulations  that  were  identified  within  that 
subregion.  Circulations  that  developed  within  the  SCS  accounted  for  30.7%  of  the  total 
tracked  vortices,  and  21.7%  of  the  numbered  storms.  The  formation  rate  was  the  smallest 
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in  this  subregion,  with  only  12.8%  of  the  potential  seedlings  developing  to  at  least 
Tropical  Depression  intensity. 

A  total  of  27  TCFAs  were  issued  for  the  23  storms  that  developed  during  1  May  - 
31  October  2002.  TCFAs  were  not  issued  for  two  of  the  developing  storms,  and  multiple 
TCFAs  were  issued  for  six  developing  storms.  Two  of  these  twice- issued  TCFAs  were 
not  available  for  analysis.  Formation  occurred  during  the  TCFA- specified  formation 
window  for  33%  of  the  developing  storms.  Of  the  remaining  16  that  failed  to  verify 
within  the  specified  window,  61%  verified  prior  to  the  TCFA- specified  formation 
window.  About  19%  of  the  developing  storms  formed  after  the  specified  formation 
window.  In  general,  when  formation  did  not  occur  by  the  end  of  the  formation  window,  a 
second  TCFA  was  issued  to  correct  the  formation  window  and  location,  and  formation 
typically  did  occur  within  the  second  specified  formation  window.  Six  false  alarm 
TCFAs  were  issued  for  circulations  that  were  never  warned  on  by  JTWC.  Due  to  the 
small  number  of  cases  available  when  analyzing  developing  storms  based  on  formation 
time  relative  to  the  TCFA  formatbn  window,  significant  systematic  NOGAPS  model 
trends  cannot  be  identified.  Before  any  conclusions  can  be  drawn  regarding  the 
performance  of  NOGAPS  in  predicting  formation  relative  to  the  TCFA- specified 
formation  window,  especially  for  those  circulations  that  formed  late,  a  study  that  includes 
more  circulations  meeting  these  criteria  is  necessary.  Analysis  based  on  formation 
region,  rather  than  time,  will  be  the  focal  point  of  further  analysis. 

A.  DEVELOPING  VORTICES 

Formation  time  for  the  circulations  that  were  later  numbered  and  warned  on  by 
JTWC  was  taken  to  be  the  time  the  first  warning  was  issued.  This  time  (designated  Fo) 
was  identified  for  each  developing  storm  from  the  2002  Annual  Tropical  Cyclone  Report 
(ATCR)  published  by  JTWC.  A  second  formation  time  called  the  first  Best-Track  time, 
(designated  F*o)  was  defined  as  the  time  the  developing  circulation  was  first  apparent  in 
the  model  analysis,  which  is  determined  by  JTWC  during  post- analysis,  and  also 
published  in  the  2002  ATCR.  The  first  Best- Track  time  always  preceded  the  first 
warning  time,  normally  by  at  least  12  hours. 
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The  timing  error  (?  ti)  between  when  a  forecast  curve  first  crossed  the  analyzed 
850-mb  relative  vorticity  threshold  (5.0  x  10'^  s'*)  and  the  first  JTWC-waming  time  was 
calculated.  On  average  the  forecast  curve  crossed  the  vorticity  threshold  16.7  h  prior  to 
Fo  for  all  forecasts,  but  only  6  h  prior  to  Fo  for  the  short-term  (<  48  h)  forecasts.  A  timing 
error  less  than  6  h  is  not  possible  due  to  the  resolution  of  the  data.  This  relatively  small 
(negative)  timing  error  indicates  that  the  NOGAPS  model  devebps  circulations  slightly 
faster  than  they  develop  in  the  environment.  Analysis  of  the  analyzed  and  forecast 
vorticity  curves  in  Figure  3.5  (upper  panels)  indicates  that  in  the  36  h  leading  up  to  the 
first  JTWC-waming  time,  the  NOGAPS  model  forecasts  of  vorticity  were  only  slightly 
greater  than  the  analyzed  vorticity  at  those  forecast  verification  times.  At  the  longer 
forecast  ranges,  the  forecast  vorticity  exceeded  the  analyzed  vorticity  by  less  than  1  x  10'^ 
s'*,  which  indicates  an  over- forecast,  but  rot  an  excessive  one.  Given  these  relatively 
accurate  forecasts,  the  NOGAPS  model  appears  to  be  handling  the  physical 
parameterizations  well  in  the  formation  environment. 

Threshold  values  at  these  formation  time  definitions  were  determined  for  five 
selected  forecast  variables:  850-mb  relative  vorticity,  deep  layer  (200-850  mb)  wind 
shear,  sea- level  pressure,  925- mb  wind  speed,  and  vapor  pressure.  Additionally,  the  size 
of  the  circulation  was  analyzed  to  determine  the  accuracy  of  the  NOGAPS  representation 
of  the  circulation,  but  was  not  treated  as  a  forecast  variable.  The  threshold  values 
corresponding  to  the  average  value  of  each  variable  at  the  first  JTWC-waming  time  is 
given  in  the  second  column  of  Table  4.1.  A  second  set  of  threshold  values  corresponded 
to  the  value  of  each  variable  at  the  first  Best- Track  time  is  given  in  the  third  column  of 
Table  4.1. 


Table  4.1.  Summary  of  threshold  values  relative  to  Fq  and  F*o.  Values  correspond  to  the 
average  analyzed  value  of  the  variable  at  the  reference  time. 


Forecast  Variable 

Analyzed  Threshold  Values 

Fo 

F*o 

850-mb  relative  vorticity 

Deep  layer  (200-850  mb)  wind  shear 

-1.8  ms'' 

-3.8  m  s'* 

Sea- level  pressure 

1006.9  mb 

1006.6  mb 

925-mb  wind  speed 

9.5  m  s'* 

Vapor  pressure  (500-700  mb  average) 

3.6  Pa 

3.4  Pa 
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The  strongest  trends  relative  to  the  formation  time  were  found  in  the  850- mb 
relative  vorticity  and  the  deep  layer  wind  shear.  When  the  analyzed  average  8 50- mb 
relative  vorticity  value  was  approximately  4  x  10'^  s'^  the  sign  of  the  analyzed  deep  layer 
wind  shear  tended  to  change  from  westerly  (positive)  to  easterly  (negative).  Classic 
research  (Gray  1968)  indicates  that  a  minimum  in  vertical  wind  shear  over  the  center  of 
the  potential  tropical  seedling  is  necessary  for  formation  to  occur,  along  with  easterly 
shear  to  the  north,  and  westerly  shear  to  the  south.  Some  research  model  studies  (Tuleya 
and  Kurihara  1981)  indicate  that  zero  shear  above  the  developing  circulation  is  less 
desirable  than  slightly  negative  (easterly)  shear,  which  is  consistent  with  this  research. 
That  is,  the  average  shear  over  these  developing  tropical  circulations  was  -1.8  ms*. 

If  it  is  assumed  that  the  hours  prior  to  It)  represent  the  time  during  which  the 
tropical  cyclone  seedling  undergoes  the  formation  process,  then  the  analyses  and 
forecasts  of  the  selected  variables  during  that  time  yield  information  about  the  formation 
environment.  Focusing  on  the  850- mb  relative  vorticity  and  deep  layer  wind  shear  in  the 
24-36  h  leading  up  to  the  first  warning  time,  it  is  apparent  that  the  over- forecasts  of 
vorticity  and  deep  layer  shear  are  consistent  with  one  another  (increased  vorticity 
forecasts  correspond  to  more  easterly  shear).  This  consistency  indicates  that  the 
NOGAPS  model  is  accurately  handling  the  physical  parameterizations. 

At  the  same  time,  NOGAPS  forecasts  of  mid- level  vapor  pressure,  which 
indicates  the  amount  of  moisture  available  to  the  developing  circulation,  were  not  only 
drier  than  analyzed,  but  also  had  a  substantially  smaller  signal  than  was  expected,  since  it 
varied  less  than  0.5  Pa  over  the  lifetime  of  the  analyzed  storms.  While  the  lack  of 
variability  of  this  variable  during  the  early  stages  of  development  and  intensification  was 
surprising,  the  vapor  pressure  outliers  were  almost  all  significantly  drier  than  the 
analyzed  value.  In  a  global  model  such  as  NOGAPS,  erroneous  vapor  pressure  forecasts 
are  expected  to  be  too  moist,  rather  than  too  dry.  This  excessive  drying  of  the  mid- levels 
within  the  NOGAPS  model  is  likely  a  function  of  the  convective  parameterization,  and 
works  against  the  formation  processes.  If,  as  in  this  case,  the  internal  dynamics  of  the 
tropical  seedling  vortex  inhibit  formation,  the  formation  process  shifts  from  a  mesoscale 
interaction  to  a  large-scale  process.  Since  the  NOGAPS  model  is  capable  of  accurately 
forecasting  tropical  cyclone  formation  with  inhibited  formation  on  the  mesoscale,  the 
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large-scale  processes  are  more  likely  accurately  represented  by  the  model.  However, 
improvements  could  be  made  in  the  modeling  of  formation  if  these  mesoscale  processes 
could  be  included. 

Sea- level  pressure  and  925- mb  wind  speed  had  less- well  developed  indicators  for 
developing  circulations,  and  thus  serve  more  as  proxies  for  the  intensity  of  the 
circulation,  rather  than  providing  an  indication  of  the  state  of  an  atmospheric  predictor 
that  either  enhances  or  inhibits  formation.  While  these  variables  taken  individually  may 
not  provide  obvious  indications  of  formation,  they  are  significant  to  the  formation 
forecast  process  in  that  they  represent  significant  features  within  the  large-scale 
atmosphere. 

Several  significant  differences  between  the  2002  and  2003  data  sets  were  noted. 
Fewer  tropical  cyclones  were  numbered  by  JTWC  during  2003,  during  which  a  weak  La 
Nina  flow  pattern  dominated  the  western  North  Pacific.  More  variation  existed  in  the 
analyzed  vorticity  curve  for  the  2003  developing  circulations  (Figure  3.37)  than  in  2002 
(Figure  3.5),  and  the  range  of  forecast  values  was  smaller.  The  average  analyzed 
vorticity  at  the  first  JTWC-waming  time  was  similar  (5.0  x  10'^  s'^  in  2002  versus  4.8  x 
10'^  s'^  in  2003),  although  the  change  in  vorticity  in  the  first  12  h  after  the  first  warning 
time  when  the  synthetic  tropical  cyclone  observations  were  inserted  into  the  NOGAPS 
model  analysis  was  larger  in  2003  than  during  2002,  which  indicates  that  the  model 
analysis  may  have  been  slightly  below  what  actually  occurred  in  the  environment. 

Similarly,  the  NOGAPS  analyses  of  deep  layer  wind  shear  were  very  similar  in 

2002  and  2003,  and  the  average  analyzed  value  at  the  first  JTWC  warning  time  was 
slightly  smaller  in  2003  (-1.6  m  s'^)  than  2002  (-1.8  m  s'^).  The  larger  (more  negative) 
value  in  2002  indicates  that  the  deep  layer  wind  shear  was  more  easterly  in  2002.  The 
weaker  analyzed  easterly  wind  shear  in  2003  may  be  attributed  to  the  weak  La  Nina  flow 
pattern  that  dominated  the  western  North  Pacific  during  that  year. 

It  is  possible  that  some  of  the  decrease  in  the  forecast  range  of  the  variables  in  the 

2003  data  set  is  associated  with  changes  made  to  the  NOGAPS  model  during  2003. 
These  changes  include  an  increase  in  the  amount  of  convective  momentum  transport  in 
the  Emanuel  cumulus  parameterization  scheme,  and  an  upgraded  data  assimilation 
system.  The  transition  from  a  multivariate  optimum  interpolation  (MVOI)  data 
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assimilation  scheme  to  the  Navy  Data  Assimilation  System  (NAVDAS),  which  has  the 
capability  to  ingest  satellite  observations  in  addition  to  traditional  surface,  radiosonde, 
and  synthetic  observations,  occurred  during  October  2003.  Changes  made  to  the  process 
by  which  synthetic  tropical  cyclone  observations  are  assimilated  into  the  NOGAPS 
model  analysis  should  be  transparent  to  the  user  (B.  Strahl,  FNMOC,  personal 
conversation),  but  yield  a  more  accurate  model  analysis.  However,  the  changes  to  the 
convective  momentum  transport  during  2003  may  have  yielded  a  decrease  in  the  easterly 
deep  layer  wind  shear,  which  decreased  the  number  of  false  alarm  circulations  developed 
by  the  NOGAPS  model. 

B.  NON -DEVELOPING  VORTICES 

The  non- developing  vortices  were  separated  into  three  main  categories:  those  that 

met  the  minimum  duration  criterion  of  24  h;  those  that  did  not  meet  the  minimum 
duration  criterion;  and  false  alarms.  The  non- developing  vortices  that  met  the  minimum 
duration  criterion  comprised  the  NDl  data  set.  The  Group  Two  subset  of  non- developing 
vortices  (NDG2)  was  extracted  from  the  NDl  data  set,  and  consisted  of  the  non¬ 
developing  vortices  that  met  the  minimum  duration  criterion,  and  also  \^ere  forecast  to 
exceed  the  850- mb  relative  vorticity  threshold  (5.0  x  10'^  s‘^)  established  for  developing 
vortices.  The  non- developing  vortices  that  did  not  meet  the  minimum  duration  criterion 
were  placed  in  the  ND2  data  set.  Non-developing  over- forecast  vortices  (NDOF)  were 
identified  from  the  ND2  data  set  as  those  that  were  forecast  to  exist  for  at  least  24  hours 
longer  than  they  actually  existed  in  the  NOGAPS  model  analyses.  After  the  NDOF 
vortices  were  extracted,  the  remaining  vortices  in  the  ND2  data  set  were  not  analyzed 
further  since  the  vortices  contained  therein  corresponded  to  short-term  convectively- 
driven  circulations  and  one-time  spurious  circulations.  The  third  main  data  set  of  model 
false  alarms  was  identified  as  vortices  that  were  forecast  to  develop  for  at  least  three 
consecutive  model  runs,  but  were  never  analyzed.  These  model  false  alarms  were  also 
compared  to  the  TCFA  false  alarms. 

Since  no  first  warning  time  existed  for  each  non- developing  vortex,  and  no  first 
Best- Track  time  had  been  established,  the  reference  time  for  each  non- developing  vortex 
(No)  that  met  the  minimum  duration  criterion  was  defined  to  be  the  time  of  the  maximum 
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analyzed  850- mb  relative  vorticity.  This  Nb  time  could  be  calculated  for  each  non¬ 
developing  vortex  h  the  NDl  and  NDG2  data  sets.  Since  the  vortices  in  the  ND2  and 
NDOF  data  sets  corresponded  to  only  one  or  two  analyses,  the  reference  time  used  for 
those  data  sets  (N*o)  corresponded  to  the  last  verifying  analysis. 

The  threshold  values  for  the  non-developing  vortices  relative  to  No  and  N*o  for 
each  of  the  five  selected  variables  (Table  4.2)  may  be  compared  to  the  threshold  values 
determined  for  the  developing  circulations  in  Table  4.1.  In  general,  the  threshold  values 
for  the  non- developing  vortices  were  smaller,  closer  to  zero,  or  less  favorable  for 
formation  than  the  threshold  values  for  the  developing  circulations. 


Table  4.2.  Summary  of  threshold  values  relative  to  No  and  N*o.  Values  correspond  to 
the  average  analyzed  value  of  the  variable  at  the  reference  time. 


Forecast  Variable 

Analyzed  Threshold  Values  (Rei 

ference  Value) 

NDl  (No) 

NDG2  (No) 

NDOF  (N*o) 

8 50- mb  relative  vorticity 

MgliiMM 

Deep  layer  (200-850  mb)  wind  shear 

4.1  ms'^ 

Sea- level  pressure 

1008.0  mb 

1007.0  mb 

1007.8  mb 

925-mb  wind  speed 

6.7  m  s'* 

Vapor  pressure  (500-700  mb  avg.) 

3.4  Pa 

3.6  Pa 

3.4  Pa 

The  850- mb  relative  vorticity  threshold  for  the  non- developing  vortices  (3.3  x  10“^ 
s‘^)  was  less  than  the  vorticity  for  the  developing  circulations  (5.0  x  10'^  s'^),  which 
indicates  that  less  organization  existed  within  the  non- developing  circulations  at  the  time 
of  maximum  analyzed  vorticity  than  in  the  developing  circulations  at  the  first  JTWC- 
waming  time.  The  analyzed  vorticity  for  the  developing  circulations  would  have  been 
equal  to  the  analyzed  value  for  these  non- developing  cases  (at  No)  roughly  96  h  prior  to 
Fq.  If  it  is  assumed  that  this  equivalence  of  vorticity  values  represents  a  similar 
developmental  stage,  this  signifies  the  maximum  intensity  of  the  developing  (non¬ 
developing)  vortices  do  (do  not)  undergo  further  development. 

The  average  analyzed  shear  at  No  (3.3  x  10'^  s'^)  was  also  less  than  the  minimum 
shear  value  (4  x  10'^  s'*)  at  which  deep  layer  wind  shear  transitioned  from  westerly  to 
easterly.  That  the  deep  shear  on  average  was  less  than  this  threshold  value  indicates  that 
wind  shear  remained  westerly  throughout  the  lifecycle  of  the  non- developing  vortex,  and 
never  transitioned  to  the  more  formation-  favorable  easterly  shear. 
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One  potential  explanation  for  the  lack  of  further  development  of  the  non¬ 
developing  vortices  is  considered  to  be  the  sign  of  the  analyzed  deep  layer  wind  shear. 
Since  the  analyzed  deep  layer  shear  remains  positive  (westerly)  throughout  the  lifetime  of 
the  non- developing  vortices,  the  atmosphere  in  which  these  circulations  exist  is 
dominated  by  westerly  shear,  which  is  not  favorable  for  tropical  cyclone  formation. 

As  was  the  case  for  the  developing  vortices,  neither  the  sea- level  pressure  nor  the 
vapor  pressure  variable  had  a  strong  signal  for  forecasting  the  development  potential  of 
these  tropical  storm  seedlings.  In  contrast  to  the  vapor  pressure  for  the  developing 
circulations,  the  vapor  pressure  for  the  non- developing  cases  did  have  more  variation 
relative  to  the  average  analyzed  value  at  Nq.  While  the  occurrence  time  of  the  analyzed 
maximum  vapor  pressure  did  not  correspond  to  the  occurrence  time  of  the  analyzed 
maximum  vorticity,  the  vapor  pressure  for  the  non- developing  cases  was  consistently  less 
than  the  threshold  value  (3.6  Pa)  determined  for  the  developing  circulations,  which 
indicates  both  a  weaker  potential  warm  core  and  that  the  formation  environment  may  be 
less  favorable  for  further  development  of  the  tropical  storm  seedling. 

Additionally,  the  outliers  of  vapor  pressure  for  the  non- developing  circulations 
were  predominantly  less  than  the  analyzed  vapor  pressure  for  this  data  set,  which 
indicates  forecasts  for  a  drier  mid- troposphere  than  was  analyzed.  Forecasts  for 
excessive  dryness  in  the  mid- levels  are  counter-intuitive,  since  it  is  expected,  in  general, 
that  the  global  models  will  forecast  excessive  moisture  as  a  result  of  convective 
parameterization  of  meso-  and  smaller- scale  processes.  As  noted  for  the  developing 
circulations,  this  excessive  dryness  indicates  that  the  internal  dynamics  of  the  tropical 
vortex  are  working  against  formation.  While  accurate  formation  forecasts  are  possible 
with  this  counter-balance,  it  is  also  assumed  that  this  helps  to  prevent  more  false  alarms 
and  over- forecasts  than  those  already  observed. 

The  analyzed  sea- level  pressure  for  the  non- developing  vortices  was  at  all  times 
higher  than  the  average  analyzed  sea- level  pressure  for  the  developing  vortices.  While 
the  sea- level  pressure  in  the  early  stages  of  development  for  both  the  developing  and  non¬ 
developing  circulations  was  similar,  the  sea- level  pressure  for  the  developing  vortices 
began  to  distinctly  decrease  about  48  h  prior  to  Fq,  and  the  sea-level  pressure  for  the  non¬ 
developing  vortices  never  showed  a  corresponding  decrease.  Rather,  the  sea- level 
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pressure  for  the  non- developing  vortices  remained  consistently  within  about  1  mb  of  the 
average  analyzed  value  at  Nq.  Since  the  non- developing  vortices  never  experienced  a 
substantial  decrease  in  sea- level  pressure,  the  magnitude  of  the  warm  core  of  the  non¬ 
developing  vortices  was  presumably  smaller  than  for  the  developing  circulations. 

Although  the  925- mb  wind  speed  variable  did  not  provide  a  strong  formation 
signal,  it  nonetheless  is  a  proxy  for  the  intensity  of  the  lo\w  level  circulation.  This 
analyzed  wind  speed  fluctuates  throughout  the  lifetime  to  indicate  pulses  of  development 
and  decay  that  occurred.  The  average  analyzed  wind  speed  for  non- developing  vortices 
was  at  all  times  less  than  the  threshold  value  determined  for  the  developing  circulations 
(10.8  m  s^).  Even  at  120  h  prior  to  It),  the  925- mb  wind  speed  of  the  developing 
circulations  was  greater  than  any  of  the  reliable  analyzed  wind  speed  values  for  the  non¬ 
developing  vortices.  Thus,  the  non- developing  vortex  never  achieved  the  same  intensity 
as  the  developing  circulations,  and  the  low-level  circulation  never  achieved  the 
magnitude  capable  of  further  development  from  air- sea  fluxes. 

1.  False  Alarms 

False  alarms  in  the  NOGAPS  model  are  a  result  of  the  model  tendency  to  over¬ 
forecast  the  development  of  tropical  circulations.  Two  types  of  false  alarms  were 
identified:  model  false  alarms,  which  included  non- developing  vortices  that  were  forecast 
to  occur  for  at  least  three  consecutive  forecasts,  but  were  never  analyzed;  and  the  so- 
called  Group  2  false  alarms  that  included  analyzed  non- developing  vortices  that  were 
forecast  to  exceed  the  relative  vorticity  threshold  determined  for  developing  circulations. 

The  forecasts  for  the  model  false  alarms  in  which  a  vortex  was  never  analyzed, 
approximately  resemble  the  early- stage  forecasts  of  developing  circulations,  with  the 
forecast  850- mb  relative  vorticity  forecast  to  exceed  the  analyzed  values.  Examination  of 
the  selected  forecast  variables  would  likely  reveal  excessive  850-mb  relative  vorticity, 
and  more  easterly  deep  layer  wind  shear.  It  may  be  that  the  forecasts  for  the  model  false 
alarms  are  similar  to  those  for  the  non- developing  circulations  that  were  forecast  to 
exceed  the  850-mb  relative  vorticity  threshold  (NDG2)  determined  for  developing 
vortices  (5.0  x  10'^  s‘^).  A  more  easterly  wind  shear  is  also  expected  in  conjunction  with 
this  vorticity  over- forecast.  This  tendency  cannot  be  confirmed  here  as  analyses  for  the 
selected  forecast  variables  are  not  available  since  these  vortices  never  existed  in  an 
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analysis  panel.  However,  incorrect  forecasts  of  these  two  variables  yields  excessive 
development  of  what  might  otherwise  be  a  short-lived  convective  event. 

The  15  model  false  alarms  that  occurred  during  1  May  -  31  October  2002 
represent  about  13%  of  the  total  number  of  tracked  circulations.  Thus,  more  than  one 
model  false  alarm  was  detected  for  every  two  developing  circulations  (24),  which 
indicates  a  problem  for  the  forecaster  in  interpreting  NOGAPS  formation  forecasts. 
During  2003,  the  number  of  model  false  alarms  jumped  to  27  (column  four  of  Table 
3.10). 

The  Group  2  false  alarms  in  which  a  non- developing  vortex  was  forecast  to 
exceed  the  850- mb  relative  vorticity  threshold  determined  for  developing  circulations 
(5.0  X  10'^  s'^)  represent  the  second  type  of  false  alarms  in  the  NOGAPS  model.  These 
false  alarms  pose  a  serious  challenge  for  operational  forecasters  because  each  vortex  must 
be  examined  since  it  exists  in  the  NOGAPS  analysis.  Examination  of  the  850- mb  relative 
vorticity  and  deep  layer  wind  shear  profiles  indicate  that  the  early  stages  of  development 
of  these  vortices  are  very  similar  to  similar  stages  in  the  lifecycle  of  the  developing 
circulations.  Short-range  vorticity  forecasts  prior  to  the  time  of  the  maximum  analyzed 
vorticity  are  very  similar  to  the  analyzed  vorticity  values,  which  indicates  that  the 
NOGAPS  model  is  accurately  handling  both  the  physical  parameterizations  of  the  vortex 
and  the  large-scale  environmental  flow.  The  mechanism  that  leads  to  the  sudden 
decrease  in  the  analyzed  850- mb  relative  vorticity  is  unclear. 

Though  the  number  of  NOGAPS  model  false  alarms  in  2003  was  greater  than  in 

2002,  there  was  a  significant  decrease  in  the  number  of  Group  2- type  false  alarms.  Since 
fewer  than  five  Group  2  false  alarms  were  forecast  in  2003,  these  were  not  analyzed  due 
to  small  sample  size.  The  decrease  in  the  number  of  Group  2  false  alarms  may  be 
attributed  to  the  increase  in  convective  momentum  transport  of  the  Emanuel  cumulus 
parameterization  scheme  (Hogan  et  al.  2004).  The  ultimate  result  of  increasing  the 
amount  of  convective  momentum  transport  was  to  decrease  the  easterly  deep  layer  wind 
shear  that  was  forecast  by  the  NOGAPS  model,  which,  in  turn,  significantly  decreased 
the  number  of  false  alarm  forecasts  that  were  generated  by  NOGAPS  during  2003. 

Despite  the  significant  decrease  in  the  number  of  Group  2  false  alarms  during 

2003,  the  number  of  model  false  alarms  (27)  actually  increased  from  2002  (15).  This 
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increase  is  attributed  to  year-to-year  environmental  variability  in  the  western  North 
Pacific  basin.  A  detailed  explanation  of  this  increase  required  a  more  extensive 
investigation  of  the  analysis  and  forecast  fields  associated  with  these  vortices.  The 
tentative  explanation  is  that  the  overall  environmental  conditions  were  less  favorable  for 
tropical  cyclone  formation  during  2003.  However,  the  NOGAPS  model  continued  to 
over- forecast  the  spin- up  of  vortices  for  at  least  three  consecutive  forecasts  that  simply 
did  not  verify  in  this  less- favorable  environment  during  2003.  Thus,  the  number  of 
model  false  alarms  can  be  larger  even  though  the  over- forecast  vorticity  tendency  of 
NOGAPS  with  the  generally  weaker  initial  vortices  during  2003  did  not  lead  to  forecasts 
exceeding  5.0  x  10'^  s'\  which  is  the  definition  of  Group  2  false  alarms  (markedly 
decreased  during  2003). 

2.  Model  Over-forecasts 

Model  false  alarms  were  distinguished  from  over- forecast  non- developing 
vortices  in  that  the  model  false  alarms  never  existed  in  a  NOGAPS  analysis,  whereas  the 
over- forecast  non- developing  vortices  existed  in  one  or  two  NOGAPS  analyses.  The 
existence  of  over- forecast  non- developing  vortices  is  another  indicator  of  the  NOGAPS 
tendency  to  over-develop  tropical  circulations.  The  reasons  for  this  over- development 
are  likely  linked  to  the  NOGAPS  model  parameterization  schemes  that  define  and 
develop  the  internal  vortex  of  the  system.  That  a  circulation  can  be  identified  in  the 
NOGAPS  model  analyses  indicates  that  a  potential  tropical  cyclone  seedling  is  present  in 
the  environment.  However,  the  actual  and  forecast  development  profiles  of  these  vortices 
differ. 

In  general,  these  over- forecast  non- developing  vortices  existed  in  two  consecutive 
analyses,  and  had  a  relatively  low  analyzed  850- mb  relative  vorticity  (2.4  x  10'^  s‘^) 
compared  to  the  vorticity  thresholds  determined  for  the  non- developing  (3.3  x  10'^  s‘^) 
and  developing  circulations  (5.0  x  10"^  s'^).  Solely  based  on  the  average  analyzed 
vorticity  at  the  reference  time  (the  last  analysis  time.  No  and  Fo,  respectively),  the 
analyzed  vorticity  for  this  subset  of  non- developers  is  lower  than  the  minimum  (4  x  10“^ 
s'*)  identified  previously  that  is  necessary  in  order  for  deep  layer  wind  shear  to  transition 
from  the  inhibitive  westerly  to  the  more  favorable  easterly. 
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However,  the  forecast  850- mb  relative  vorticity  always  exceeded  the  analyzed 
vorticity,  and  was  generally  forecast  to  continue  to  increase  in  time.  The  erroneous 
vorticity  forecasts  produced  by  the  NOGAPS  model  for  these  non- developing  vortices 
had  an  associated  erroneous  easterly  deep  layer  wind  shear.  Furthermore,  the  sea- level 
pressure  forecasts  indicated  a  more  favorable  environment  for  formation  than  was 
analyzed.  Why  these  vortices  were  so  significantly  over- forecast  is  unclear.  Since  these 
vortices  formed  preferentially  in  association  with  the  monsoon  trough  (in  the  South  China 
Sea  and  Philippine  Sea),  it  may  indicate  that  they  may  have  been  short-lived  convective 
events,  or  weak  circulations  associated  with  the  monsoon  trough  or  depression. 
Inaccurate  representation  of  either  the  seedling  disturbance  or  the  formation  environment 
may  have  yielded  incorrect  forecasts  of  development  by  the  NOGAPS  model. 

Since  the  number  of  over- forecast  non- developing  vortices  decreased  dramatically 
during  2003,  the  modifications  made  to  enhance  the  convective  momentum  transport  of 
the  Emanuel  cumulus  parameterization  scheme  during  late  2002  were  effective. 
Increasing  the  amount  of  convective  momentum  transport  resulted  in  a  decrease  in  the 
amount  of  easterly  wind  shear,  and  evidently  made  the  forecast  environment  less 
favorable  for  formation.  Weaker  easterly  shear  decreased  the  number  of  spurious 
convective  events,  and  thus  the  number  of  non- developing  vortices  that  were  forecast  to 
develop  to  tropical  depression  strength  or  greater. 
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V.  RECOMMENDATIONS  EOR  EURTHER  STUDY 


The  use  of  the  TCVTP  in  assessing  model  performance  in  the  prediction  of 
tropical  cyclone  formation  should  be  continued  in  other  studies.  Specifically,  subsequent 
studies  in  the  following  list  would  be  beneficial  to  tropical  cyclone  formation  prediction 
assessment  while  also  documenting  the  real-time  potential  of  the  TCVTP: 

•  Further  analysis  of  the  selected  variables  based  on  a  multiple  year  sample 
of  circulations  in  the  western  North  Pacific. 

o  Identify  predictors  for  objective  assessment  of  potential  tropical 
cyclone  formation. 

o  Stratification  of  the  larger  data  set  by  examining  circulations  that 
formed  in  the  South  China  Sea,  Philippine  Sea,  and  east  of  the 
monsoon  trough,  and  in  addition  to  separating  the  Philippine  Sea 
into  the  monsoon  trough  and  a  new  subregion  north  of  the 
monsoon  trough. 

o  Further  analysis  of  the  larger  sample  by  examining  circulation 
formation  time  relative  to  the  TCFA  formation  window. 

•  Further  study  that  applies  the  TCVTP  in  the  western  North  Pacific  to  other 
global  numerical  models  for  comparison  to  NOGAPS. 

o  Comparative  studies  of  formation  thresholds  for  circulations 
forming  in  different  tropical  cyclone  basins. 

o  Application  and  use  of  the  TCVTP  in  real-time  operations. 

•  Further  analysis  of  model  false  alarm  vortices  and  non- developing  over¬ 
forecast  vortices,  including  assessment  of  the  large-scale  environment  and 
vortex  characteristics  throughout  the  forecast  lifetime  of  the  circulation. 

•  Further  study  using  the  first  Best- Track  time  as  the  formation  time,  as  well 
as  the  first  time  the  TCVTP  captured  the  circulation  in  the  analysis. 
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APPENDIX  A:  TCFA  CRITERIA  CHECKLIST  -  WESTERN 
NORTH  PACIFIC  &  NORTH  INDIAN  OCEANS 

The  JTWC  uses  the  following  checklist  to  determine  whether  or  not  a  Tropical 
Cyclone  Formation  Alert  (TCFA)  should  be  issued  for  a  suspect  circulation  center.  Items 
denoted  with  an  asterisk  are  progressive  in  nature.  For  example,  if  any  wind  speed  within 
5  degrees  of  the  circulation  center,  as  determined  by  satellite  imagery  and/or  hand 
analysis,  is  35  knots,  then  a  total  of  6  points  will  be  contributed  by  items  E,  F,  and  G.  A 
TCFA  is  issued  when  the  sum  of  points  is  30  or  greater. 


Point  Value 


2 

4 

4 

3 


1 

2 

3 

1 

3 

1 

2 

3 


1 

2 


1 

3 

1 


Item: 

I.  SFC/  Gradient  Level 

A.  A  circulation  is  evident  in  the  wind  field 

B.  A  circulation  has  been  evident  for  24  hours 

C.  Environmental  MSLP  -  CNTR  SEP  =  4MB  (EST) 

D. *  Westerly  SEC/Gradient  level  winds  of  at  least  10  kt 

south  of  the  circulation  center 

n.  Disturbance,  and  within  05  degrees  of  center 

E. *  Any  wind  associated  with  center  is  at  least  20  kt 
E.*  Any  wind  associated  with  center  is  at  least  25  kt 

G. *  Any  wind  associated  with  center  is  at  least  30  kt 

H. *  24  hour  pressure  decrease  at  nearby  station  =  2  MB 
L*  24  hour  pressure  decrease  at  nearby  station  =  3  MB 

J. *  EST.  MSLP  of  TD  is  ^1008  MB 

K. *  EST.  MSLP  of  TD  is  ^1006  MB 

L. *  EST.  MSLP  of  TD  is  ^1004  MB 

III.  500  mb 

A.  There  is  evidence  of  at  least  a  trough 

B.  There  is  evidence  of  a  closed  circulation 

rV.  200  mb 

A.  TUTT  to  the  northwest  of  the  TD 

B.  Evidence  of  an  anticyclone  over  the  center  of  the  TD 

C.  200  mb  wind  in  the  center  <  25  kt 
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V. 


SST 


A.  SST>28°C 
VI.  Satellite  Data: 

A. *  The  TD  has  persisted  for  at  least  24  hours 

B. *  The  TD  has  persisted  for  at  least  48  hours 

C. *  The  TD  has  persisted  for  at  least  72  hours 

D. *  Dvorak  classification  of  at  least  TO.O 

E. *  Dvorak  classification  of  at  least  Tl.O 

F. *  Dvorak  classification  of  at  least  T2.0  (warning 

should  already  be  issued) 

Vn.  Miscellaneous: 

A.  Double  vortex  interaction  (cross- equatorial)  exists 

B.  Tropical  disturbance  is  within  72  hours  of  a  DoD 

resource 

C.  Synoptic  circulation  and  satellite  fix  are  consistent  in 

location  (within  02  degrees) 

D.  20  kt  synoptic  wind  reports  within  3  degrees  of  the 

satellite  fix  (does  not  apply  to  winter  gales) 


Total  Points  Possible 
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APPENDIX  B:  GENERATION  OE  BOX  PLOTS 


Box  plots  for  each  selected  variable  indicate  the  range  of  specified  forecast 
variable  values  at  each  forecast  time  relative  to  the  Fb,  F*o,  or  No,  which  are  the  first 
JTWC  warning  time  for  numbered  storms,  the  first  Best- Track  analysis  time  for 
numbered  storms,  and  the  time  corresponding  to  the  maximum  analyzed  850- mb  relative 
vorticity  for  non- developing  vortices,  respectively.  Each  figure  consists  of  four  panels, 
which  represent  the  forecast  intervals  of  24-h,  48-h,  72-h,  and  96-h  (see  Figure  B.l). 
Although  the  features  within  each  figure  are  similar,  the  reference  time  used  to  calculate 
the  values  varies  depending  on  the  data  set  analyzed.  The  title  on  each  panel  indicates 
the  reference  time  (e.g.,  Fo),  the  year  and  basin  in  which  the  vortex  occurred  (e.g.,  2002 
WPAC),  the  numerical  forecast  model  (e.g.,  NOGAPS),  the  data  set  to  which  the  vortex 
belonged  (e.g..  Developing  DV  or  Non- developing  NDl),  and  the  forecast  time  (e.g., 
-1-24).  The  abbreviations  used  to  identify  the  data  set  used  are:  DV,  NDl,  NDG2,  ND2, 
NDOF,  where  DV  represents  those  vortices  that  developed  into  storms  that  were  later 
numbered  by  JTWC,  NDl  represents  the  non- developing  vortices  that  met  the  minimum 
duration  criterion,  NDG2  represents  the  non- developing  vortices  that  met  the  minimum 
duration  criterion  and  were  forecast  to  exceed  the  ?o  threshold,  and  ND2  represents  the 
non- developing  vortices  that  did  not  meet  the  minimum  duration  criterion.  NDOF 
represents  the  non- developing  vortices  that  did  not  meet  the  minimum  duration  criterion, 
but  did  meet  the  over-  forecast  criteria. 

Separate  charts  for  each  formation  region  (SCS,  PS,  and  EMT)  were  also  created 
for  the  850- mb  vorticity  variable  for  the  DV,  and  NDl  data  set.  The  thin  horizontal  lines 
drawn  across  each  regional  chart  indicate  the  mean  forecast  variable  value  for  that  region. 
For  example,  the  mean  vorticity  value  at  Fq  for  all  developing  vortices  is  5.0  x  10'^  s‘\ 
while  the  value  for  those  vortices  developing  in  the  SCS  is  4.8  x  10'^  s'\  5.5  x  10'^  s'*  in 
the  PS  region,  and  4.7  x  10'^  s'*  in  the  EMT  region. 

A.  ANALYSIS  DATA 

The  heavy  black  line  in  each  panel  represents  the  average  analyzed  value  (at  -i-OO) 
of  the  indicated  variable,  and  thus  represents  a  threshold  value  for  either  developers  or  for 
non- developers  as  appropriate.  The  thin  green  whiskers  that  extend  from  the  average 
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Figure  B.l.  Average  analyzed  and  forecast  850- mb  relative  vorticity  for  developing 
vortices  relative  to  the  first  JTWC  warning  time  (Fq).  Heavy  black  line  represents 
analyzed  average  850- mb  relative  vorticity  (at  -i-OO),  and  red  line  represents  average 
forecast  850- mb  relative  vorticity  for  the  indicated  forecast  time.  Average  850- mb 
relative  vorticity  at  Fo  (light  solid  black  line)  was  5.0  x  10'^  s'\  which  is  then  a  reference 
for  the  vorticity  magnitudes  before  (to  the  left)  and  after  (to  the  right)  the  first  warning 
time.  The  sample  sizes  available  for  comparison  of  the  forecast  values  at  each  forecast 
verification  time  are  listed  near  the  top  of  each  panel  for  24-h  (upper- left),  48-h  (upper- 
right),  72-h  (lower- left),  and  96-h  (lower-right)  forecasts  of  these  developing  storms. 


analyzed  value  represent  the  standard  deviation  (i.e..  Table  B.l)  of  the  selected  variable 
at  each  time  interval.  The  thin  black  horizontal  line  extending  from  -120  h  to  -1-120  h 
cross  each  plot  is  the  mean  value  of  the  selected  variable  for  the  developing  vortices  at 
the  -1-00  analysis  time. 
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B.  FORECAST  DATA 

Separate  box  plots  were  created  for  each  forecast  verification  time  relative  to  Fo, 
F*o,  and  Nb  from  +00  to  +120  h  if  the  forecast  can  be  verified.  For  example,  the  24-h 
forecast  in  the  upper- left  panel  of  Figure  B.l  represents  all  of  the  24-h  850  mb  relative 
vorticity  forecasts  for  a  sample  of  developing  vortices  that  verified  at  various  times 
relative  to  Fq.  Similarly,  the  upper- right  panel  in  Figure  B.l  represents  all  48-h  850- mb 
relative  vorticity  forecasts  that  verified  at  various  times  relative  to  Fq.  The  heavy  red 
line  in  each  panel  represents  the  average  forecast  value  at  the  time  indicated  at  the  end  of 
the  title  for  each  panel  (i.e.,  +24,  +48,  +72,  or  +96).  The  total  number  of  forecasts 
included  in  each  column  for  each  forecast  time  relative  to  Fo,  F*o,  or  No  are  indicated  at 
the  top  of  the  column.  The  range  of  forecast  values  for  each  variable  that  can  be  verified 
at  that  time  relative  to  Fq  is  indicated  at  each  time  interval  by  the  blue  box  plots.  The  top, 
middle,  and  bottom  horizontal  lines  that  make  up  each  box  represent  the  upper  quartile, 
median,  and  lower  quartile.  The  dashed  whiskers  extending  from  the  box  indicate  the 
values  that  are  not  within  the  lower  and  upper  quartiles,  and  the  outliers  are  indicated  by 
an  open  circle.  The  dashed  horizontal  line  in  the  lower  portion  of  the  plot  is  the  mean 


Table  B.l.  Summary  of  standard  deviations  of  850- mb  relative  vorticity  (x  10'^  s‘^)  for 
each  forecast  time  for  +96  hour  forecasts  of  developing  vortices. 


Forecast  Time 

Std  Dev 

Forecast  Time 

Std  Dev 

-120 

0.71 

+00 

2.21 

-96 

1.54 

+12 

2.00 

-84 

1.34 

+24 

2.08 

-72 

1.81 

+36 

2.63 

-60 

1.65 

+48 

3.00 

-48 

1.84 

+60 

2.93 

-36 

1.58 

+72 

3.33 

-24 

1.92 

+84 

3.25 

-12 

1.59 

+96 

3.42 

— 

— 

+120 

3.69 

value  of  the  selected  variable  at  the  F*o  or  No  times,  or  the  region- specific  average  value, 
as  indicated.  Data  from  non- developing  vortices  were  not  included  in  the  plots  of  Fq  and 
F*o,  which  refer  only  to  the  developing  vortices. 
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